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is shown in black 
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5.13 Ratio of the number of D'^ to as a function of reconstructed D 
momentum at i?cm = 4170 MeV. The differences in the branching 
fractions were included in the ratio, while the efficiencies were take 
to be constant over the momentum range. The MC includes only 
D*Dtt multi-body production and uses a simple ISR model that 
assumes all events are produced from the '?/'(4160) resonance, with 
a width of 78 MeV. The large discrepancy between the data and the 
MC between 200 and 400 MeV/c clearly indicates that something 
is incorrect in our assumptions used to generate the MC sample. . 

5.14 Ratio of D*^ to D*^, both using the 7r° channel, with the subse- 
quent D-meson decays K^tt^ and K^tt^ti^, respec- 
tively. The ratio includes the branching fractions but excludes the 
efficiencies, which are constant over the momentum range. The 
MC includes only D*D7c multi-body production and uses a simple 
ISR model that assumes all events are produced from the ?/'(4160) 
resonance, with a width of 78 MeV. The discrepancy at 650 MeV/c 
indicates that ISR is important and that it is incorrectly handled 
in our MC sample 
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5.15 The two-body cross sections used in EVTGEN for DD^ D*D, and 
D*D* for the MC samples with improved treatment of ISR. A sim- 
ple implementation of the cross sections were made using straight- 
line approximation for the real, Born-level, two-body, cross sec- 
tions. The main assumption is that the shape of the cross section 
will be preserved after applying radiative corrections, that is the 
observed cross section shape is the same as the Born-level cross 



section shape Il23 

5.16 The two-body cross sections used in EVTGEN for DgDs, D*Ds, and 
D*D* for the MC samples with improved treatment of ISR. A sim- 
ple implementation of the cross sections were made using straight- 
line approximation for the real. Born-level, two-body cross sections. 
The main assumption is that the shape of the cross section will be 
preserved after applying radiative corrections, that is the observed 
cross section shape is the same as the Born-level cross section shape. 1 124 
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5.17 Sideband-subtracted momentum spectrum for D'^ K^tx^ at 
4170 MeV. Data are shown as points with error bars which are 
fit to the improved ISR MC (histograms). The fit uses a spin- 
averaged phase-space model MC representation of D*D7t, shown 
in dark red. The sohd histograms correspond to i?^ < 1 MeV 
(no ISR) while the dashed histograms correspond to -E-^ > 1 MeV 
(with ISR). The total fit result is shown by the sohd black line 
{X^/NDF = 248/132) 

5.18 Sideband-subtracted momentum spectrum for — >■ K^-n-^n^ at 
4170 MeV. Data are shown as points with error bars which are 
fit to the improved ISR MC, the histograms. The fit uses a spin- 
averaged phase-space model MC representation of D*D7r, shown 
in dark red. The solid histograms correspond to i?^ < 1 MeV 
(no ISR) while the dashed histograms correspond to -E^ > 1 MeV 
(with ISR). The total fit result is shown by the solid black line 
ixVNDF = 182/132) 
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5.19 Sideband-subtracted momentum spectrum for 07r+ at 4170 MeV. 

Data are shown as points with error bars which are fit to the im- 
proved ISR MC (histograms). The sohd histograms correspond to 
< 1 MeV (no ISR) while the dashed histograms correspond to 
> 1 MeV (with ISR). The total fit result is shown by the solid 



black line {x^/NDF = 165/124) 

5.20 The ratio of the cross section results from the K'tt~^ fit to 
the results from the —>■ K^n^-ri^ fit. Since the fits are inde- 
pendent, the agreement in the ratios demonstrates the reliability 
of the fitting procedure in extracting the mode-by-mode cross sec- 
tions. The ratio, in all cases, should be one, indicated by the solid 
black line. All event types are in good agreement with what is 
expected 

5.21 Comparison between MC and data at 4170 MeV for D° i^'"7r+. 
Based on the momentum spectra fit results, an indepedent MC 
sample with the same statistics as the 178 pb^^ collected at 4170 
MeV was generated using all the knowledge gained about ISR, 
angular distributions (see the App.) and cross sections 
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5.22 Comparison between MC and data at 4170 MeV for L>+ R-ti+ti^ 
Based on the momentum spectrum fit results an independent MC 
sample with the same statistics as the 178 pb^l collected at 4170 
MeV was generated using all the knowledge gained about ISR, an- 
gular distributions (see the App.) and cross sections 

5.23 Comparison between MC and data, at 4170 MeV, for — > ^vr"*". 
Based on the momentum spectrum fit results an independent MC 
sample with the same statistics as the 178 pb~l collected at 4170 
MeV was generated using all the knowledge gained about ISR, 
angular distributions (see the App.) and cross sections 

5.24 Ratio of D*^ to D*+ in 4170 MeV data and MC as a function of 
the reconstructed D* momentum. Both D* states are detected 
with the 7r° channel, with the subsequent decays K~7i~^ and 

—>■ K^n^n^, respectively. The ratio is corrected for branching 
fractions, while efficiencies were assumed to be constant. The MC 
includes D*D-n multi-body and uses an updated ISR model that 
incorporates individual two-body cross sections as described in the 
text. The updated MC sample is in great agreement with the data 
at 4170 MeV 
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5.25 Ratio of the number of to as a function of reconstructed D 
momentum at E^ra = 4170 MeV. The ratio is corrected for branch- 
ing fractions, while efficiencies were assumed to be constant. The 
MC includes D*Dn multi-body and uses an updated ISR model 
that incorporates individual two-body cross sections as described 
in the text. The updated MC sample is in great agreement with 
the data at 4170 MeV 

5.26 Observed cross sections for e~^e~ DD, D*D and D*D* as a 
function of center-of-mass energy. They are determined by fitting 
the momentum spectrum for K^ti^ and K^ti^h^ 
with the updated ISR MC. Systematic errors are not included. . . 

5.27 Observed cross sections for e~^e~ D^Ds, D*Ds and D*D* as a 
function of center-of-mass energy. They are determined by fitting 
the momentum spectrum for (pir^ with the updated ISR 
model MC. A branching fraction of 3.6% was also used in the de- 
termination of the cross section. Systematic errors are not included. 
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5.28 Observed cross sections for e~^e~ D*D7t and D*D*tx as a func- 
tion of center-of-mass energy. They are determined by fitting the 
momentum spectrum for —>■ K^tt~^ and K^tt^tt^ with 
the updated ISR model MC. The multi-body contribution was 
modeled in MC using spin-averaged phase space. Systematic errors 
are not included 

5.29 Plot of the total charm cross section as calculated by each of the 
three methods described in the text, in addition to the results 
from the momentum-spectrum fit. The reason for the discrepancy 
between the inclusive and exclusive method is due to the presence 
of multi-body background and ISR. These two effects are taken into 
account when a fit to the momentum spectrum is performed. There 
is good agreement between all three methods, the inclusive, the 
hadron counting and the momentum fits. The previous exclusive 
result (ignoring multi-body) is shown, in addition, for comparison. 

5.30 Momentum-spectrum fit for the improved 178 pb^^ MC sample for 

K^n^. The MC sample was treated just as the data and 
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used as a cross-check of the momentum-spectrum fit procedure. . Il39 
5.31 Momentum-spectrum fit for the improved 178 pb~^ MC sample for 
K^n^Ti^. The MC sample was treated just as the data and 
used as a cross-check of the momentum-spectrum fit procedure. . Il40 
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5.32 Momentum-spectrum fit for tlie improved 178 pb^^ MC sample for 

07r+. The MC sample was treated just as the data and used 
as a cross-check of the momentum-spectrum fit procedure 

5.33 Missing-mass spectrum using K^n^ at 4170 MeV. The fit 
function consists of a Gaussian and an Argus function 

5.34 Missing-mass spectrum using D'^ — > K^tt^ti^ at 4170 MeV. The 
fit function consists of a Gaussian and an Argus function 

5.35 Missing-mass spectrum using —>■ K~7r~^, —>■ K^tt^tt^, and 

— >■ K^TT^n^TT^ at 4260 MeV. The fit function consists of two 
Gaussians, one for each peak, and an Argus function 

6.1 The fit result after adjusting the D*D cross section from its nom- 
inal value to one that is fiat between 4060-4170 MeV. The change 
results in a poor fit most notable at 600 MeV/c 

6.2 The fit result after adjusting the angle between the D in the rest 
frame of the D* with respect to the momentum of D* in the lab 
from its nominal value of 0.8 to 2.0. The fit is poor most notably 
between 500-600 MeV/c 
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7.1 The radiatively corrected inclusive cross sections using all methods 
are compared to the BES and CB -urfs-subtracted R data. Note 
that the corrected cross sections (including BES and CB) include 
the effect of vacuum polarization as defined in Eq. 17.41 as well as 
the bubbles due to n and r leptons |165 

7.2 The radiatively corrected inclusive cross sections using the KF in- 
terpolation method on the CB data. The error bars include the 
systematics as discussed in the text 1166 

7.3 R measurement obtained from these inclusive charm cross sections, 
radiatively corrected, using the KF interpolation method on the 
CB data. The error bars include the systematics as discussed in 
the text 

8.1 The production cross section times branching ratio for K~it~^ 
and K^-k^-k^ as a function of energy for this analysis, as 
compared to previous measurements. The results from the CLEO-c 
scan agree quite nicely with previous experiments 

8.2 The exclusive cross sections for e'^e^ D^D*^ (top) and e+e~ 
D*+D*- from Belle [6Tj. The dashed lines correspond to the ?/'(4040) 
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?/'(4160), Y(4260), and ^/'(4415). 
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8.3 The exclusive cross sections for e~^e~ — DD* (top) and e^e^ — >■ 
D*D*, after applying radiative corrections, for the work presented 
in this dissertation. The dashed lines correspond to the '?/'(4040), 
■^(4160), Y(4260), and tp{AA15). There is a factor of two difference 
between these results and those of Belle in Fig. 18.21 

8.4 Results for the non-strange D-meson cross sections. The error bars 
include both statistical and systematic errors 

8.5 Results for the Dg-meson cross sections. The error bars include 
both statistical and systematic errors 

8.6 The fit results from Ref. [66] which used preliminary CLEO-c 
results. The data points that are excluded in the fit are shown 
with circles {DD at 4010 and 4015 MeV and D*D* at 4015 MeV). 
The fit are to dimensionless rate coefficients Rk where Rk oc (Xfc 

(top left) corresponds to DD, R2 (top right) corresponds to 
DgDs, i?3 (bottom left) corresponds to DD*, and R4 (bottom right) 
corresponds to D*D*) 

8.7 Fit results, assuming a single resonance [6^, using the cross sec- 
tions presented in this dissertation. The fit are to dimensionless 
rate coefficients Rk where Rk oc cr^ (i?i (top left) corresponds to 
DD, R2 (top right) corresponds to DgDg, R^ (bottom left) corre- 
sponds to DD*, and R4 (bottom right) corresponds to D*D*). 
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.8 Comparison of the total charm cross section using the 3 different 

methods. This is a graphical representation of Table 18. 4[ 1X85 

.9 Results for the multi-body cross sections. The error bars include 

both statistical and systematic errors Il86 

.10 Sideband-subtracted momentum spectrum for K~tt~^ at 
3970 MeV. The data is the points with error bars and the his- 
tograms are MC Il99 

.11 Sideband-subtracted momentum spectrum for K~7t~^7[~^ at 
3970 MeV. The data is the points with error bars and the his- 
tograms are MC 1200 

.12 Sideband-subtracted momentum spectrum for — ^ K^n^ at 
3990 MeV. The data is the points with error bars and the his- 
tograms are MC 1201 

.13 Sideband-subtracted momentum spectrum for K^tt^tt^ at 
3990 MeV. The data is the points with error bars and the his- 
tograms are MC 1202 

.14 Sideband-subtracted momentum spectrum for K^ir^ at 
4010 MeV. The data is the points with error bars and the his- 
tograms are MC 1203 
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8.15 Sideband-subtracted momentum spectrum for K^ti^ti^ at 
4010 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.16 Sideband-subtracted momentum spectrum for —>■ K^n^ at 
4015 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.17 Sideband-subtracted momentum spectrum for K^-k^-k^ at 
4015 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.18 Sideband-subtracted momentum spectrum for K~7t^ at 
4030 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.19 Sideband-subtracted momentum spectrum for K~7t~^7[~^ at 
4030 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.20 Sideband-subtracted momentum spectrum for K^tc^ at 
4060 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.21 Sideband-subtracted momentum spectrum for K^n^ii^ at 
4060 MeV. The data is the points with error bars and the his- 
tograms are MC 
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8.22 Sideband-subtracted momentum spectrum for D'^ K^tt^ at 
4120 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.23 Sideband-subtracted momentum spectrum for K^ir^n^ at 
4120 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.24 Sideband-subtracted momentum spectrum for K^n^ at 
4140 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.25 Sideband-subtracted momentum spectrum for K^-k^-k^ at 
4140 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.26 Sideband-subtracted momentum spectrum for K'n^ at 
4160 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.27 Sideband-subtracted momentum spectrum for K^-k^-k^ at 
4160 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.28 Sideband-subtracted momentum spectrum for —>■ K^tt^ at 
4180 MeV. The data is the points with error bars and the his- 
tograms are MC 



8.29 Sideband-subtracted momentum spectrum for K^ti^ti^ at 
4180 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.30 Sideband-subtracted momentum spectrum for —>■ K^n^ at 
4200 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.31 Sideband-subtracted momentum spectrum for K^-k^-k^ at 
4200 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.32 Sideband-subtracted momentum spectrum for K~7t^ at 
4260 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.33 Sideband-subtracted momentum spectrum for K~7t~^7[~^ at 
4260 MeV. The data is the points with error bars and the his- 
tograms are MC 

8.34 Sideband-subtracted momentum spectrum for 07r+ at 3970 
MeV. The data is the points with error bars and the histograms 
are MC 

8.35 Sideband-subtracted momentum spectrum for </)7r+ at 3990 
MeV. The data is the points with error bars and the histograms 
are MC 



8.36 Sideband-subtracted momentum spectrum for Df 07r+ at 4010 
MeV. The data is the points with error bars and the histograms 
are MC 

8.37 Sideband-subtracted momentum spectrum for 07r+ at 4015 
MeV. The data is the points with error bars and the histograms 
are MC 

8.38 Sideband-subtracted momentum spectrum for — > 07r+ at 4030 
MeV. The data is the points with error bars and the histograms 
are MC 

8.39 Sideband-subtracted momentum spectrum for — > 07r"^ at 4060 
MeV. The data is the points with error bars and the histograms 
are MC 

8.40 Sideband-subtracted momentum spectrum for (f)7i~^ at 4120 
MeV. The data is the points with error bars and the histograms 
are MC 

8.41 Sideband-subtracted momentum spectrum for 07r+ at 4140 
MeV. The data is the points with error bars and the histograms 
are MC 

8.42 Sideband-subtracted momentum spectrum for </)7r+ at 4160 

MeV. The data is the points with error bars and the histograms 

are MC 
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8.43 Sideband-subtracted momentum spectrum for Df 07r+ at 4180 
MeV. The data is the points with error bars and the histograms 
are MC 

8.44 Sideband-subtracted momentum spectrum for 07r+ at 4200 
MeV. The data is the points with error bars and the histograms 
are MC 

8.45 Sideband-subtracted momentum spectrum for — > 07r+ at 4260 
MeV. The data is the points with error bars and the histograms 
are MC 

A.l Generator-level angular distribution for D*D events. 6 is defined 
as the angle between the D* and the beam axis 

A. 2 Generator-level angular distribution for D in D* for D*D 
events. 9' is defined as the angle between the D in the rest frame 
of the D* and the momentum vector of the D* in the lab 

A. 3 Generator-level angular distribution for D in D* Dn for D*D 
events. 6' is defined as the angle between the D in the rest frame 
of the D* to that of the momentum vector of the D* in the lab. . 

A. 4 Generator-level angular distribution for D*D* events. 6 is defined 
at the angle between the D* and the beam axis 
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A. 5 Generator-level angular distribution for D in D* ^ D'y for D*D* 
events. 9' is defined as the angle between the D in the rest frame 
of the D* to that of the momentum vector of the D* in the lab. . 

A. 6 Generator-level angular distribution for D in D* — > Dir^ for D*D* 
events. 9' is defined as the angle between the D in the rest frame 
of the D* to that of the momentum vector of the D* in the lab. . 

A. 7 Momentum distribution for in D*D* events. The different his- 
tograms correspond to different helicity amplitudes, and therefore 
different values for a in Eq. IA.21[ The solid histogram has a = 0.65 
and a = —0.24 for the pion and gamma decay respectively. The 
dashed histogram has a = 5.0 and a = —0.71 for the pion and 
gamma decay respectively. 

A. 8 This has the same parameters and values as Fig. IA.7I but has been 
convoluted with a gaussian of 10 MeV to account for detector res- 
olution. This distribution should be compared to Fig. I8.18[ . . . . 

A. 9 Generator level scatter plot showing the expected dependence 

for e+e" ^ 7* ^ D*D where D* Dn 

B. l a, defined in 1 -|- acos^^, as a function of center-of-mass energy 

from MC. The solid black lines correspond to the MC input value. 

Top plot presents the results of D in DD and the bottom plot of 

D* in D*D* 
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B.2 a, defined in 1+a cos^ 9, as a function of center-of-mass energy from 
Data. The red and blue lines correspond to the average value for 
—>■ K'TT^ and —>■ K^tt^tt^ respectively. Top plot presents 
the results of D in DD and the bottom plot of D* in D*D*. . . . 1260 

B.3 Plot of the production angle for D* in D*D* events for data at 
4170 MeV. The solid line shows the best fit to the data using the 
fit function 1 + acos^(6'), where 9 is the production angle. Based 
on the fit a = -0.62 ± 0.01 |263 

B.4 The momentum spectrum. The points with error bars is the 
data collected at 4170 MeV and the solid line is the best fit result. 
The fit function consisted of Eq. IA.21I convoluted with a gaussian 
resolution function. The result of the fit yields a = 0.80 ± 0.06. . |266 

B.5 Plot of the production angle for D* at 4170 MeV. Data is shown 
along with a generated MC sample using the helicity amplitudes 
described in the text 1268 

B.6 Plot of the angle between the D in the rest frame of the D* and 
the momentum of the D* in the lab for D*D* events at 4170 MeV. 
Data is shown along with a generated MC sample using the helicity 
amplitudes described in the text 1269 
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Chapter 1 
Introduction 



This dissertation is devoted to the study of charm-meson production in e'^e~ 
annihilations at thirteen center-of-mass energies between 3.97 GeV and 4.26 GeV. 
Specifically, we have used the CLEO-c detector to measure exclusive cross sections 
for several final states and the inclusive cross sections for the production of the 
charmed mesons D^, and Dg. To provide the background for this research I 
first review some of the foundations of elementary particle physics, including the 
Standard Model of quarks and leptons. 

During the past century, our understanding of nature and structure of mat- 
ter has been dramatically transformed, beginning with the discovery that atoms 
are divisible structures composed of subatomic particles. With this discovery 
and others the modern picture began to form of atoms consisting of a positively 

charged nucleus surrounded by negatively charged electrons. 
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Elementary particle physics is the study of the fundamental constituents of 
matter and their interactions. It began with the discoveries of the electron (Thom- 
son), atomic nucleus (Rutherford) and the neutron (Chadwick), developing slowly 
at the same time that the theoretical foundations of quantum mechanics were 
being laid. Meanwhile, investigation of cosmic rays and the development of 
accelerator-based experiments in the 1940's, 1950's, and the 1960's led to the 
discovery of many more particles, particles that do not exist in ordinary mat- 
ter. With the seemingly endless additions of new particles, the sentiment grew 
that there must be another level of substructure to explain this particle "zoo" . 
This sentiment ultimately led to the invention, by Gell-Mann and Zweig [1], of 
a model that describes most observed particles as being composed of the more 
fundamental particle known as the quark. 

1.1 The Standard Model 

The modern framework which describes the fundamental particles and interac- 
tions is called the Standard Model. The Standard Model incorporates the quarks 
and the leptons (the electron and its relatives) into a successful framework that 
has proven to be very rugged and reliable over the years. The quarks and leptons 
are said to be fundamental because they are structureless: they are point-like 
and indivisible. Even though they are structureless, however, they still possess 
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intrinsic properties such as spin, charge, color, etc. All particles but one can be 
arranged into three groups: quarks, leptons, and force carriers or mediators. The 
remaining particle, the Higgs boson, has a special job in the model. The Higgs is 
the last remaining particle of the Standard Model yet to be discovered. It plays 
the key role in explaining the mass of the other particles, specifically the large 
mass difference between the photon, the vector bosons and the quarks. 

The leptons come in six types which are listed in Table 11.11 along with their 
corresponding masses and charges. The lightest charged lepton is the familiar 
electron. The muon (/x) and tau (r) have the same general properties as the 
electron, but they have larger masses. These heavier versions are unstable and 
therefore not found in ordinary matter. Each charged lepton is accompanied by 
a weakly interacting electrically neutral neutrino: z/g, z/^, and Vr- Together each 
lepton and its accompanying neutrino form a family, sometimes referred to as a 
generation. All leptons have intrinsic angular momentum, or spin, of ^fi and are 
therefore fermions. 

Similarly, the quarks also come in six types, which can also be split into three 

generations. The first generation of quarks consists of the up (m) and down {d). 

Together with the electron the first generation of quarks make up all the ordinary 

matter around us. The second generation (strange (s) and charm (c) quarks) 

and the third generation (bottom (6) and top (t) quarks) make-up the rest of the 

known quarks. Since they have intrinsic angular momentum of all quarks 
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Table 1.1: The six types of leptons in the Standard ModeL The masses are taken 
from Ref. [2]. The electric charges are given in units of |e|, where e is the charge 
of an electron (|e| = 1.602 x 10~^^ Coulombs). The neutrino's masses are given 
at 95% confidence level. 



Family 


Name 


Electric 


Mass 






Change |^ 




I 


e 


-1 


511 keV 









<3 eV 


II 




-1 


106 MeV 









<0.19 MeV 


III 


r 


-1 


1.78 GeV 









<18.2 MeV 



are fermions. Some properties of the quarks are shown in Table II. 2[ Quarks 
are never individually observed but rather are seen only in combinations called 
hadrons. There are two types of hadrons: mesons, which are bound states of a 
quark and an anti-quarkj^ and baryons, which are bound states of three quarks 
or anti-quarksc 

The interactions among these fundamental particles are the four forces: strong, 
weak, electromagnetic, and Gravitational. Of these. Gravity is by far the weakest, 
plays no role in the interactions described in this dissertation, and is not part of 
the Standard Model. In the Standard Model the three remaining forces are rep- 
resented as an exchange of gauge bosons between interacting particles. All gauge 
bosons have integer spin, and some of their properties are shown in Table 11.31 



^For example, a tt"*" meson is a bound state of an up-quark and an anti-down-quark (ud) 
^For example, a proton (p) is a bound state of two up-quarks and a single down-quark (uud). 
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Table 1.2: The six types of quarks in the Standard Model. The masses are taken 
from Ref. [2]. The electric charges are given in units of |e|, where e is the charge 
of an electron (|e| = 1.602 x 10~^^ Coulombs). The quark masses are approximate 
since the complicated strong interactions that take place inside the hadrons make 
it difficult to define and measure the individual quark masses. However, because 
of its short lifetime, the top quark does not have time to bind into hadrons, and 
therefore it has the smallest fractional uncertainty on its mass. 



Family 


Name 


Electric 
Change |^ 


Mass 


I 


u 
d 


^3 
1 


1-4 MeV 
4-8 MeV 


II 


c 
s 


4 


1.15-1.35 GeV 
80-130 MeV 


III 


t 
b 


i 

3 


174 GeV 
4.1-4.4 GeV 



Of the fundamental fermions, only the quarks interact through the action of all 
three forces. The charged leptons experience only the weak and electromagnetic 
forces, and the neutrinos interact only via the weak force. 



Table 1.3: The gauge bosons in the Standard Model. The masses are taken from 
Ref. [2]. The electric charges are giving in units of |e|, where e is the charge of 
an electron (|e| = 1.602 x 10~^^ Coulombs). 



Force 


Name 


Electric 


Mass 






Change 




strong 


gluons g 








weak 


W-Bosons 


±1 


80.4 GeV 




Z-Boson 





91.2 GeV 


electromagnetic 


Photon 7 









The electromagnetic force is responsible for the binding of electrons to the 
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nucleus inside the atom. The mediator of this force, the photon, couples to 
electrically charged particles. The quantum theory describing the electromagnetic 
interaction is called Quantum Electrodynamics (QED) and has been extensively 
tested. QED is the most accurate physical theory constructed. 

The weak interaction was first discovered in nuclear decays. All fundamen- 
tal fermions, including the neutrinos, participate in the weak interaction. The 
mediators of the weak force are the vector bosons, the and the Z, which 
are massive, explaining its very short range, since the range is proportional to 
1/M. Interactions involving the charged VT-bosons involve changes the flavor of 
the particles involved. These changes can lead to spontaneous decays of the type 
fj,~ e~h'e^fi- Since the decay involves the charged VT-boson, these types of 
interactions are known as "charged current" interactions. 

The strong force involves the exchange of gluons between particles possessing 
"color charge" [3], which is analogous to the electric charge. The only fundamental 
particles possessing color charge are the quarks and gluons. Therefore, the leptons 
do not participate in strong interactions. Each quark carries one of the three 
colors, usually denoted red (r), green (^f), and blue (6), while the anti-quarks 
carry anti-color: f, g, and b. The quantum theory which describes the strong 
force is known as quantum chromodynamics (QCD). 

One interesting aspect of QCD is color confinement, the feature that no free 

quarks exist in nature. Rather, quarks must bind into the composite particles 
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known as hadrons, are color-neutral. Therefore the anti-quark of the meson will 



have the anti-co 
a different color 



or of its quark partner, while for a baryon each quark will possess 
3 The origin of the confinement is related to the fact that gluons 
themselves carry color and therefore interact with each other. As two quarks 
are pulled apart, the gluon-gluon interactions form a narrow tube resulting in a 
constant force. As the separation of the two quarks increases, so does the energy, 
until at some distance it becomes energetically favorable to "pop" out quarks from 
the vacuum to "dress" the quarks into hadrons. Confinement is not a feature of 
QED since the photon is electrically neutral. That is, as two charged particles 
get pulled apart, the force between them decreases allowing atoms to ionize. 

Another interesting aspect of QCD is asymptotic freedom [4J. As the momen- 
tum transferred in an interaction increases, the strength, or coupling, in the inter- 
action decreases. Therefore, in these types of interactions, perturbative methods 
can be applied to the calculations, as they are in QED. However, because the 
strength of the interaction increases at lower momenta, QCD calculations are 
very difficult in most cases. 

Since this thesis is devoted to the production of charm mesons, a dedicated 
review of charm follows in the next section. 



While the color of the quarks has nothing to do with color in visual perception, the similarity 
of the rules of combinations to those of color theory make the label an apt one. 
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1.2 Charm 

The quark model, when it was introduced by Gell-Mann and independently by 
Zweig in 1964, explained all hadrons matter as being combinations of three kinds 
of quarks: up (m), down (rf), and strange (s). Initially, there was no experimental 
reason for any additional quarks, but Bjorken and Glashow [5J, as a way of making 
nature more symmetrical at a time when there were three known quarks and 
four known leptons, predicted the existence of the yet-to-be discovered charm (c) 
quark. The Glashow-Iliopoulos-Maiani (GIM) mechanism [6l[7] made the case for 
charm more compelling, explaining the experimentally unobserved strangeness- 
changing neutral currents (SCNC) in the theory by adding charm in a particular 
way to the weak hadronic current which would otherwise be expected. 

In 1963, Cabibbo had noticed that the decay rates of strange particles differed 
from those of non-strange particles [8]. For example, the decay rate for K~ — >■ 
/x^z/^ was different from tt^ even after correcting for phase space. Since 

the quark content of the particles was different, ud for the vr" as compared to us 
for the K~ , Cabbibo hypothesized that, in the weak interaction, the u couples 
to a d quark, which is a superposition of the physical d and s quarks. To be 
explicit, d! = dcosOc + ssin^c, where 6c ~ 13° is known as the Cabibbo angle. 
The transition probability of a d quark changing to a m quark is proportional to 
cos^ 9c whereas it is sin^ 9c, for a s quark changing to a m quark. This leads to 
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Figure 1.1: The decay of f^^fJ'^ and the GIM mechanism. The only 

difference between a) and b) is the virtual quark, that is the u quark is replaced 
by the c when going from a) to b). The decay vanishes with the addition of the 
charm quark, thus eliminating strangeness-changing neutral currents as demanded 
by the data. 



the difference in the decay rates. 

At the time, and assuming no charm quark, the decay /^^A*" "was 

calculated to have a rate much larger (oc sin 6'c cos ^c) than the experimental 
limits would allow {Kg /i^A^^ < 3.2 x 10^'' However, with the addition of 
the charm quark, another diagram, or term, enters into the calculation, canceling 
the decay in the limit of SU(4) symmetryo Therefore, not only do the s and 
d quark couple to the u quark, but they now also couple to the c quark, where 
s' = s cos 6c — dsinOc- Now, the transition probability of a s quark changing to 
a c quark is proportional to sin^ 9^, whereas it is cos^ 9^ for a s quark changing 



'*The cancellation is not exact because the masses of the u and c quarks are not the same. 
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to a c quarklfl Notice the c quark diagram contribution (oc — sin Of. cos 9^) to the 
amphtude of Kg fJ'~^f^~ , meaning that SCNC are canceled with the addition of 
the c quark. Figure [LT] shows the two diagrams for the decay /^^/U^Jj In 

addition to the charm quark explaining the cancellation of SCNC, it also explained 
the smallness of the Ks-Kl mass difference, which arises because of mixing. This 
difference, gives one an idea, in broken SU(4), of the mass of the charm quark, 
about 1.5 GeV [lOj. Now, assuming the binding energy is small relative to the 
mass of the charm quark, a meson composed of a charm quark should also have 
a mass close to the charm quark mass. 

The so-called "November Revolution" followed the 1974 discovery of the J/ip 
resonance, and the subsequent exploration of the charm sector resoundingly vali- 
dated the quark model. The computational power of the quark model was further 
verified in that the closed-charm (cc) states J/ ip and ip' had widths in agreement 
with the Okubo-Zweig-Iizuka (OZI) suppression hypothesis [TT\ . 

1.2.1 D-mesons 

The theory of weak decays of quarks, as formulated by Cabibbo and extended 
by GIM, predicted that most of the charmed particles should have a K meson 
in the final state. This can be viewed as the charm quark transforming into a 



^This idea has been expanded to incorporate all three generation of quarks by Kobayashi 
and Maskawa [9]. 

nice review can be found in D. Grifhths "Introduction to Elementary Particles" [7]- 
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strange quark followed by the dressing of the quarks into mesons. Gaillard, Lee, 
and Rosner [12] proposed that charmed mesons, either D^{cu) or D'^{cd), could 
be observed by looking for peaks in the invariant-mass spectrum of K'tt^ or 
K^n^n^. Diagrams of these decays are shown in Figs. 11.21 and [Ol 



Figure 1.2: Quark diagram for D K . The light u quark does not partici- 
pate in the decay and is referred to clS cl "spectator" . 

The first detection of a D-meson decay was made by the Mark I collaboration 
[13] in the modes K-n+, D° ^ K-n+n-n+, and D+ K-n+n+. In 

addition to these charm mesons, the Mark I collaboration also observed the lowest 
charm-meson excited states of D*° and D*+ [H]. 

In addition to a charm quark combining with an up or a down quark to form 
a meson, it can also combine with a strange quark to form what is referred to 
as a Dg meson. The decay of this type of meson can be characterized by the 




K 
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Figure 1.3: Quark diagram for K 7r+7r+. The light d quark does not 

participate in the decay and is referred to as a "spectator" . 

conversion of the charm quark into a strange quark with an emission of a virtual 
W^-boson. The simplest possible final products result when the newly-formed 
strange quark pairs with a strange quark to form a or an meson and the 
virtual VT-boson forms a charged pion. Therefore, a search for an excess of rjTi 
or </)7r production would be evidence of Dg mesons. The discovery of Dg (initially 
called F) was presented by the CLEO collaboration in 1983 by observing a peak 
in the invariant mass spectrum for the decay mode (pir [15] , following a period of 
considerable experimental uncertainty |ltij . 
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1.3 Motivation for this Measurement 

The years following the charm discovery provided few opportunities for its detailed 
study, especially the Ds- The CLEO-c experiment was proposed in 2001 as to 
provide the first high-statistics study of charmed particles with a state-of-the-art 
detector. Initial CLEO-c running beginning in late 2003, focused on studies of 

and with the expectation that Dg studies would follow. 

From August until October 2005, the CLEO-c collaboration carried out a scan 
of the center-of-mass energy range from 3970 to 4260 MeV. The main purpose 
of this scan was to determine the optimal running point for CLEO-c studies 
of Dcj-meson decays. Secondary objectives included detailed measurements of 
the properties of charm production in this region, tests of previous theoretical 
predictions, and confirmation and additional studies of the y (4260) state reported 
by the BaBar collaboration [17]. In this dissertation, measurements with the 
scan data of the cross sections for inclusive hadron production, charmed-hadron 
production, and the production of specific final states including charmed mesons 
are presented. We present results both without and with radiative corrections. 
We also report initial studies with detailed features of these events and present 
comparisons with past measurements and theoretical expectations. All of which 
will add to our present understanding for systems containing heavy and light 
quarks. 
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1.3.1 Previous Experimental Results 

The total hadronic cross section in e~^e~ annihilations is generally presented in 
terms of R, which is defined as follows: 

a(e~^e~ —>■ hadrons) 

R= \ + _ 

Incorporating the well-known /x-pair cross section with its 1/s energy dependence 
gives 

/ + - N + - +-N ..86.8 nb GeV^ 
0"(e e hadrons) = Ra[e e /i /i ) = R — . (l-^j 

cm 

If strong interactions are ignored: 

R = 2,Y,i (1-3) 

i 

where the summation is over all kinematically allowed quark fiavors at the energy 
of interest, qi is the quark's corresponding charge (either —1/3 or 2/3), and the 
factor are 3 is due to the fact that there is 3 different colors available. This 
approximation only holds in energy regions well above the qq bound states. 

Equation 11.11 counts the number of possible quarks that are available at a 
particular energy. In the energy region below, the J/ifj resonance, Mj/^ = 3.097 
[2], only the m, and s quarks contribute to R\ 

R = ^[{\f + {\f + {\Y] = '^ (1-4) 

Above this energy the c quark can also contribute to the ratio. 
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R has been measured over a very wide energy range by many experiments [2], 
including recent measurements with the Beijing Spectrometer (BES) [IB] in the 
energy range of interest for CLEO-c. Fig. II .41 shows the current state of knowledge 
of R at center-of-mass energies between 3.2 and 4.4 GeV. The data in Fig. II. 41 has 




Figure 1.4: (T(e^e~ — > hadrons) from the 2005 PDG [2j. The solid line is a - fit 
to the data and the dotted line is the extension of the fit to higher energies. The 
vertical lines indicate the various thresholds and resonances present in the region. 



been radiatively corrected. The correction depends on the center-of-mass energy 

and the behavior of the cross section at lower energies. There is a rich structure 

in this energy region, reflecting the production of cc resonances and the crossing 
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of thresholds for specific charmed-meson final states. Some of these "landmarks" 
are highlighted with vertical lines in Fig. II. 4[ The specific energies corresponding 
to these thresholds are tabulated in Table II. 4[ 



Table 1.4: Thresholds of interest in the e+e~ center-of-mass energy range of the 
CLEO-c scan. Threshold values are computed from the PDG |2|, and where 
needed charged and neutral states of the D were averaged. 



Center-of-Mass 
Energy 


3875 MeV 


3937 MeV 


4017 MeV 


4081 MeV 


4224 MeV 


State 


D*D 




D*D* 




d:d: 



There are two interesting features in the hadronic cross section between 3.9 
and 4.2 GeV. There is a large enhancement at ~ 4 GeV corresponding to the 
D*D* thresholdjj Next, there is a fairly large plateau that begins at the D*^Dj 
threshold. There is considerable theoretical interest and little experimental infor- 
mation about the specific composition of these enhancements. 

Prior to the CLEO-c scan run there were insufficient data on Ds production 
for an informed decision about the best energy at which to undertake CLEO-c 
studies of Ds decays. BES measured the inclusive Ds production cross section 
times the Df — > 07r+ branching ratio at the center-of-mass energy 4030 MeV to 
be (11.2 ± 2.0 ± 2.5) pb [19j. Since there is only one accessible final state with 
Ds at this energy, this measurement suggested a cross section of about 0.3 nb for 
D^D~ . The Mark III collaboration previously measured the same quantity at 



^Throughout this paper charge-conjugate modes are imphed 
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Figure 1.5: The production cross section times branching ratio for —>■ K n 
and K^n^n^ as a function of energy. 
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a center-of-mass energy of 4140 MeV to be (26 ± 6 ± 5) pb [20]. Dg production 
had previously been demonstrated to be dominated at this energy by Dl^D~ 
[2T] . The abihty to do Ds physics with CLEO-c depends both on the quantity of 
Dg production and on the complexity of the events. It was therefore essential to 
measure all accessible final states and carefully assess the physics reach for future 
Dg studies under the conditions prevailing at each energy. 

Studies of Dg production, when combined with measurements of and 
production, would constitute a comprehensive analysis of all charm production 
in the region just above threshold. There are more previous measurements of 
and production than of Dg, but here too the information is limited. 
BES and MARK II made measurements of cross section time branching ratio 
for D° K-n+ and D+ K-n+n+ [221 [221 [21 ES] , which are shown in Fig. 
11.51 Interpretation of these data points is complicated by the presence of several 
possible final states: DD, D*D, and D*D* in both charged states. Strong inter- 
action theory provides predictions of the overall cross sections and proportions of 
the various final states, some of which are described in the next section. Detailed 
measurements at several points would allow more rigorous testing of models of 
charm production in the region above cc threshold. 
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1.3.2 Theoretical Predictions 

Soon after the discovery of the D meson relative production ratios were predicted 
by counting available spin states [26] for the possible reactions: 

• 1" 0"0~, as in e+e" ^ Y ^ DD 

• I- ^ as in e+e" ^ 7* ^ DD* 

• I- ^ 1-1", as in e+e- 7* ^ D*D* 
This argument gives thefollowing ratios: 

(^DD '■ (^D*D '■ (^D*D* = 1:4:7 (1-5) 

This naive expectation disagrees with the experimentally observed ratios obtained 
by the Mark I experiment [H] at 4028 MeV: 

(^DD ■ (^D-D ■ (^D*D* = 0.2 ± 0.1 : 4.0 ± 0.8 : 128 ± 40, (1.6) 

where the phase space factor has been removed. The severe disagreement 
between what is expected and what is experimentally observed suggests that 
major additional effects are present. 

A serious theoretical calculation of the charm cross sections was first at- 
tempted by Eichten et al. in 1980 [27j with a coupled-channel potential model. 
Their predictions for the excess AR above uds production are shown in Figs. 11.61 

and ll.7[ According to these predictions the large enhancement in the cross section 
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at 4030 MeV is dominated by D*D and D*D* production. A detailed investi- 
gation of this energy region could definitively confirm or refute this prediction. 




*i7 fcO IjO 1-1 t-C Kl^a M 4.S 

Figure 1.6: The charm contribution to R as calculated in the coupled channel 
model of Eichten et al. [27] for events. 

The mass of the Dg meson used in Eichten's 1980 prediction is incorrect, as 
indicated by use of the older notation F for that particle. To update Eichten's 
prediction, it is therefore necessary to shift the cross sections downward in en- 
ergy by 144 MeV and 100 MeV for DfD~ and D*~^Dj, respectively. After this 
correction, Eichten predicts that the largest Dg yield should be at ~ 4160 MeV, 
where the dominant source is D^D*~ events. 
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Figure 1.7: The Dg contribution to R as calculated in the coupled channel model 
Eichten et al. [27| for of events. In the fi gure Ds is called F. 
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More recently, T. Barnes [28] has presented calculations, using the phenomeno- 
logical model [22], at 4040 and 4159 MeV, which are summarized in Table 11.51 
Of these two energy points, Barnes predicts that 4159 MeV is the better place for 



Table 1.5: Partial widths, in units of MeV, of charm production as predicted by 
Barnes for two center of mass energies. 



Center-of-Mass 
Energy 


DD 


D*D 


D*D* 






SUM 


Exp. 


4040 MeV 


0.1 


33 


33 


7.8 




74 


52±10 


4159 MeV 


16 


0.4 


35 


8.0 


14 


74 


78±20 



Ds physics, with a total cross section for Dg production that is three times larger 
than that at 4040 MeV. In addition, he finds the enhancement in R at ~4 GeV to 
be due to an equal mixture of D*D and D*D* events. While the experimentally 
measured summed rate differs by ~ 2a from Barnes's prediction at 4040 MeV, 
the precision is insufficient for a definitive conclusion. Precise measurements of 
the partial widths with CLEO-c would be decisive in testing this model. 
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Chapter 2 



Experimental Apparatus 

2.1 CESR - The Cornell Electron Storage Ring 

The Cornell Electron Storage Ring, or CESR, is located in central New York on 
the Cornell University campus. As shown in Fig. 12. ![ CESR consists of three 
basic parts: a linear accelerator (linac), a synchrotron, and the storage ring. The 
storage ring and synchrotron are housed in a circular tunnel which has a diameter 
of 244 meters. The linac is located in the inner part of the ring. The CLEO-c 
detector is located at the south end of the tunnel. 

Electrons for acceleration are created by boiling them off a heating filament. 
Once the energy is sufficient for them to escape from the surface, they are collected 
and bunched together for acceleration. The linac consists of oscillating electric 
fields which accelerate the electrons down the length of the linac. By the end, 
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1600799-005 

CESR 




CHESS CLEO CHESS 



Figure 2.1: The CESR accelerator facility. The main components of the facility 
are the linac, the synchrotron, and the storage ring. The linac converter is a piece 
of tungsten used to induce an electromagnetic shower so as to produce positrons. 
CHESS refers to the Cornell High Energy Synchrotron Source. 
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the electrons have an energy of approximately 200 MeV. 

Positrons are created by inserting a tungsten target, a converter, halfway down 
the linac. The target is used to generate electromagnetic showers consisting of 
electrons, positrons, and photons. The positrons in the shower are selected out 
and accelerated the rest of the way down the length of the hnac. By the end of 
the acceleration, the positrons reach an energy of 200 MeV. 

The bunches of the electrons and positrons from the linac are injected sepa- 
rately and in opposite directions into the synchrotron. The synchrotron consists 
of a series of dipolc magnets and four three-meter-long linear accelerating cavities. 
The dipole magnets steer the electrons and positrons around the ring while ac- 
celerating cavities increase the particles' energies to about 2 GeV. As the energy 
of the particles increases the magnetic fields of the dipoles are increased to keep 
the particles in the ring. Once accelerated to 2 GeV, the electrons and positrons 
are transferred to the storage ring. 

The storage ring consists of a series of a dipole magnets which steer the elec- 
trons and positrons around the ring, in addition to quadrupole and sextupole 
magnets which focus the beams. As the particles traverse the ring they lose en- 
ergy due to synchrotron radiation. The energy is replaced by superconducting 
radio-frequency (RF) cavities which operate at a frequency of 500 MHz. These 
RF cavities are similar to those used in the linac and synchrotron, except that 

they do not significantly accelerate the beams, but primarily replace the radiated 
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energy. 

In current running conditions, CESR operates with nine bunch trains each for 
the electrons and positrons. Each train consists of as many as five bunches. To 
avoid unwanted interactions with the two counter-rotating beams four electro- 
static horizontal separators are used. These separators set up what are known as 
pretzel orbits and ensure that the electrons and positrons miss each other when 
they pass through the unwanted intersecting locations, sometimes referred to a 
parasitic crossing. A picture showing an exaggerated view of the pretzel orbits is 
shown in Fig. 12.21 At the interaction point, which is surrounded by the CLEO-c 
detector, the two beams are steered into each other. However, the two beams do 
not collide head-on, but rather at a small crossing angle of 2.5 mrad (~ 1.7°). 
This allows for bunch-by-bunch interactions of the electron and positron trains. 

Between 1979 and 2003, CESR operated at the T(4S') resonance, which cor- 
responds to a center-of-mass energy of 10.6 GeV or beam energy of 5.3 GeV. The 
CLEO-c charm program is carried out at center-of-mass energies between 3 and 5 
GeV, which required major changes to CESR. The rate of synchrotron radiation, 
energy radiated by the beams due to acceleration by the bending magnets, is pro- 
portional to [30]. The decrease in the amount of synchrotron radiation affects 
storage ring performance through two important beam parameters. One is the 
damping time with which perturbations in beam orbits caused by injection and 

other transitions decay away. In particular, a particle that is off the ideal orbit 
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0140895-023 




Figure 2.2: Schematic of CESR showning the "pretzel" orbits. The pretzel or- 
bits are used to separate the electron and positron beams at parasitic crossing 
locations. 
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because of a larger energy radiates slightly more energy, whereas a particle with 
a lower energy radiates slightly less. As a whole, the energy spread among the 
particles becomes reduced which shows up as a damping of the oscillations |31j . 
In the original CESR design the radius was chosen to ensure adequate damping 
through this mechanism at a 5-GeV beam energy. Since the amount of radia- 
tion is dependent upon the energy of the corresponding particles in the beam, 
with a reduction in beam energy to ~ 2 GeV, the damping time becomes too 
long for effective operation. The other parameter is the horizontal beam size, or 
horizontal emittance, which measures the spread of particles in the bend plane. 
The horizontal beam size, which results from the betatron oscillation in addition 
to the rate of quantum fluctuations in the synchrotron radiation, decreases with 
decreasing beam energy thereby limiting the particle density per bunch |31j . 

These effects limited the achievable collision rate (luminosity) of CESR at 
the lower center-of-mass energies of CLEO-c to unacceptable levels. To increase 
the luminosity CESR accelerator physicists proposed to increase the amount of 
radiation through the use of wiggler magnets [30] . 

A wiggler magnet is a series of dipole magnets with high magnetic fields. Each 

successive dipole has its direction of magnetic field flipped. Therefore, when a 

particle passes, it will oscillate, which results in emission of additional synchrotron 

radiation without changing the overall path of the particle around the ring. As a 

result, the damping time is decreased while the horizontal beam size is increased, 
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thereby increasing the luminosity. 

CESR has installed twelve superconducting wiggler magnets for low-energy 
running. Each wiggler consists of eight dipole magnets with a maximum field 
strength of 2.1 Tesla ^30j. These wigglers decrease the damping time by a factor 
of 10, while increasing the beam emittance by a factor of 4 to 8 [32]. In addition, 
the energy resolution has increased, as compared to no wigglers, by a factor of 4 
to as/E = 8.6 X 10"^ 

The center-of-mass (CM) energy of the collision is an important quantity 
needed to put the cross section measurements into context. In order to determine 
the CM energy of the colliding electrons and positrons, one needs to know the 
energy of the beams. The energy of the beams, to first order, can be determined 
by the following [55] : 



where c is the speed of light in a vacuum and e is the charge of an electron. 
The summation in Eq. 12.11 is over all dipole magnets in the ring, where Bi, 6i, 
and pi are the magnitude of the magnetic field, the bending angle, and bending 
radius of curvature of the i-th dipole magnet, respectively. The result of Eq. 12.11 
needs to be corrected for shifts in RF accelerating cavities, for the currents of the 
the focusing and steering magnets, and for the horizontal separators. The total 
uncertainty in the CM energy is of order 1 MeV. 




(2.1) 
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Another quantity needed in determining the cross sections is the luminosity, 
which quantifies the rate of e'^e~ coUisions. The number of events expected for a 
particular process is given by 

Ni = ai J Cdt, (2.2) 

where C is the instantaneous luminosity and (Xj is the cross section for that process. 
The integral of the instantaneous luminosity, / Cdt = L, is the quantity needed 
for this analysis and is referred to as the integrated luminosity or just luminosity. 
In CLEO-c, three final states are used to obtain the luminosity. The processes 
e'^e~ — > e+e~, fJ^'^fi', and 77 are used since their cross sections are precisely 
determined by QED. Each of the three final states relies on different components 
of the detector, with different systematic effects [31]. The three individual results 
are combined using a weighted average to obtain the luminosity needed for this 
analysis. 

2.2 The CLEO-c Detector 

The electron and positron beams are steered into each other at the central location 
known as the interaction point (IP) of the CLEO-c detector. CLEO-c is a general- 
purpose detector designed to identify and measure relatively long-lived charged 
and neutral particles. The charged particles detected are the electron, muon, 

pion, kaon, and proton. While some neutral particles are easy to identify, like the 
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photon, others hke the neutrinos are impossible. 

The remainder of this chapter is dedicated to describing the detection systems 
that make up the CLEO-c detector. As shown schematically in Figs. 12.31 and \2A\ 
CLEO-c is a cylindrically symmetric detector aligned along the beam line, the 
z-axis, that covers 93% of the solid angle. Starting from the IP and moving out, 
the detector is composed of an inner drift chamber [30], an outer or main drift 
chamber [301 ES], a Ring- Imaging Cherenkov detector [201 ESI EZ], denoted by 
the acronym RICH, a crystal calorimeter [301 138] , and finally a muon detector 
[5U1 EHl EH]- All systems, except for the muon detector, are within a uniform 
1-Tesla magnetic field, produced by a superconducting solenoid, which is aligned 
with the beam axis. 
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Figure 2.3: The CLEO-c detector. 
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Figure 2.4: A slice in the r-z plane of the CLEO-c detector. 
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2.2.1 Inner Drift Chamber 

The inner drift chamber (ZD) is the innermost detector of CLEO-c, located right 
outside the beam pipe. The inner drift chamber replaced the silicon vertex detec- 
tor of CLEO III because multiple scattering degrades the momentum resolution 
for soft tracks, which are more prevalent at low center-of-mass energies. Therefore, 
minimizing the material is crucial. The goal of the detector is to detect charged 
particles with | cos 6*1 < 0.93, where 9 is defined as the angle of the particle with 
respect to the beam pipe, the z-axis. It consists of helium-propane-gas-filled vol- 
ume segmented into 300 cells (half-cell size of 5 mm), where each cell consists of 
a sense wire held at +1900 V. These cells are surrounded by field wires, held at 
ground, which neighboring cells share. When a charged particle travels through 
a cell the gas is ionized. The resulting ionized electrons travel away from the field 
wire and toward the sense wire. Since the electric field increases close to the sense 
wire, the primary electrons will ionize other atoms in the gas, thereby creating 
an avalanche of electrons at the sense wire. The time of the resulting electric 
pulse seen on the sense wire, which is synchronized with the bunch crossing, is 
converted using the drift velocity of the gas to a distance of closest approach to 
the sense wire. This information can then be used to map out the trajectories of 
the charged particles through the drift chamber. 
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2.2.2 Outer Drift Chamber 

Located directly outside the inner drift chamber is the main drift chamber, some- 
times referred to as the outer drift chamber (DR). The main drift chamber has 
a similar design to the ZD, except that it is larger in size both overall, with 47 
layers of field and sense sires as compared to 6 for the ZD, and in the size of the 
cell, with a half-cell size of 7 mm as compared to 5 mm. In addition, the field 
wires are held at a higher potential, +2100 V rather than +1900 V. 

The energy lost by a charged particle in the drift chamber is used in identify- 
ing what particle traversed the volume. The energy lost per unit length, dE/dx, 
is related to the particle's velocity. By constructing a x^-like variable, the consis- 
tency of the actual energy lost per unit length with what is expected for different 
particle hypotheses can be assessed: 

^ _ dE/dx~''''''^ - dE/dxT"""'"^ ^2 3) 

where i = e, /x, vr, K, or p. The quantity a is the uncertainty in the measurement, 

usually approximately 6%. Fig. 12.51 shows the measured dE/dx as a function of 

particle momentum for a large number of charged particles detected in CLEO-c. 

The figure clearly shows good n — K separation for momenta below 500 MeV. At 

momenta above 500 MeV, dE/dx is quite limited and the additional separation 

power of the RICH is needed. 

As mentioned earlier, the ZD and DR are contained inside a 1-Telsa magnetic 
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Figure 2.5: dE/dx as a function of momentum for a large population of charged 
particle measured in CLEO-c. The tt, and p bands, from left to right, are 
clearly seen. 
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field oriented along the beam. This field causes charged particles to travel in 
helical paths as they travel through the detector volume. In other words, the 
particle's path is circular in the x-y plane and moves with constant velocity along 
the z-axis. Pattern recognition computer programs group the hits into tracks 
and fitting programs are used to obtain the parameters [IQ]. The CLEO-c fitter 
is an implementation of the Billoir or Kalman algorithm, which incorporates 
the expected energy loss of a particular particle to optimize the determination 
of its momentum and trajectory. At 1 GeV/c the charged-particle momentum 
resolution is approximately 0.6%. 

2.2.3 RICH - Ring Imaging Cherenkov Detector 

Directly outside the main drift chamber is the Ring Imaging Cherenkov Detector, 
commonly referred to as the RICH. Radiation is emitted when a charged particle's 
velocity is greater than the speed of light in the medium through which the particle 
is traveling; this radiation is known as Cherenkov radiation. The radiation is 
emitted in a cone with a characteristic opening angle known as the Cherenkov 
angle. The Cherenkov angle is related to the velocity of the particle by 



where (3 is the velocity in units of c and n is the index of refraction of the medium. 



cos6' 



1 



(2.4) 



Using (3 



^, in addition to E"^ 



m? + p"^, Eq. 12.41 can be rewritten in terms of 
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the particle's mass and momentum as follows: 

cose = -Jl + —. (2.5) 
nV 

This shows that one can identify the type of particle by using its momentum from 
the fitted track and the observed Cherenkov angle. 

The RICH detector is shown schematically in Fig. 12.61 It covers approximately 
83% of the full 4:71 solid angle. Cherenkov photons are produced when a charged 
particle passes through entrance windows fabricated from LiF crystals, as shown 
in Fig. l2.6[ There are a total of 14 rows, or rings, of these radiator crystals. All but 
the central four have flat surfaces, while the remainder have a "sawtooth" surface 
to reduce loss of photons by total internal reflection. The photons, which have a 
typical wavelength A = 150 nm, travel through an expansion volume filled with 
nitrogen gas which is effectively transparent. After traveling through the expan- 
sion volume they pass through a CaF2 window and into a multi-wire proportional 
chamber (MWPC) filled with a methane-TEA (tri-ethyl amine) mixture. Here 
photoelectrons are created which are collected in the same manner as described 
above for the drift chambers. Examples of the Cherenkov rings are shown in Fig. 

Using Cherenkov photon images, like those in Fig. \2.7\ one can construct the 
likelihood for a particular particle hypothesis. A x^-like variable for identification 
can be constructed to discriminate between two different particle hypotheses: 
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gure 2.6: A section of the CLEO-c RICH detector shown in r — cross section. 
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Figure 2.7: Cherenkov rings produced by a charged track for the flat radiator 
(right) and sawtooth (left). The central hits of each image are due to the charge 
track passing through the wire chamber. The other hits arc due to the Chernkov 
photons that are produced. Half of the image of the fiat radiator is missing 
because it was trapped in the radiator by total internal reflection. The sawtooth 
imagine is distorted by refraction. 
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X,'-X- = -21nL, + 21nL,-, (2.6) 

where Lj and Lj are the hkehhoods for particle hypotheses i and j, respectively. 
An illustration of the power of the RICH detector is shown in Fig. 12.81 Requiring 
x|- ~ < 0, that is that the particle is more kaon-like than pion-like, and that 
momentum is greater than 700 MeV/c, leads to a kaon identification efficiency of 
92% with a pion-fake rate of 8%. Fig. 12.91 shows particle separation as a function 
of momentum for different particle hypotheses above their respective thresholds, 
where the threshold is determined by the index of refraction of the LiF radiator, 
n = 1.4. 



41 




Figure 2.8: Kaon identification efficiency (solid points) and pion fake rate (open 
circles) as a function of various cuts on Xk — xi- Momentum is restricted to be 
above 700 MeV/c. 



42 



0140501-006 




0.5 1.0 1.5 2.0 

P (GeV/c) 

Figure 2.9: Particle separation as a function of momentum for different particle 
hypotheses in the RICH detector. Both particles are restricted to be above their 
respective threshold, where the threshold is determined by the index of refraction 
of the LiF radiator, n — 1.4. 
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2.2.4 Calorimetry 

Located just outside the RICH and just inside CLEO-c's superconducting magnet 
is the electromagnetic crystal calorimeter (CC). The CC consists of about 7,800 
thallium-doped cesium iodide crystals. About 80% of the crystals are arranged 
in a projective geometry in the barrel region, defined by |cos^| < 0.8. The 
remainder are in two end-caps, covering 0.85 < |cos^| < 0.93. The transition 
regions between the barrel and end-caps provide substandard performance and 
are generally not used. 

When charged particles or photons enter these highly dense crystals they 
interact and lose energy through various mechanisms: ionization, bremsstrahlung, 
pair conversion, and nuclear interactions. Electromagnetic interactions with these 
high-Z nuclei are very effective in stopping electrons and photons. Hadrons, on 
the other hand, lose energy less quickly in the CC electromagnetically, and their 
hadronic showers extend into the steel of the flux return of the superconducting 
solenoid. Muons and noninteracting hadrons, which deposit only a small fraction 
of their energy inside the calorimeter, and are referred to as minimum ionizing 
particles (MIPS). While hadrons generally are absorbed in the steel, the muons 
continue to lose energy only by ionization and travel through the magnet and into 
the proportional chambers that comprise the muon detector. 
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Electrons and photons lose energy through the successive generation of brems- 
strahlung photons and e^e~ production, together known as an electromagnetic 
shower. These showers produce numerous low-energy electrons which are then 
captured by the thallium atoms. The photons emitted by the de-excitation of 
thallium, A = 560 nm, are invisible to the rest of the crystal. This means, they 
can propagate through the rest of the 30 cm-long crystal and be collected by the 
photo-diodes mounted on the back of the crystal. The energy resolution is about 
4.0% at 100 MeV and 2.2% at 1 GeV. The resolution of the total shower energy 
is increased when more than one crystal is used in the reconstruction. This can 
be seen in Fig. 12.101 The center of the shower is then determined by an energy- 
weighted average of the blocks used in the sum. The number of crystals used is 
logarithmic, and ranges from 4 at 25 MeV to 13 at 1 GeV. 
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Figure 2.10: Shower energy resolution as a function of the blocks used in the 
shower reconstruction for the CLEO II detector. The barrel calorimeter has not 
changed since the installation of CLEO II. The smooth lines are from Monte-Carlo 
(MC) simulation of the 100 McV photons. The points are experimental data from 
the 100 MeV transition photon in T(3S') 7X;,j(2P). The arrow indicates the 
number of crystals used in the shower reconstruction of 100 MeV photons. 
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Figure 2.11: M(77) resolution in CLEO III. 
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2.2.5 Magnetic Field 

CLEO-c's 1-T magnetic field is provided by a large liquid-helium-cooled super- 
conducting solenoid with a diameter of 3 m and length of 3.5 m. The resulting 
field produced is uniform to ±0.02% over the entire tracking volume. The iron 
flux return of the magnet is used in muon detection as an absorber, which is 
described next. 

2.2.6 Muon Detection 

The muon detector exploits the fact that muons do not participate in the strong 
interaction. Since muons are much heavier than electrons, they lose energy much 
more slowly as they travel through material. Therefore, muons can penetrate 
much greater depths of material. The detector consists of three layers of gas-filled, 
wire-proportional tracking chambers in between 36 cm iron absorbers surrounding 
the detector (see Fig. 12.121) . 

The muon detector's layers provides information on a particle traveling differ- 
ent interaction lengths. The interaction length is the average distance a charged 
hadron has to travel before having an interaction. The three layers of the detector 
are located at approximately 3, 5, and 7 interaction lengths. Information from 
this detector is useful for identifying muons above about 1 GeV. Its applicability 
to the studies with CLEO-c reported in this thesis is limited. 
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Figure 2.12: Cross section of the muon chamber. It consists of three layers of 
8-cell proportional counters. The middle counter layer is offset by one half a cell 
width to improve acceptance. 
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2.2.7 Data Acquisition 

During the CLEO-c scan, bunch crossings happened at a rate on the order of 1 
MHz. However, the actual rate of interesting physics events was much smaller, 
on the order of 1 Hz. 

The CLEO-c detector includes hundreds of thousands of sensitive components 
and associated electronic channels. At every crossing each one of these compo- 
nents can deliver a signal representing the traversal of a particle produced in the 
e^e~ annihilation, although in typical events only a small fraction of these events 
have valid data. In general data are read out locally and "sparsified" to eliminate 
uninteresting channels. The nontrivial information is read out to on-line comput- 
ers in tens of microseconds. During the data-gathering process, the detector can 
not acquire any new events. Because of this, it is essential to record only those 
events that contain interesting physics. This amount of "dead-time" is defined as 
the time between the trigger signal and the end of the digitization process. The 
maximum readout-induced dead-time is < 3% [30] . 

The selection of these interesting physics events is achieved with a multi- 
layered trigger system [501 ISl 1121 113] • A schematic view of the CLEO-c trigger 
system is shown in Fig. I2.14[ Currently, there are eight trigger lines used, listed in 
Table 12. 1[ Data from the DR and the CC are received and processed in separate 
VME crates and yield basic trigger parameters. These parameters are tracks 
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Table 2.1: Definitions of the trigger lines. 



Name 


Definition 


Hadronic 


(iVaxW > l)Sz{NcBlov, > 0) 


Muon Pair 


two back-to-back stereo tracks 


Barrel Bhabha 


back-to-back high showers in CB 


End-cap Bhabha 


back-to-back high showers in CE 


Electron track 


{Na^ial > 0)L{NcBmed > 0) 


Tau 


(Nstereo > l)HNcBlow > 0) 


Two Track 


(Naa^ial > 1) 


Random 


random 1 kHz source 



counts, or multiplicity, topology in the main drift chamber, and the number of 
showers and topology in the calorimeter. The information from both systems is 
correlated by a global trigger which generates a "pass" signal every time a valid 
trigger condition is satisfied. The trigger system consists of two tracking triggers, 
one using information provided by the axial wires and the other using the stereo 
wires of the DR, a CC trigger, and a decision and gating global trigger system. 

The data acquisition system (DAQ) consists of two equally important partsll 
The DAQ system is responsible for the transfer of the data from the front end elec- 
tronics to the mass storage device, while the slow control monitors environmental 
conditions, the quality of the data, and the status of detector components. 

For each acceptable trigger from the CLEO-c detector, about 400,000 chan- 
nels are digitized. The front-end electronics provide the data conversion in par- 
allel with a local buffer and waits for an asynchronous readout by the DAQ. 
^The review of the data acquisition system is based upon Reference |30j. 
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Figure 2.13: Block diagram of the CLEO III data acquisition system (DAQ). The 
only difference in the between CLEO III and CLEO-c detectors, and in turn the 
DAQ system is the replacement of the silicon vertex detector with the inner drift 
chamber. 
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A dedicated module, the Data-Mover, in each front-end crate, assures transfer 
times of the data are below 500 fis, in addition to providing another buffer. The 
Data- Mover moves the data to the Event-Builder which reconstructs the accepted 
events for transfer to mass storage. Also, a fraction of reconstructed events are 
analyzed on-the-fly by the CLEO monitoring system, commonly referred to as 
Online-Passl, to quickly discover problems and check the quality. 

In addition to Online-Passl, there is an Of f line-Passl, which is commonly 
referred to as Caliper (which stands for CALIbration and PERformance monitor- 
ing). Caliper, unlike Online-Passl, can be run over all the data in an efficient 
manner by applying harsh cuts to get at the interesting physics events for data 
quality and monitoring. As a result. Caliper was used extensively during the 
scan running used for the analysis in this thesis. 

A block diagram of the DAQ for the CLEO III detector is shown in Fig. 12.131 
The only difference between the CLEO III and the CLEO-c detector, and that of 
the DAQ, is the replacement of the silicon vertex detector within the inner drift 
chamber. 
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Figure 2.14: Overview of the trigger system. 
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Chapter 3 



CLEO-c Ds Scan 



The CLEO-c project description (Yellow Book) [30] includes measurements of 
Dg branching fractions and other properties among the principal goals of the 
program. It makes a specific proposal that a scan run be undertaken to determine 
the running energy above D^D^ threshold that will maximize the sensitivity for 
Dg physics. It also acknowledges that such a data sample would include non- 
strange charmed mesons in the form of DD, D*D and D*D* events produced in 
different quantum states from those at the ?/'(3770). It is on this scan data that 
this thesis is based. 
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3.1 Data and Monte Carlo Samples 

The scan run was designed to provide maximum information in the available 
running period of August-October, 2005. The objective at each energy point was 
a measurement of the cross sections for all accessible final states consisting of 
a pair of charmed mesons. At the highest energy the possibilities include all of 
the following: DD, D*D, D*D*, D^Ds, D^D^, and D*D*, where the first four 
include both charged and neutral states. The original plan for the scan included 
ten energy points, with an integrated luminosity target for each point of ~ 5pb"^. 
It was recognized that specific energies might reveal themselves as unpromising 
with less than this luminosity. The plan was therefore designed to be flexible, with 
the option of adding or revisiting points to the extent that the data suggested 
and time allowed. In the end two points were added to the original list, including 
one at 4260 MeV to investigate the y(4260). The center-of-mass energies and 
integrated luminosities for the twelve scan points are given in Table 13.11 and Fig. 
13.11 An additional point was added as a result of the scan and corresponds to 
the location that maximizes the Dg yield. This point, 4170 MeV, is added to this 
analysis and its larger data sample was essential in understanding the nature of 
charm production throughtout this energy region. 

Numerous Monte Carlo (MC) samples have been generated in the develop- 
ment of the procedures and the determination of efficiencies and backgrounds 
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Table 3.1: Center-of-mass energies and integrated luminosities for all data points 
in the CLEO-c scan. 



^ (MeV) 


jCdt (nb-i) 


3970 


3854.30 ± 8.48 


3990 


3356. 15 ± 7.73 


1010 


5625.(35 ± 10.20 


4015 


1470.35 ± 5.32 


4030 


3005.55 ± 7.48 


4060 


3285.65 ± 7.79 


4120 


2759.20 ±7.29 


4140 


4871.85 ± 9.77 


4160 


10155.40 ± 13.96 


4160 


178942.15 ± 59.37 


4180 


5666.90 ± 10.51 


4200 


2809.90 ±7.59 


4260 


13107.60 ± 16.45 
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Figure 3.1: Integrated luminosity as a function of energy: The two largest inte- 
grated luminosity points correspond to 4160 and 4260 MeV. A point correspond- 
ing to 4170 MeV is left off, but has a luminosity of about 178 pb^^ 



for this analysis. The goal in specifying these MC samples was not to repro- 
duce reality precisely, but to include all relevant final states in sufficient quantity 
to develop selection criteria and assess the potential for cross-feed backgrounds. 
Thirteen 100 pb~^ samples (one for each center-of-mass energy), both signal and 
continuum, were generated. At each energy the signal MC sample includes all 
kinematically allowed D^^^D^^!^ final states. The continuum samples included only 
the uds background. All MC samples for the scan were generated on a system 
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of ~ 50 computers operated by the Minnesota CLEO-c group ("MN MC farm"). 
A breakdown of the signal MC samples is shown in Table 13.11 Additional sam- 
ples, described later, were subsequently produced to aid in assessing the potential 
contributions of "multi-body" production. 

Table 3.2: Composition of the 100 pb^^ signal MC samples. The assumed total 
charm cross section was 11 nb for Ecm = 4015 — 4200 Mev and 5 nb for Ecm = 
3970 - 4015 and 4260 MeV. 



Event 


3970 


4015 


4020-4080 


4080-4200 


4260 




11.6% 


6.9% 


6.9% 


6.6% 


6.3% 


D+D- 


11.6% 


6.9% 


6.9% 


6.6% 


6.3% 




35% 


13.7% 


13.7% 


13.2% 


12.6% 


D*+D- 


35% 


13.7% 


13.7% 


13.2% 


12.6% 






54.9% 


27.4% 


26.4% 


25.1% 


D*+D*- 






27.4% 


26.4% 


25.1% 




7% 


4% 


4% 


2.5% 


1% 










5% 


2% 


Dl+D*- 










9% 



3.2 Decay Modes and Reconstruction 

Eight Dg modes are used to measure the Dg production at each energy point. 
These modes are listed in Table 13.31 In addition to these Dg decays, several 
and decay modes, given in Table [3741 were used to investigate the amount of 
DD, D*D, and D*D* produced at each energy point. 

In selecting these decays the CLEO-c standard DTAG [H] code was used with 
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Table 3.3: The Dg branching fractions including the updated branching fractions, 
in 10-2. 



Modes 


Branching Fraction 


07r+, 10 MeV cut on the Invariant K+K' Mass [44J 


1.98 ±0.15 


K*<^K+,K*'^ K-TT- ig 


2.2 ±0.6 


7]7r+,r] ^ 77 [21111] 


0.58 ±0.07 


rip^, T] — i> 77, p"*" 7r+7r° [2j 


4.3 ± 1.2 


rj TT"*", ?7 ^ 7r+7r-?7, ^ 77 [21 IHj 


0.7 ±0.1 


7] p^,v[ — s> 'K~^7T~ri,ri 77, — * tt+tt" [2] 


1.8 ±0.5 


^ p+ ^ 7r+7r° [2] 


3.4 ±1.2 


KsK+,Ks TT+TT- [21111 


1.0 ±0.07 



the following modifications to the usual criteria: 

• The ^ requirements for charged pions and kaons were relaxed from 3cr to 
3.5a. 

• The Kg mass requirement was tightened from 4.5cr to 3o". 

• The /S.E cut for tag selection was relaxed from |Ai?| < 0.1 GeV to \^E\ < 
0.5 GeV. 

• The Mbc cut for tag selection was relaxed from 1.83 < Mbc < 2.0 GeV to 
1.7 < Mbc < 2.14 GeV. 

In addition to the above, the following cuts on intermediate-particle masses rela- 
tive to nominal values were applied for the modes: 

• \Mp - M^^^l < 150 MeV 
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Table 3.4: Decay modes used to select and D'^ in the scan data. The branching 
fractions are from a published CLEO-c analysis [121 HB] . 





JJl CLliL-liliiti ± IdL'ljl'^li 


rlpppiAA rnoHp 






3.91 ±0.12% 




14.94 ± 0.56% 




8.29 ± 0.36% 


D'^ decay mode 






9.52 ± 0.37% 




6.04 ±0.28% 




1.55 ±0.08% 




7.17 ±0.43% 




3.2 ±0.19% 




0.97 ±0.06% 



• \Mk*-M^P^\ < 75 MeV 

• \M^ - M^^^l < 10 MeV 

• \M.r,'-M^^^\ < 10 MeV 

MC invariant-mass plots for the selection of intermediate states used in the re- 
construction of Dg are shown in Fig. 13. 2[ 



3.3 Event Selection and Kinematics 

We begin by assuming that charm production in the threshold region is dominated 
by final states with two charmed mesons and no other particles: e^e~ 
For these events, the energy and momentum of the primary charmed mesons are 
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Figure 3.2: Invariant- mass distributions for the intermediate states involved in Ds 
reconstruction. The plots are for 4160 MeV MC: a) ^ K'^K~, b) rj' Tr'^n'r), 
c) K* K-n+, and d) p+ 7r+7r°. 



well defined. In general, for e'^e — > XY in the center-of-mass frame we have 



Ex — 



s + Ml-M'^ 



2^s 



(3.1) 



and 



IP;.! = = ipl = ^ El - M 



(3.2) 



where s — A:E'^^^^. The energy and momentum of one of the two charmed mesons 

is therefore sufficient to assign an event to one of the possible two-body processes. 
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In practice, we use familiar forms of these variables for this classification: the 



candidate's beam-constrained mass (Mbc = y E^e^m " l-^P ) energy deficit 

relative to the beam {AE = -Ebeam — Ed). 

As the center-of-mass energy increases above DD and D^D~ thresholds, it 
becomes possible to produce the "starred" states, D*^ , D*^ and -D*"*". These are 
not fully reconstructed in this analysis, since momenta and energies are suffi- 
cient to identify the origin of the reconstructed -D(s)- Reconstructed D and Dg 
candidates from these sources do not have a well-defined momentum, since they 
are daughters of starred parents and exhibit Doppler broadening. This Doppler 
broadening manifests itself through smeared distributions in both energy and mo- 
mentum. Some properties of the intermediate starred states are summarized in 
Table 



Table 3.5: and D*+ decay modes [2] 



Modes 


Branching Fraction 


D*^ decays mode 






61.9 ±2.9% 




38.1 ±2.9% 


D*^ decays mode 






67.7 ±0.5% 




30.7 ±0.5% 


L)+7 


1.6 ± 0.4% 


Z^*+ decays mode 






94.2 ± 2.5% 




5.8 ±2.5% 
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To illustrate the separation of events we show the momentum spectrum of 
candidates within 15 MeV of the nominal mass in Fig. 13.31 for the center-of-mass 
energy 4160 MeV. The top plot in Fig. 13.31 is from the MC described in Sect. 13.11 
and the bottom plot is from the 10.16 pb^^ of data. There are three distinct 
concentrations of entries near 0.95, 0.73 and 0.5 GeV/c, corresponding to DD, 
D*D, and D*D* production, respectively. Similar distributions for Dg candidates 
are shown in Fig. 13. 4[ In this case there are two distinct peaks in the MC at 0.675 
and 0.4 GeV/c, corresponding to DfD~ and D*^Dj, respectively, which are the 
only two accessible final states at this energy. Only the peak at 0.4 GeV/c is 
visible in data, demonstrating that the cross section for DfD^ at this energy is 
consistent with zero. At center-of-mass energies above D*^D*^ threshold, such 
as 4260 MeV (not shown), there are three accessible final states, so three peaks 
in momentum are possible. 

For this analysis we need to separate clearly and measure the cross sections for 

the nine event types: D^D^, D*^D^, D*^D*^, D+D', D*+D-, D*+D*-, DfD~, 

D*^D^, and D*^D*^. We want to choose variables that are as orthogonal as 

possible for the purpose of this separation. For D'^D^, D^D^ and D^D^, the 

variables used were the D candidate's energy {AE) and its momentum (in the 

form of Mbc). The separation of Z^-meson candidates at 4160 MeV into the 

possible event types is illustrated in the AE vs. Mbc plot in Fig. 13. 5[ The 

corresponding plot for Dg is given in Fig. 13. 6[ The quantitative task is to count the 
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Figure 3.3: The momentum spectrum at 4160 MeV for D° — > ir~7r"'" candidates 
within 15 MeV of the nominal mass. The top plot is from the MC described in 
the text and the bottom plot is from the 10.16 pb~^ of data at this energy. The 
three distinct concentrations of entries near 0.95, 0.73 and 0.5 GeV/c correspond 
to DD, D*D, and D*D* production, respectively. 



65 




P„ (GeV/c) 



Figure 3.4: The momentum spectrum at 4160 MeV for Df — > (/)7r+ candidates 
within 15 MeV of the nominal mass. The top plot is from the generic MC dis- 
cussed in the text and the bottom plot is from the 10.16 pb~^ of data at this 
energy. The peak at ~ 0.675 McV/c is missing from the data, giving advanced 
notice that the cross section for DfDj at this energy is consistent with zero. 
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events and determine the cross section for each event category while controlhng 
contributions from backgrounds and cross-feed from other two-charm final states. 



+ ContinuuTi 

- d:d: 




1.9S 



Figure 3.5: AE vs. Mtc for D° R-tt^ at 4160 MeV using the 100 pb^^ MC 
sample. The plot illustrates the clear separation that is achieved for this choice 
of variables. 



To obtain the number of signal candidates for each event type, a signal region 

must first be defined. The signal region in Mbc for D^D^, D^D~ and D^D~ is 

±9 MeV around the respective particle masses. For the other event types {D*^D^, 

D*^D*^, D*+D-, D*+D*-, DI+D-, and D*+D*-) the requirement is ±15 MeV in 

invariant mass. To estimate combinatoric and other backgrounds, a sideband is 
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Figure 3.6: AE vs. Mbc for (f)7i+ at 4160 MeV using the 100 pb"^ MC 

sample. The plot illustrates the clear separation that is achieved by the choice of 
variables not just for Ds but also for since the branching ratio for (j)7r~^ 
is nonzero. 
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defined on either side of tlie signal region. In every case tlie sidebands are spaced 
from the nominal particle mass by 5 cr. The sizes vary from mode to mode because 
of differing resolutions and the need to exclude potential peaking contributions 
(such as decay modes that are common between and D^). These sidebands 
are all chosen to be significantly larger than the signal region to minimize the 
statistical uncertainty of the background subtraction. MC is used to determine 
the sideband normalization, which is defined as the total sideband yield divided 
by the MC-tagged background contribution in the signal region. In almost all 
cases the normalization given by MC is consistent with the ratio of the sizes of 
the signal and sideband regions. The background procedure is illustrated with 
one example for Mbc {D^ K^-k^tx^ in D^D^) and one for invariant mass 
{D+ K-71+7T+ in D*+D*-) in Figs. O and ESI respectively. 

For the other event types (i.e. those involving one or two "starred" charmed 
mesons) the variables used for the separation were Mbc (candidate momentum) 
and invariant mass, which is a combination of momentum and energy. Regardless 
of the origin of a D candidate, the invariant mass peaks at the D mass. Since 
invariant mass does not differentiate between event types, Mbc provides all of the 
event-type separation for these events. The separation can be seen in Figs. 13.91 

The difference in variable choices between unstarred and starred events is a 

matter of convenience. We use the same procedure that has been used with great 
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Figure 3.7: Plot of Mbc for K-TT+n+ in D+D' events in MC at 4160 MeV. 

The red indicates the signal region and the blue the sideband. 
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Figure 3.8: Plot of the invariant mass for K~tt^7c~^ in D*~^D*~ events in 

MC at 4160 MeV. The red indicates the signal region and the blue the sideband. 
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Figure 3.9: Invariant mass vs. Mbc for K^n^ at 4160 MeV using the 100 

pb^^ MC sample. The pfot illustrates the clear separation that is achieved with 
this choice of variables. 
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Figure 3.10: Invariant mass vs. Mbc for — >■ ^tt"*" at 4160 MeV using the 100 
pb~^ MC sample. The pfot illustrates the clear separation that is achieved with 
this choice of variables, not just for Ds but also for D~^, since the branching ratio 
for 07r+ is nonzero. 
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success in the CLEO-c analysis of ^^(3770) DD in Refs. |l5l HB]. Mbc has 
a practical advantage over the raw momentum in that Mbc changes much more 
slowly with beam energy. For a given center-of-mass energy, the expected values 
of Mbc and AE in e+e^ — > XY are given by 



As the center-of-mass energy increases from 3970 to 4180 MeV, Mbc for the D in 
DD* changes by ~ 0.01%, while the momentum changes by ~ 80%. 

As discussed above, the method for D^D^, D^D^ and D^D^ is to cut on 
{/S.E < 15 MeV) and use the Mbc distribution to determine the yield (Fig. 13. 7p . 
For all other event types the method is to cut on Mbc and use the invariant-mass 
distribution to determine the yield (Fig. 13.81) . The cut on Mbc was determined 
by kinematics, since Mbc is center-of-mass energy dependent in addition to being 
dependent on the nature of the decay of the starred state. In order to choose the 
cut range we assumed that D*^ and D*+ decay 100% of the time by jD^ 
and D*~^ —>■ respectively, since the pion transitions will fall in the cut 

window for the 7 decay. We assumed that D*^ decays only to D^n^ and D^n^, 
since the branching ratio for D*^ D^'j is only 1.6%. The equation that 




(3.3) 



and 



AE 



lMl-_ 
2 



M2 



(3.4) 
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determines the Mbc cut window for D*D and D*D* with D* — > DX is 

M^t = ^EL.m -Ml- {lEo ± hPD)\ (3.5) 

with 7 = /57 = ^ and = ^"*2aS"^^ ' ^^^^^ Mfl,^^ and are 

the mass, energy, and momentum of the daughter D in the rest frame of the D*, 
and and P^* are determined by Eqs. 13.11 and 13. 2^ respectively. For D*^ and 
D*"*", Mx = M7 = 0, since only the 7 decay is considered in the calculation and 
for Mx = M^o = 135 MeV. After calculating the maximum and minimum 
values using the above equation and assumptions, we expand the interval by 
5 MeV on each end to account for resolution effects. 

An important point is that there is a certain energy at which an overlap 
between the D^D*^ and D*^£)*^ will occur. This is only a problem for D^D*^ and 
D*^D*^, and not D^D*^ and D*^D*^, because only the vr decay was considered. 
This means these event types will have a significantly smaller cut window in Mbc 
as compared to neutral D events. If this occurs, the high Mbc cut of D^D*^ will 
be set equal to the low Mbc cut of D*^D*^ so that events cannot pass both event- 
type cuts and be double-counted. This is done because r)*^D*^ will have twice as 
many D mesons populating this area, assuming equal rates, therefore there will 
be less contamination from D^D*^ in D*o/)*o ]yy this method as compared to its 
inverse. 
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Table 3.6: Selection criteria, in units of Mbc, for measuring the cross sections for 
D*D and D*D*. A 15 MeV cut on was made when selecting DD. 



E,ra (McV) 






D+D*- (^) 


D*+D*- (^) 


3970 


(1.902,1.972) 


— 


(1.933,1.958) 


— 


3990 


(1.901,1.976) 


— 


(1.933,1.961) 


— 


4010 


(1.899, 1.98) 




(1.934, 1.963) 




4015 


(1.899,1.981) 


(1.994,2.01) 


(1.934, 1.963) 




4030 


(1.898, 1.984) 


(1.987,2.02) 


(1.934,1.965) 


(2.003,2.018) 


4060 


(1.896,1.98) 


(1.98,2.03) 


(1.935,1.968) 


(2.002,2.022) 


4120 


(1.893, 1.974) 


(1.974,2.044) 


(1.935, 1.974) 


(2.003,2.029) 


4140 


(1.892,1.972) 


(1.972,2.048) 


(1.935,1.976) 


(2.004,2.032) 


4160 


(1.891,1.97) 


(1.97,2.052) 


(1.935,1.978) 


(2.004,2.034) 


4170 


(1.89,1.97) 


(1.97,2.054) 


(1.935,1.979) 


(2.005,2.035) 


4180 


(1.89,1.969) 


(1.969,2.056) 


(1.935,1.98) 


(2.005,2.036) 


4200 


(1.889,1.968) 


(1.968,2.06) 


(1.935,1.982) 


(2.006,2.038) 


4260 


(1.887,1.965) 


(1.965,2.071) 


(1.935,1.988) 


(2.008,2.044) 



Table 3.7: Selection criteria, in units of Mbc, for measuring the cross section for 
D*Ds and D*D*. A 15 McV cut on AE was made when selecting DsDg. 



Ecm MeV 




7~)*+ 7~)*~ GeV 


3970 






3990 






4010 






4015 






4030 






4060 






4120 


(2.015,2.06) 




4140 


(2.011,2.067) 




4160 


(2.009,2.071) 




4170 


(2.008,2.073) 




4180 


(2.007,2.076) 




4200 


(2.005,2.08) 




4260 


(2.001,2.087) 


(2.087,2.132) 
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Chapter 4 



Cross Section Calculation 



4.1 Determination of Cross Sections 
4.1.1 Exclusive Cross Sections 

The cross section for from any of the nine possible event types can be com- 
puted with the following equation: 

AT (signal) 

a^(.)5(.)(A«)) = ^^, (4.1) 

where iV(signal) is the number of signal events, B is the branching ratio for the 
particular D decay being used (Table E31 13 IB] and Table 1^:1 [151115] ) . C is the 
integrated luminosity, and e is the detection efficiency (determined with MC). 

iV(signal) is obtained by counting candidates in the signal region and subtract- 
ing sideband-estimated backgrounds with normalizations determined by MC. This 
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procedure and the definitions of tlie signal and sideband regions for each mode 
appear in Sect. 13.31 

The cross sections for the DD, D*D and D*D* production modes are calcu- 
lated as follows: 

(t{D*^D*^) = ^{ao^D* P°) - 2i3(D*+ ^ D\+)a{D*+D*-)), (4.3) 

a{D*^D'^) = 1((Tb,5(Z}0) - B{D*+ D\+)a{D*+ D-)), (4.5) 

-(D-D-) = (4.6) 

and 

aiDD) = (4.7) 

where (D"*") is the cross section of produced in D*D* events, a£,*i)*{D^) 

is the cross section of produced in D*D* events, a£)*£){D^) is the cross section 
of produced in D*D events, aD*jj{D^) is the cross section of produced in 
D*D events, £,(/)+) is the cross section of produced in DD events, and 
aj;)f){D^) is the cross section of D° produced in DD events. 

To determine cr{D^) and a{D~^), weighted averages of the three D^ modes and 
five D^ modes are calculated, with weights defined as — . 

For DgDs, D*Ds and D*D*, a weighted sum technique is used to combine the 

eight Dg decay modes and obtain the cross sections. The weights that minimize 
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the error are given by the following: 



(A i 



where fi = jf, Ni is the yield, ai is the uncertainty on the yield for mode i. 



Z^i=0 ffN, 

It is perhaps counterintuitive that for modes with equal precision the weight 
is inversely proportional to the yield, thereby suppressing the weight of modes 
with high yields. This feature guarantees that the weighting of different modes 
in the cross section measurement is determined by the precision of the yield 
measurement rather than its magnitude. (The conclusion can be easily verified 
by considering a "toy" example of two modes, each with 10% precision and yields 
of 100 and 1000, respectively. Eq. 14.81 properly assigns roughly equal weighting 
and achieves an uncertainty of ~ 7%, rather than the 10% that would follow if 
the "bigger" mode were allowed to dominate.) 

The weights were determined by MC for each of the eight modes for each of the 
three possible event types, DgDs, D*Ds, D*D*. The weights were then averaged 
across these possible event types and used in calculating the cross sections. The 
equation for determining either the DgDs, D*Ds, D*D* cross sections is as follows: 



where e is determined with MC by the following: 



(4.9) 



e = ^^y^'- (4.10) 

-'^generated 

79 



where the 8 in the summation refers to the eight Dg modes that were used during 



the scan. 



4.1.2 Efficiencies for Exclusive Selection 

Efficiencies for exclusive selection of all accessible event types and center-of-mass 
energies were determined by analyzing the MC samples described in Sect. 13.11 
The efficiency for detecting a particular decay mode is defined as follows: 



-^detected obtained by applying the same selection and background-correction 
procedures to MC as are used to determine the signal yields in data. The effi- 
ciencies for each of the twelve energies and event types are listed in Tables 14.11 
14.21 (-D"*") and 18.51 {Df - note that this multi-page table appears on P.193 
after the References). 

4.1.3 Cross Section Results 

Following the procedure laid out in the previous sections, the production cross 
sections were measured at thirteen center-of-mass energies. The results are shown 



e = 




(4.11) 



where the error is determined by binomial statistics: 




(4.12) 



in Tables ES\ and EM and Figs. O and S^l 
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Table 4.1: Efficiencies (units of 10 ^) at each scan-energy point for selection of 

— > K~7r^ decays in the three exchish'e event tY]:)es. 



Ecm (McV) 


DD 


D*D 


D*D* 


3970 


62.42 ± 0.72 


61.78 ±0.36 


- ± - 


3990 


63.04 ±0.72 


62.37 ±0.36 


- ± - 


4010 


61.75 ±0.73 


62.83 ±0.36 


- ± - 


4015 


62.29 ±0.59 


63.31 ± 0.36 


64.23 ±0.24 


4030 


63.75 ± 0.63 


62.68 ±0.39 


64.60 ± 0.24 


4060 


60.17 ±0.64 


61.24 ±0.39 


64.42 ± 0.24 


4120 


58.97 ±0.66 


59.75 ± 0.40 


64.23 ±0.25 


4140 


59.49 ± 0.47 


57.90 ±0.29 


64.84 ±0.17 


4160 


57.94 ±0.66 


58.12 ±0.40 


64.80 ±0.25 


4170 


58.12 ±0.66 


56.53 ±0.40 


65.23 ±0.25 


4180 


59.30 ±0.66 


54.72 ± 0.41 


66.00 ±0.25 


4200 


59.06 ±0.66 


52.52 ±0.41 


66.88 ±0.24 


4260 


56.29 ± 1.01 


41.36 ±0.61 


70.02 ± 0.36 



Table 4.2: Efficiencies (units of 10 ^) at each scan-energy point for selection of 

— >■ K^Ti^n^ decays in the three exclusive event types. 



Ee„ MeV 


DD 


D*D 


D*D* 


3970 


56.07 ±0.49 


52.82 ± 0.35 


- ± - 


3990 


55.33 ±0.48 


52.78 ±0.35 


-±- 


1010 


55.59 ± 0. 19 


53.13 ±0.35 


- ± - 


4015 


53.60 ±0.39 


53.37 ±0.34 


- ± - 


4030 


56.79 ± 0.43 


52.33 ±0.37 


53.60 ±0.38 


4060 


53.93 ±0.43 


52.53 ±0.37 


53.12 ±0.38 


4120 


54.42 ± 0.44 


52.55 ±0.38 


51.64 ±0.38 


4140 


53.83 ±0.31 


51.18 ±0.27 


52.47 ±0.27 


4160 


54.12 ±0.44 


51.10 ±0.38 


50.98 ±0.38 


4170 


54.15 ±0.44 


51.15 ±0.38 


51.31 ±0.38 


4180 


54.25 ± 0.44 


51.39 ±0.38 


52.46 ±0.38 


4200 


53.38 ± 0.44 


51.37 ±0.38 


52.47 ±0.38 


4260 


54.51 ±0.66 


49.15 ±0.58 


51.29 ±0.59 
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Figure 4.1: Observed cross sections for e'^e" DD, D*D and D*D* as a function 
of center-of-mass energy. Errors are just statistical. 
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Figure 4.2: Observed cross sections for e+e DgDs, D*Ds and D*D* as a 
function of center-of-mass energy from the scan. Errors are just statistical. 
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These cross sections were then used in Eqs. l4.2ti4.7l to obtain the results shown in 
Tables EM and EM 

The weights and efficiencies are shown in Tables S]3] and |431 By summing the 
individual exclusive cross sections, one arrives at the total observed charm cross 
section which is shown, along with the exclusive cross sections, in Table 18.81 and 

El 



Table 4.3: The efficiencies in units of 10 ^ for detecting D^Ds, D*Ds, and D*D* 
at each energy point in the CLEO-c scan. 



Ecm MeV 






<d*+d:-) 


3970 


1.14 ±0.02 


- ± - 


- ± - 


3990 


1.11 ±0.02 


- ± - 


- ± - 


4010 


1.11 ±0.02 


- ± - 


- ± - 


4015 


1.13 ±0.02 


- ± - 


- ± - 


4030 


1.20 ±0.02 


- ± - 


- ± - 


4060 


1.09 ±0.02 


- ± - 


- ± - 


4100 


1.12 ±0.03 


1.00 ±0.02 


- ± - 


4120 


1.04 ±0.03 


1.04 ±0.02 


- ± - 


4140 


1.04 ±0.02 


1.06 ±0.02 


- ± - 


4160 


1.09 ±0.03 


1.06 ±0.02 


- ± - 


4170 


1.09 ±0.03 


1.06 ±0.02 


- ± - 


4180 


1.09 ±0.03 


1.10 ±0.03 


- ± - 


4200 


1.09 ±0.03 


1.07 ±0.03 


- ± - 


4260 


1.09 ±0.03 


1.10 ±0.03 


1.05 ±0.03 



84 



Table 4.4: For DgDs, D*Ds, and D*D*, a weighted sum technique that minimize 
the error is used to determine the cross sections. The weights for each mode 
are shown below. The mode with the largest weight is (pn'^, which is by far the 
cleanest of all the modes. 



Mode 


Weight 




0.14 




0.06 




0.23 


K*K+ 


0.12 




0.05 


7] 7T^ 


0.16 


1 + 

V P 


0.10 




0.13 



4.1.4 Two Other Methods for Measuring the Total Charm 

Cross Section 
Inclusive D Method 

In addition to measuring the separate cross sections for all expected charm event 
types, one can perform inclusive measurements to obtain the total observed charm 
cross section. 

As for the exclusive measurements, the efficiencies for inclusively selecting 
events with charmed mesons are determined with the MC samples described in 
Sect. 13.11 They are given in Table 1^31 

For and D'^, the event-type requirements on \AE\ and Mbc are lifted and 
the invariant mass is used to extract the yields. The inclusive K^tt^ and 
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Table 4.5: Efficiencies at each scan-energy point for inclusive selection of D' 
K-TT+, D+ K-Tr+'K+, and Df -> K~K+tt+ (units of IQ-^). 



F, MpV 






t\Ug j\ j\ n ) 


ov t u 


UO.ZiU ZL U.OZi 


p:o qy 1 n 90 


p;4 Qi _|_ n qn 




(^1 qq _L n oo 

Ul. c/c/ ZL U.OO 


c^o oq _L q 90 


c:q n4 + n QD 
00. zn u.c/u 


4m n 


62 52 + 32 


53 57 ± 28 


53 62 ± 90 


4015 


63.79 ±0.19 


51.91 ±0.26 


52.10 ±0.74 


4030 


64.52 ±0.20 


53.47 ±0.23 


52.87 ±0.81 


4060 


63.24 ±0.20 


51.95 ±0.23 


48.83 ±0.80 


4120 


62.49 ±0.20 


52.33 ±0.23 


50.90 ±0.59 


4140 


62.61 ±0.14 


51.79 ±0.16 


49.84 ±0.42 


4160 


62.54 ±0.20 


51.17±0.23 


52.75 ±0.59 


4170 


62.26 ±0.20 


51.39 ±0.23 


52.93 ±0.59 


4180 


62.32 ±0.20 


51.95 ±0.23 


53.67 ±0.59 


4200 


63.13 ±0.20 


52.00 ±0.23 


51.63 ±0.59 


4260 


60.82 ±0.31 


51.35 ±0.35 


53.74 ±0.69 



— i> K^7r+7r+ invariant-mass spectra are shown for the 4160 MeV data sample 
in Figs. 14.31 and 14.41 For Dg the event-type requirements are preserved because 
of the need to suppress the background for the high-yield mode K~K'^n^. At 
energies above 4100 MeV, for all candidates that pass the selection requirements 
for DfD^, D*^D^ and D*^D*' (the last only for 4260 MeV), the invariant mass 
is used to determine the inclusive yield. At energies below 4100 MeV, Mbc is used 
to determine the yield since these energies are below D*^Dj and D*^D*^ thresh- 
olds. The inclusive Df K^K^ir^ invariant mass in 4160 MeV data is shown 
in Fig. 14.51 Each histogram is fitted to a function that includes a Gaussian signal 
and an appropriate background function. For K^n^, K^n^n^ and 
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— s> K^K^TT^ above 4100 MeV, the background function is a second-order 
polynomial. For — > K^K^tc^ below 4100 MeV, the background function is 
an Argus function [38] . The results of the fits are shown in Table 14.61 




Figure 4.3: The invariant mass of ^ K 7r+ in data at Ecm = 4160 MeV. 



The observed inclusive cross section can be determined as follows: 

/ J- N A^(signal) 

a(e+e- ^ D(,)X) = (4.13) 

where A^(signal) is the signal yield from the inclusive fit (Table 14. 6p . e is the 

efficiency for detecting the signal (Table 14.51) . B is the branching ratio for the 

particular D decay in question, and C is the integrated luminosity. The cross 
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Figure 4.4: The invariant mass of — > K tt'^tt^ in data at Ecm = 4160 MeV. 
The peak at ~ 2 GeV is fully reconstructed D*+ L>°7r+, D° K-n+. 
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Figure 4.5: The invariant mass of — > K^K tt"'" in data at Ecm = 4160 MeV. 
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Table 4.6: Yields from inclusive fits for D^, and 



F, MpV 


\^ I 






'?Q7n 

ov t u 




«Q« An _|_ QO QQ 


QC, 00 _(_ 1 4 /I7 




«/lQ qi _|_ on Qc; 


uuu.tji zn ou.ou 


o^.u^ zn 1 .o\j 


zLDI n 






«q cfl + 1 1 7n 
UtJ. tju zn 1 1 . 1 u 




oyn 07 _)_ 01 Qc; 
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4030 


1192.53 ±37.29 


844.87 ±36.18 


43.42 ±9.14 


4060 


1191.66 ±37.68 


859.62 ±36.77 


0.00 ±0.00 


4120 


881.87 ±32.89 


576.56 ±30.62 


77.30 ± 14.50 


4140 


1512.51 ±43.45 


1006.17 ±42.40 


181.21 ±21.29 


4160 


3130.11 ±62.79 


2028.33 ±59.92 


448.12 ± 32.75 


4170 


54958 ± 262.702 


34050.8 ± 247.879 


8379.97 ± 146.051 


4180 


1668.05 ±46.10 


1088.37 ±44.12 


243.47 ±23.80 


4200 


727.11 ±30.61 


445.36 ±27.59 


97.80 ± 15.19 


4260 


1691.56 ± 51.13 


1370.82 ±52.60 


330.49 ± 34.93 



section times branching ratio for D^, and Dg and the cross sections determined 
from these with the branching ratios in Ref. [45] and Ref. [2] are shown in Tables 
14.71 and lTSl respectively. 

From the inclusive cross sections, the total charm cross section can be ob- 
tained. Since all D mesons are produced in pairs the total cross section for 
production of charm events is given by 



a[e e 



DDX) 



a^o + (J£,+ + a 



(4.14) 



at each energy point. These results are given in Table I^l9l 
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Table 4.7: The cross section, in nb, times branching ratio of D , and Df. 





a ■ B(D^ ^ K-7T+) 






3970 


0.244 ± 0.012 


311 ± 016 


017 ± 007 


3990 


309 ± 014 


342 ± 017 


019 ± 004 


4010 


0.364 ± 0.012 


0.417 ± 0.014 


0.023 ± 0.004 


4015 


0.395 ±0.023 


0.454 ± 0.030 


0.021 ±0.007 


4030 


0.615 ±0.019 


0.526 ±0.023 


0.027 ±0.006 


4060 


0.574 ±0.018 


0.504 ±0.022 


0.000 ± 0.000 


4120 


0.511 ±0.019 


0.399 ±0.021 


0.055 ±0.010 


4140 


0.496 ±0.014 


0.399 ±0.017 


0.075 ± 0.009 


4160 


0.493 ±0.010 


0.390 ±0.012 


0.084 ±0.006 


4170 


0.493 ± 0.002 


0.370 ±0.003 


0.088 ±0.002 


4180 


0.472 ±0.013 


0.370 ±0.015 


0.080 ±0.008 


4200 


0.410 ±0.017 


0.305 ±0.019 


0.067 ±0.010 


4260 


0.212 ±0.006 


0.204 ± 0.008 


0.048 ± 0.005 



Table 4.8: The inclusive cross section of and 



^cm MeV 


cr(£,o) nb 


a(D+) nb 


nb 


3970 


6.25 ± 0.30 


3.27 ±0.17 


0.31 ±0.13 


3990 


7.92 ±0.36 


3.59 ±0.18 


0.34 ±0.08 


4010 


9.31 ±0.30 


4.38 ±0.15 


0.41 ±0.07 


4015 


10.10 ±0.58 


4.77 ±0.31 


0.38 ±0.13 


4030 


15.73 ± 0.49 


5.52 ±0.24 


0.49 ±0.10 


4060 


14.67 ±0.46 


5.29 ±0.23 


-±- 


4120 


13.08 ±0.49 


4.19 ±0.22 


0.99 ±0.19 


4140 


12.68 ±0.36 


4.19 ±0.18 


1.34 ±0.16 


4160 


12.60 ±0.25 


4.10 ±0.12 


1.50 ±0.11 


4170 


12.51 ±0.06 


3.90 ±0.03 


1.59 ±0.03 


4180 


12.08 ±0.33 


3.88 ±0.16 


1.44 ±0.14 


4200 


10.48 ±0.44 


3.20 ±0.20 


1.21 ±0.19 


4260 


5.43 ±0.16 


2.14 ±0.08 


0.86 ±0.09 
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Table 4.9: The total charm cross section as obtained by the inclusive method 



E^rr, MeV 


(j(e+e" DDX) nb 


3970 


4.910 ± 0.183 


3990 


5.926 ± 0.219 


4010 


7 050 ± 180 


4015 


7.623 ± 0.360 


4030 


10.870 ± 0.294 


4060 


9.979 ±0.258 


4120 


9.132 ±0.281 


4140 


9.105 ±0.212 


4160 


9.103 ±0.147 


4170 


9.094 ±0.066 


4180 


8.701 ±0.193 


4200 


7.448 ± 0.250 


4260 


4.212 ±0.093 



Hadron-Counting Method 

Still another method and cross-check can be done. This method involves counting 
the multihadronic events in the data at all thirteen energy points and using data 
collected below cc threshold at = 3671 MeV to subtract uds continuum 
production. Except for one difference to be discussed later, this method is the 
same as that used by CLEO-c [19], [50], EI] to determine the cross section of e+e^ 
'ip{3770) — ^ hadrons at i?cm = 3770 MeV. The Standard Hadron (SHAD) cuts as 
discussed in Refs. [HI [51] are used. 

This method first starts by calculating the number of hadronic continuum 
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events at = 3671 MeV: 



A^(3671) = {Nohs-N,^{2S)-Nj/^-ej/^-Nee-eee-N^^-e^^-Nrr-err)/econt, (4.15) 

where the A^'s are the numbers of events of different types as determined from 
data or by calculating C ■ a, and the e's are the efficiencies for passing the SHAD 
cuts of the same event types. The hadronic efficiency for events from the uds 
continuum is econt- The quantity N^{2S) in Eg. I4.15l is the contribution due to the 
Breit-Wigner tail of ip{2S). It is estimated from the data collected at Ecm = 3671 
and 3686 MeV as follows: 

7V(7r+7r / + / j3686 

where N{tlj{2S)) 3686 is the number of hadronic events in 'ip{2S) decays at = 
3686 MeV. The values used for iV(7r+7r-/+/-)3686 and iV(7r+7r-/+/-)367i in Eg. 14151 
were obtained by CLEO-c f55l and give a scale factor of = = 

■J L J o N(7r+-!r l+l )3686 30462 

0.0073. The number of hadronic events in ipi^S) decays at Ecm = 3686 MeV is 
determined by 

iV(V^(2S))3686 = iV(3868)obs - S ■ iVobs(3671), (4.17) 

where the scale factor, S = ■ (f|||)^ = 0.139. In Eg. liTTl we neglect the 
contamination of ip{2S) in the off-resonance data and QED events, which is small 
compared to the large number of ip{2S) decays present in the sample. 
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After determining A^(3671)obs, it can be used in determining the number of 
hadronic events at each scan point as follows: 

N{X) = (A^obs(-^) - S ■ A^obs(3671) - N^(2S) ■ ev(25) - A^v(3770) ■ ev(3770) - 

Nj/^ ■ ej/^ - A^ee ■ Gee " ■ e^^, - N^r " e^^),(4.18) 

where X stands for the energy point in question (3970, 3990, 4010, 4015, 4030, 
4060, 4120, 4140, 4160, 4170, 4180, 4200, and 4260 MeV), and S is the scale 
factor given by ■ (^f^^)^. 

^ £3671 ^ X MeV 

In determining iV(3671) CLEO-c's method was followed exactly; that is the 
data collected at the ipi^S) resonance was used in determining N^(2S) in Eq- 
14.151 However, in regards to N{X) we used the calculated production cross 
section in determining the amount of ip{2S), in addition to the amount of J/ip 
and ?/'(3770) present at each energy point. The calculated production cross section 
is a convolution of a (5-function-approximated Breit-Wigner and an ISR kernel: 

a(e+e--7X) = ^^-/(.,x), (4.19) 

where the ISR kernel f(x,s) is defined in Eq. 28 of Ref. and reproduced here 
as Eq. 17.21 in Sect. [3, with x = ^ and Mx referring to the mass of the ip{2S), 

J/ip or '?/'(3770). The calculated production cross sections for ip{2S), J/ip, and 
V'(3770) are shown in Table H?TOl It should be noted that we could use the results 
from our data samples, the number of e^e~ 7^/^(25') — > •yn^n^ J/ip events at 
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each of the scan energies, to improve this result in the future. Also, in calculating 
N{X) it was assumed that effects due to interference between e"'"e~ 'V^(2S') — >• 
7* — > hadrons and the continuum were negligible at all the scan energies, and so 
were only included in the calculation of iV(3671) using the method described in 
Ref. [501 EI]- The numbers determined and used in this method are shown for 
all energy points in Tables 18.61 and 18. 7[ 



Table 4.10: Calculated production cross sections, in nb, for J/ip, i/j{2S), and 
'?/'(3770) following the procedure discussed in the text. 



E,^ (MeV) 




a(e+e- ^ 7^/'(25)) 


cT(e+e~ ^ 7^/'(3770)) 


3970 


0.70 


0.92 


0.13 


3990 


0.68 


0.85 


0.12 


4010 


0.66 


0.79 


0.11 


4015 


0.66 


0.78 


0.11 


4030 


0.64 


0.74 


0.10 


4060 


0.62 


0.68 


0.09 


4120 


0.57 


0.57 


0.07 


4140 


0.56 


0.54 


0.07 


4160 


0.54 


0.52 


0.06 


4170 


0.54 


0.51 


0.06 


4180 


0.53 


0.49 


0.06 


4200 


0.52 


0.47 


0.06 


4260 


0.48 


0.41 


0.05 



Wide-angle Bhabha events were generated with BHLUMI [52] using a cut off 
angle of 21.57°; /i-pairs and r-pairs were generated with FPAIR [53] and KORALB 
[51] , respectively. These calculations provide the production cross sections needed 
for Eqs. 14.181 and 14.181 and were used to produce MC samples that were used to 
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determine the SHAD selection efficiencies. The QED production cross sections are 
shown in Table 14.111 and their efficiencies are given in Table 14.121 



Table 4.11: Calculated QED production cross sections at twelve energy points. 



E,m (MeV) 


e+e-(ft™„ = 21.57°) (nb) 


fi'^fj, (nb) 


r+r (nb) 


3670 


448.2 


8.11 


2.1 


3970 


383.47 


6.99 


3.32 


3990 


379.59 


6.98 


3.38 


4010 


376.39 


6.83 


3.4 


4015 


374.01 


6.86 


3.4 


4030 


372.34 


6.81 


3.44 


4060 


366.5 


6.69 


3.45 


4120 


355.95 


6.5 


3.51 


4140 


352.72 


6.47 


3.53 


4160 


348.78 


6.4 


3.53 


4170 


346.56 


6.36 


3.54 


4180 


344.77 


6.33 


3.55 


4200 


344.71 


6.29 


3.56 


4260 


332.4 


6.17 


3.57 



Once the number of supposed pure charm decays N{X) has been obtained, 
the total charm cross section can be calculated at each energy point as follows: 

a,(X) = (4.20) 

The results using this hadr on- counting method are shown in Table 14.131 A com- 
parison of all methods is shown in Fig. 14.61 
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Table 4.12: Efficiencies (units of 10 ^) for various event types to pas SHAD liadronic 
event selection criteria. 





e+e 








V(25) 


qq 


Signal 


3670 


0.03 ±0.01 


0.12 ±0.02 


22.0 ±0.1 


41.9 ±0.6 


68.3 ±0.1 


59.0 ±0.2 


-±- 


3970 


0.02 ±0.01 


0.12 ±0.02 


22.8 ±0.8 


39.3 ±0.3 


66.3 ±0.3 


63.1 ±0.2 


80.9 ±0.1 


3990 


0.04 ±0.01 


0.12 ±0.02 


20.9 ±0.8 


39.0 ±0.3 


65.5 ±0.3 


62.8 ±0.2 


80.9 ±0.1 


4010 


0.03 ±0.01 


0.12 ±0.02 


21.6 ±0.8 


38.5 ±0.3 


65.9 ±0.3 


63.3 ±0.2 


81.0 ±0.1 


4015 


0.03 ±0.01 


0.12 ±0.02 


22.8 ±0.8 


38.2 ±0.3 


65.6 ±0.3 


63.3 ±0.2 


82.0 ±0.1 


4030 


0.05 ±0.01 


0.12 ±0.02 


21.8 ±0.8 


38.1 ±0.3 


65.7 ±0.3 


63.6 ±0.2 


81.9± 0.1 


4060 


0.02 ±0.01 


0.12 ± 0.02 


22.8 ±0.8 


39.2 ±0.3 


65.7 ±0.3 


63.7 ±0.2 


82.7 ±0.1 


4120 


0.05 ± 0.01 


0.12 ±0.02 


22.8 ±0.8 


35.1 ±0.3 


65.3 ±0.3 


64.1 ±0.2 


82.6 ±0.1 


4140 


0.04 ±0.01 


0.12 ±0.02 


20.9 ±0.8 


35.0 ±0.3 


64.8 ±0.3 


64.5 ±0.2 


82.8 ±0.1 


4160 


0.04 ±0.01 


0.12 ±0.02 


21.7±0.8 


34.2 ±0.3 


65.7 ±0.3 


64.8 ±0.2 


82.8 ±0.1 


4170 


0.04 ±0.01 


0.12 ±0.02 


21.9 ±0.8 


34.3 ±0.3 


65.4 ±0.3 


64.9 ±0.2 


82.6 ±0.1 


4180 


0.04 ±0.01 


0.12 ±0.02 


23.6 ±0.8 


33.6 ±0.3 


65.0 ±0.3 


65.1 ±0.2 


82.7±0.1 


4200 


0.03 ±0.01 


0.12 ±0.02 


21.6 ±0.8 


32.2 ±0.3 


64.5 ±0.3 


65.1 ±0.2 


82.8 ±0.1 


4260 


0.05 ± 0.01 


0.12 ±0.02 


21.9 ±0.8 


28.8 ±0.3 


64.7 ±0.3 


65.7 ±0.2 


82.7 ±0.1 



Table 4.13: The total charm cross section as obtained by the Hadron Counting 
Method. Only statistical errors are shown. 



Ecm MeV 


(T(e+e — >• Charm) nb 


3970 


4.91 ±0.13 


3990 


5.87 ±0.14 


4010 


7.21 ±0.12 


4015 


7.88 ±0.18 


4030 


11.30 ±0.15 


4060 


9.98 ±0.14 


4120 


9.43 ±0.15 


4140 


9.58 ±0.12 


4160 


9.62 ±0.09 


4170 


9.44 ±0.09 


4180 


9.07 ±0.12 


4200 


8.37 ±0.14 


4260 


4.34 ± 0.08 



97 



12 



10 



S 6 



$.95 



; 


: 1 1 

i 

i 


* A 

• • ▲ 


4 






i 


i 

± 




■ 






* 




■ 

■ 


Hadron Counting 






i 


■ 


D Meson Counting Exclusive 




■ 


• 

1 1 1 1 1 1 1 1 1 


D IVIeson Counting Inclusive 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 


, 1 , 



4- 



4.05 ^4.1 „4.15 4.2 

Ecm (GeV) 



4.25 



Figure 4.6: Plot of the total charm cross section as calculated in each of the 
three methods described in the text. The reason for the discrepancy between the 
inclusive and exclusive method is due to the presence of multi-body production, 
which is described in the text. Only statistical errors are shown. 
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Chapter 5 



Momentum Spectrum Analysis 



5.1 Multi-Body, Initial-State Radiation, and 
Momentum-Spectrum Fits 

Throughout this analysis we have assumed that all charm production is through 
two-body events, that is e^e" — > With sufficient energy, however, there 

is no reason that final states like e+e" DD^*^n or any other energetically 
allowed combination with extra pions should not be produced. Prom here on, 
these types of events are referred to as the multi-body production or just multi- 
body. For production of non-strange D states, there is no a priori expectation 
for the amount of multi-body. For Ds we expect it to be small, if not zero, since 
DgDsTT violates isospin. In both cases it is appropriate to examine our data for 
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evidence of such multi-body processes. 

The first step is to determine whether multi-body exists. Assuming that it 
does, we then need to develop and apply procedures to determine its composition: 
DDti, D*Dn, or D*DTr7r. 

To determine if multi-body exists we can apply tests of consistency between 
our measurements and the expectations for pure D|*jz)|*j states. One observable 
is the ratio of as a function of energy, which is shown in Fig. 15.11 The 
bold horizontal lines are the predicted ratios as determined by the decays of the 
D* using the information in Table 13.51 It is evident that the observed ratio 
deviates from that expected for D*D* events. The only candidate explanation 
for this observation is multi-body production. The cut window in Mbc is different 
for neutral and charged D*. The charged window is quite a bit smaller than 
the neutral, since the 7 decay is excluded when determining the cut window. 
Therefore, the cuts can select different amounts of the multi-body, which leads 
to a result which is not consistent with that expected based on known branching 
fractions. The fact that the D*D events have the correct ratio will help in pin- 
pointing the make-up of the background. 

The ratio of gives a hint that there might be multi-body background 

present at the energies of interest, but it is far from conclusive. Another indication 

that points to a multi-body background is the noticeable difference in the total 

charm cross section between the exclusive and inclusive methods. The difference 
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Figure 5.1: The measured ratio of D'^ to as a function of energy. 
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between these two methods as a function of energy can be seen in Fig. I4.6[ The 
best way to prove the presence of multi-body is by looking for D or D* mesons in 
a kinematically forbidden region. That is, we look for D mesons in a momentum 
region where one would expect none under the assumption of pure two-body 
events. Fig. 15.21 shows a plot of the invariant mass of — > K~Ti^ at the center- 
of-mass energy 4260 MeV for candidates with momenta less than 300 MeV/c. 
Under the assumption that only pure D*D*, D*D, and DD are present, no 
are allowed below ~ 500 MeV/c. The figure shows a clearly defined peak at the 
mass, demonstrating the multi-body contribution. 

This demonstrates that charm is produced through more than the two-body 
event categories, but it sheds no light on the composition. Are the multi-body 
events DDtt, DDtttt, D*D7t, etc., or a combination of all possible types? To 
help answer this question a similar study to the one above, was performed for 
Fig. [53] shows a plot of the AM spectrum for pr)*o < 500 MeV/c. If 
the assumption of only pure two-body events is made, then no D*^ decays are 
kinematically allowed below ~ 700 MeV/c. The figure clearly shows that there 
exist multi-body events of the type D*Dn in the data, evident in the well-defined 
peak located at AM = M^.o - Mdo = 0.142 MeV/c^ 

Looking in the kinematically forbidden region has given clues to the possible 

composition of the multi-body events, but still does not provide a definitive and 

quantitative breakdown. One reason is that initial-state radiation (ISR), can 

102 




Mass (GeV/cO 



Figure 5.2: The measured D'^ — > K~tt~^ invariant mass distribution for pno < 
300 MeV/c. If the assumption of pure two-body events is made, then no 
decays are kinematically allowed below ~ 500 MeV/c at 4260 MeV. 



lead to D mesons smeared outside of the two-body kinematic regions. A more 
definitive test is to reconstruct a D*, add a charged or neutral tt, and look at the 
missing mass (recoil mass) of the event. A peak in the missing mass at the D 
would clearly demonstrate the presence of D*Dt: in the data. We concentrate on 
multi-body events of the type 0^*"^^ D^t:"^ ^ which by isospin are twice as likely to 
occur as D'^*^^D^t^^ and D(*)+iJ~7r°. In addition to the factor of two from isospin, 
the multi-body events with a charged pion should be cleaner than those with a 
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Figure 5.3: The measured D*^ D°7r° AM distribution for D° K-tt+ satisfy- 
ing pd*o < 500 MeV/c. If the assumption of pure two-body events is made, then 
no D*^ decays are kinematically allowed below ~ 700 MeV/c at 4260 MeV. 
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neutral pion, since the additional pion will be soft, < 150 MeV/c. With this 
method the observation of multi-body cannot be obscured by ISR, because the 
presence of the radiative photon will prohibit peaks in the missing-mass spectrum. 

For this study we use the high-statistics data sample collected at 4170 MeV, 
consisting of 179 pb'^ Using D*+ D+7r° with D+ K-7t+7t+, and D*^ 
D°7r° with K-TI+, R-n+n^, or D° K-n+n+n-, in addition to 

DTAG-like cuts for the charged pion, we obtain the missing-mass spectra shown 
in Figs. 15.41 and 15.51 Fig. 15.41 is the invariant mass spectrum of X in e~^e~ — *• 
D*'^7i^X, while Fig. 15.51 is for X in e+e^ D*^n'^X, where for the latter decay 
the charge of the Z)-daughter kaon is used to obtain the correct combination of 
the neutral D* and the charged pion. For both cases we define the signal region as 
a 6 MeV window centered on the previously measured (PDG) mass difference. A 
cut on the reconstructed D* momentum of 400 MeV/c was used to exclude two- 
body events. These figures provide conclusive evidence that multi-body events 
of the form D*D7i exist in the data. We must now determine the amount of 
multi-body present to assess the effect on the two-body cross sections that were 
determined earlier. 

We performed a similar study with the data collected at 4260 MeV. The 

corresponding missing-mass spectrum, only for D*^ D^tt^ with —>■ K^n^, 

K-n+n'^, or R-n+n+n-, is shown in Fig. for the 13.11 pb^^ 

collected at 4260 MeV. A cut of 500 MeV/c was made on the reconstructed D* 
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Figure 5.4: The mass spectrum of X, in e~^e~ — > D n^X , using the 179 pb~ 
collected at 4170 MeV. The peak at the D mass shows conclusive evidence that 
events of the form D*Dt: are present in the data. A cut on the reconstructed D* 
momentum of 400 MeV/c was used to exclude two-body events. 
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Figure 5.5: The mass spectrum of X, in e+e" D*^tt'^X, using the 179 pb"^ 
coUected at 4170 MeV. The peak at the D mass shows conclusive evidence for 
production of D^Dtt events. A cut on the reconstructed D* momentum of 400 
MeV/c was used to exclude two-body events. In addition, the charge of the D- 
daughter kaon is used to obtain the correct combination of the neutral D* and 
the charged pion. 
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momentum to select entries from the multi-body region. In addition, the charge 
of the -D-daughter kaon is used to obtain the correct combination of the neutral 
D* and the charged pion. It is clear from the plot that D*D*tt multi-body events 
exist, in addition to D*Dn, at 4260 MeV. This is shown by the clear peaks at 
both the D and the D* masses in Fig. 15.61 



25 



20 



15 



> 

0) 

m 

00 
(A 

8 10 



5- 



n 



_L 



1.75 1.8 1.85 1.9 1.95 2 
Missing l\/lass (GeV/c^) 



2.05 



2.1 



Figure 5.6: Mass spectrum for X, in e~^e^ D*'^'k^X, at 4260 MeV. The charge 
of the daughter kaon is used to obtain the correct combination of the neutral 
D and the charged pion. The peak at the D mass is additional evidence of the 
existence of D*D7i, whereas the peak at the D* mass is evidence of the production 
of D*D*7T. A cut of 500 MeV/c was made on the reconstructed D* momentum 
to select entries from the multi-body region. 
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Having demonstrated unambiguously that events of the form D*Dt: are present 
in our energy region, it remains to be determined if other types of events, hke 
DDn, also contribute. We performed a similar study to the one above using a D 
rather than a D* to investigate this possibility in the 4170 MeV data sample. Us- 
ing only — s> K~TT~^ and — * K~7i^7i~^, in addition to another charged pion, 
we obtained the missing-mass spectra in Figs. l5.7 land l5.8[ Fig. 15. 71 shows the mass 
of X in e^e^ D^n^X, while Fig. 15.81 gives the mass of X in e~^e^ D^n^X. 
For the latter decay we use the charge of the D-daughter kaon to obtain the cor- 
rect combination of the neutral D and the charged pion. In both plots, the signal 
region is defined to be a 30-MeV window centered on the known D mass. A cut 
on the reconstructed D momentum of 250 MeV/c was made to exclude two-body 
events. In both Figs. l5.7l and l5.8l clear peaks at the D* mass give further evidence 
for D*Dtt multi-body events. The absence of a peak at the D mass indicates that 
there is no evidence for DD-k production, although such a contribution may still 
be present in the data and suppressed by the D-momentum cut. 

Now that the components of the multi-body have been identified, a method 

to measure its yield is needed so that we can correct the exclusive D*D* and 

total exclusive charm cross sections for the portion of the multi-body that was 

excluded. To accomplish this, we use a two-body MC representation of the various 

exclusive channels and a spin-averaged phase-space model MC representation of 

D*Dtt, and D*D*7r, and fit the sideband-subtracted momentum spectrum for 
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Figure 5.7: The mass spectrum of X, in e+e" — > D'^n^X, using the 179 pb~^ 
collected at 4170 MeV. A cut on the reconstructed D momentum of 250 MeV/c 
was made to exclude two-body events. The peak at the D* mass shows, again, 
conclusive evidence for events of the form D*D7i being present in the data. The 
lack of a peak at the D mass shows the lack of DDtt type events in the data. 
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Figure 5.8: The mass spectrum of X, in e+e~ — ^ D'^tt^X, using the 179 pb^^ 
collected at 4170 MeV. A cut on the reconstructed D momentum of 250 MeV/c 
was made to exclude two-body events. In addition, the charge of the £)-daughter 
kaon is used to obtain the correct combination of the neutral D* and the charged 
pion. The peak at the D* mass shows further evidence for events of the form 
D*Dti. The absence of a peak at the D mass indicates the lack of DDti events 
in the data. 
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In the sideband-subtraction method the signal region, in invariant mass, is defined 
to be ±3(7, where a is 5 MeV, about the PDG mass. The low-side and high- 
side sidebands start at 7.5 a away from the peak and extend outward by 3 a. 
The reconstructed D momentum for each region is plotted and the resulting 
histograms are subtracted. The fit made at 4170 MeV is shown in Fig. 15.91 The 
fit seems to replicate the overall structure reasonably well, but is clearly lacking 
in some regions. 

Besides fitting the — >■ K^tt'^ momentum spectrum, a fit to the — »• 
K^Ti^Ti^ can be a good cross-check and test of the method. Fig. 15. 101 shows this 
test using all the assumptions made up to this point. The fit is not as good as 
the and it is clear that that something is missing. The disagreement between 
the fit and the data is dramatic between 600 and 750 MeV/c. 

We considered three possible explanations for this discrepancy: (1) the branch- 
ing fraction for the gamma decay channel is wrong, (2) DDn is present in the 
data; and (3) something more fundamental is missing. The first possibility seems 
unlikely, because enhancements are not seen in the region around D*^D*^ outside 
that of the 7r° enhancement, like those seen in the —>■ K~Tr~^ spectrum. The 
second possibility seems more likely since all the searches for multi-body thus far 
only looked in momentum regions that limited the allowable amount of DDtt. If 

we take Fig. 15.101 as evidence for the presence of DDn, we can modify our fit 
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Figure 5.9: Sideband-subtracted momentum spectrum for K^n^ at 

4170 MeV. Data are shown as the points with error bars which are fit to the 
MC (histograms). The fit uses a spin-averaged phase-space model MC represen- 
tation of D*Dtt, shown in dark red. The fit seems to rephcate the overall structure 
reasonable well, but is clearly deficient in some regions. 
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Figure 5.10: Sideband-subtracted momentum spectrum for — > K~ti^ti^ at 
4170 MeV. Data are shown as points with error bars which are fit to the MC 
(histograms). The fit uses a spin-averaged phase-space model MC representation 
of D*D7^, shown in dark red. There is a large discrepancy between what is 
expected and what is seen around 600-750 MeV/c. This discrepancy suggests 
that DDtt may be present at this energy. 
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procedure to allow a DD-n contribution. This greatly improves the fit quality, 
which can be seen in Fig. 15 .111 for K^-k^ and Fig. 15. 12] for — > K^-k^-k^ . 
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Figure 5.11: Sideband-subtracted momentum spectrum for K^n^ at 

4170 MeV. Data are shown as points with error bars which are fit to the MC 
(histograms). The fit uses a spin-averaged phase-space model MC representation 
for D*D7i, shown in dark red, and DDn is shown in black. 



While the fits are improved, we find large discrepancies in the D^/D^ ratio 
for D*Dti using the composition obtained from the fits. The ratio of a{D^)jj*i)^ 
to a{D^)fjtf)^ should be 2, while the result from the fit is about 5. Also, the ratio 
of a{D^)jjj^^ to a{D^)j^i)^ should be 1, but the result from the fit is about 2. So 
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Figure 5.12: Sideband-subtracted momentum spectrum for — > K^-k^-k^ at 
4170 MeV. Data are shown as points with error bars which are fit to the MC 
(histograms). The fit uses a spin-averaged phase-space model MC representation 
for D*D7i, shown in dark red, and DDn is shown in black. 
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even though the fits look good, the corresponding results are not consistent and 
further investigation is warranted. 

Using all the data collected at 4170 MeV, it is possible to break down the 
ratio of to D+ as a function of energy. The ratio can be seen in Fig. 15.131 
Notice that the ratio is not corrected for efficiency, which is found with MC to be 
constant across the momentum range. Fig. 15. 131 shows the ratio both for the 4170 
MeV data set and for a MC sample that only assumes D*D7i multi-body and 
implements a simple model for the effects of ISR. The ISR model treats all charm 
production as coming from the '?/'(4160) resonance, which has a width of 78 MeV 
[2j. There is a large discrepancy between the MC and data in the momentum 
range 200-400 MeV/c. This discrepancy would be even worse if DDtt were added 
to the MC. While this particular model does solve the problem of the discrepancies 
seen in the momentum fits, it does suggest the need for improvements in our MC. 

To investigate further, we similarly computed the ratio of D*^ to D*^, again 

not correcting for efficiency, which is found to be constant across the momentum 

range. The result is shown in Fig. I5.14[ The same MC sample with a simple 

resonance model is used. The most surprising result can be seen at 650 MeV/c 

in the region between the peaks for the two-body events D*D* and DD*. In 

the limit of perfect resolution, there should be two delta functions, one for D* in 

D*D events and one D* in D*D* events, where the 650 MeV/c bin falls right in 

the middle. This is true even accounting for resolution, and therefore, the ratio 
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should be zero. That this ratio is zero for MC but not data demonstrates that 
our simple ISR MC is incorrect and a better model is needed. 

Our improved approach is to use EVTGEN's model of EvtVPHOtoVISR [57], rep- 
resenting the cross sections for each event type of two-body events (Figs. 14.11 
and 14.21) with simple linear interpolations between our measured data points, as 
is shown in Figs. 15.151 and 15.161 For multi-body constant cross sections are as- 
sumed. Momentum fits to these models give the results for K^-k^ and 
K^Ti^-K^ shown in Fig. 15.171 and Fig. I5.18[ respectively. In addition to 
fitting and one can fit the Dg momentum distribution using the same 
procedure. This result, only for D"^ ^tt^, is given in Fig. 15.191 

Momentum-spectrum fits for data collected at other energy points are shown 
in Figs. [HIinHH32] for and and Figs. [H31HH3S1 for D^. It is interesting to see 
the emergence of the various final states as the center-of-mass energy is increased. 

Since the fit for is independent of , a ratio of these fit results, calculated 

for each event type, can be a good check of the method. The ratios for all event 

types and all energies are shown in Fig. 15.201 At all energies the ratios are in good 

agreement with what is expected (the solid horizontal lines) . As a further check, 

a MC sample with the same statistics as the 178 pb~^ of data was generated. 

Figs. I5.2T| 15.221 and 15.23] show the scaled momentum spectrum for K^n^, 

K^n^7r~^, and 07r+, respectively, for this sample overlaid with the 

data. The agreement is quite exceptional. 
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For another confidence test of tlie procedure we used tliis one-times sample 
to verify the determination of the production ratios. Comparisons of D*^ to D*~^ 
and to for the updated MC and data can be seen in Fig. 15.241 and Fig. 
15.251 Overall, the previously observed discrepancies between the data and MC 
have been resolved with the use of a more accurate description of ISR in the MC. 

The final measurements of multi-body production at energies above 4030 MeV 
are summarized in Table 15.11 The exclusive cross sections for all modes studied 
from the momentum-spectrum fits are shown in Tables [8. 101 and [8. Ill as well as in 
Figs. 15.261 15.271 and 15.281 The total charm cross section for all three methods, as 
discussed in the text, are shown in Fig. 15.291 It should be noted, unless otherwise 
stated, that the weighted sum technique from Sect. 14. II is used for the individual 
Ds cross sections whose result is used in the determination of the total exclusive 
charm cross section, e.g. Fig. I5.29[ 
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Table 5.1: The amount of multi-body present at energies above 4060 MeV ob- 
tained by fitting the sideband-subtracted momentum spectra for D° K~7r~^ 
and D'^ — > K~tt~^7i~^ with the two-body MC representation of the various ex- 
clusive channels and a spin-averaged phase-space model MC representation of 
multi-body. There is no evidence for the DDn final state at any energy. Only 
statistical errors are shown. 



Ecm MeV 


a{e+e- -> D*D7r) nb 


a{e+e- D*D*n) nb 


4060 


0.14 ±0.09 


-±- 


4120 


0.05 ± 0.08 


-±- 


4140 


0.41 ± 0.09 


-±- 


4160 


0.39 ±0.06 


-±- 


4170 


0.44 ±0.01 


-±- 


4180 


0.58 ±0.09 


-±- 


4200 


0.74 ±0.13 


- ± - 


4260 


0.64 ±0.09 


0.32 ±0.07 
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Figure 5.13: Ratio of the number of to as a function of reconstructed 
D momentum at E^ra = 4170 MeV. The differences in the branching fractions 
were included in the ratio, while the efficiencies were take to be constant over 
the momentum range. The MC includes only D*D7t multi-body production and 
uses a simple ISR model that assumes all events are produced from the ^p{4160) 
resonance, with a width of 78 MeV. The large discrepancy between the data and 
the MC between 200 and 400 MeV/c clearly indicates that something is incorrect 
in our assumptions used to generate the MC sample. 



121 



3 



+ 

Q 

T 

+ 

* 

T 



OQ 



^.5 

o 

Q 

T 

Q 

cE 

OQ 

* 

S.5 

T 

Q 

£ 1 



0.5 



2- 



lySpb"' DATA 



▼ ISR MC, only D'Dti MB 



0.1 



0.2 



0.3 



0.4 0.5 

.2x 



0.6 



0.7 



0.8 



Pp. (GeV/cO 



Figure 5.14: Ratio of D*^ to D*^, both using the 7r° channel, with the subsequent 
D-meson decays K^n^ and K^tx^ti^ , respectively. The ratio 

includes the branching fractions but excludes the efficiencies, which are constant 
over the momentum range. The MC includes only D*Dtt multi-body production 
and uses a simple ISR model that assumes all events are produced from the 
'0(4160) resonance, with a width of 78 MeV. The discrepancy at 650 MeV/c 
indicates that ISR is important and that it is incorrectly handled in our MC 
sample. 
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4.1 4.15 



Figure 5.15: The two-body cross sections used in EVTGEN for DD, D*D, and D*D* 
for the MC samples with improved treatment of ISR. A simple implementation 
of the cross sections were made using straight-line approximation for the real, 
Born-level, two-body, cross sections. The main assumption is that the shape of 
the cross section will be preserved after applying radiative corrections, that is the 
observed cross section shape is the same as the Born-level cross section shape. 
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Figure 5.16: The two-body cross sections used in EVTGEN for DgDs, D^Dg, and 
D*D* for the MC samples with improved treatment of ISR. A simple implemen- 
tation of the cross sections were made using straight-line approximation for the 
real, Born-level, two-body cross sections. The main assumption is that the shape 
of the cross section will be preserved after applying radiative corrections, that is 
the observed cross section shape is the same as the Born-level cross section shape. 
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Figure 5.17: Sideband-subtracted momentum spectrum for —>■ K^tt^ at 
4170 MeV. Data are shown as points with error bars which are fit to the im- 
proved ISR MC (histograms). The fit uses a spin-averaged phase-space model 
MC representation of D*Dn, shown in dark red. The sohd histograms corre- 
spond to < 1 MeV (no ISR) while the dashed histograms correspond to 
> 1 MeV (with ISR). The total fit result is shown by the solid black line 
ix^/NDF = 248/132). 
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Figure 5.18: Sideband-subtracted momentum spectrum for K^tt^tt^ at 

4170 MeV. Data are shown as points with error bars which are fit to the improved 
ISR MC, the histograms. The fit uses a spin-averaged phase-space model MC 
representation of D*D7r, shown in dark red. The sohd histograms correspond 
to < 1 MeV (no ISR) while the dashed histograms correspond to -E^ > 1 
MeV (with ISR). The total fit result is shown by the solid black line {x^/NDF = 
182/132). 
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Figure 5.19: Sideband-subtracted momentum spectrum for —>■ (p-n^ at 

4170 MeV. Data are shown as points with error bars which are fit to the im- 
proved ISR MC (histograms). The sohd histograms correspond to -E^ < 1 MeV 
(no ISR) while the dashed histograms correspond to E"^ > 1 MeV (with ISR). 
The total fit result is shown by the solid black line {x^/NDF = 165/124). 
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Figure 5.20: The ratio of the cross section results from the K^n^ fit to the 

results from the —>■ K^-k^-k^ fit. Since the fits are independent, the agreement 
in the ratios demonstrates the reliability of the fitting procedure in extracting the 
mode-by-mode cross sections. The ratio, in all cases, should be one, indicated by 
the solid black line. All event types are in good agreement with what is expected. 
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Figure 5.21: Comparison between MC and data at 4170 MeV for K^tx^ . 

Based on the momentum spectra fit results, an indepedent MC sample with the 
same statistics as the 178 pb^^ collected at 4170 MeV was generated using all 
the knowledge gained about ISR, angular distributions (see the App.) and cross 
sections. 
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Figure 5.22: Comparison between MC and data at 4170 MeV for — > K^tv^tv^. 
Based on the momentum spectrum fit results an independent MC sample with 
the same statistics as the 178 pb~l collected at 4170 MeV was generated using all 
the knowledge gained about ISR, angular distributions (see the App.) and cross 
sections. 
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Figure 5.23: Comparison between MC and data, at 4170 MeV, for — » 07r+. 
Based on the momentum spectrum fit results an independent MC sample with 
the same statistics as the 178 pb~l collected at 4170 MeV was generated using all 
the knowledge gained about ISR, angular distributions (see the App.) and cross 
sections. 
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Figure 5.24: Ratio of to in 4170 MeV data and MC as a function of the 
reconstructed D* momentum. Both D* states are detected with the vr" channel, 
with the subsequent decays —>■ K^n^ and K^ti^ti^, respectively. The 

ratio is corrected for branching fractions, while efficiencies were assumed to be 
constant. The MC includes D*Dtt multi-body and uses an updated ISR model 
that incorporates individual two-body cross sections as described in the text. The 
updated MC sample is in great agreement with the data at 4170 MeV. 
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Figure 5.25: Ratio of the number of D to as a function of reconstructed D 
momentum at = 4170 MeV. The ratio is corrected for branching fractions, 
while efficiencies were assumed to be constant. The MC includes D*Dn multi- 
body and uses an updated ISR model that incorporates individual two-body cross 
sections as described in the text. The updated MC sample is in great agreement 
with the data at 4170 MeV. 
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Figure 5.26: Observed cross sections for e+e^ — > DD, D*D and D*D* as a 
function of center-of-mass energy. They are determined by fitting the momentum 
spectrum for D^ —> K~7i~^ and D~^ —>■ K~n^7i^ with the updated ISR MC. 
Systematic errors are not included. 
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Figure 5.27: Observed cross sections for e+e^ DsDg, D*Ds and D*D* as a 
function of center-of-mass energy. They are determined by fitting the momentum 
spectrum for 07r+ with the updated ISR model MC. A branching fraction 

of 3.6% was also used in the determination of the cross section. Systematic errors 
are not included. 
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Figure 5.28: Observed cross sections for e~^e 0*0-71 and D*D*tt as a function 
of center-of-mass energy. They are determined by fitting the momentum spectrum 
for K'TT+ and D+ K-TT+n+ with the updated ISR model MC. The 

multi-body contribution was modeled in MC using spin-averaged phase space. 
Systematic errors are not included. 
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Figure 5.29: Plot of the total charm cross section as calculated by each of the three 
methods described in the text, in addition to the results from the momentum- 
spectrum fit. The reason for the discrepancy between the inclusive and exclusive 
method is due to the presence of multi-body background and ISR. These two 
effects are taken into account when a fit to the momentum spectrum is performed. 
There is good agreement between all three methods, the inclusive, the hadron 
counting and the momentum fits. The previous exclusive result (ignoring multi- 
body) is shown, in addition, for comparison. 
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5.2 Momentum Spectrum and Multi-Body Cross- 
check 

Using the improved 178 pb^^ MC sample that was generated to check the neutral- 
to-charged D-meson ratio, one can check the momentum-spectrum fit procedure. 
The momentum-spectrum fit results for K^n^, K^n^n^, and 

— s> 07r+ are shown in Figs. 15.301 15.311 and 15.321 respectively. The comparison 
between what was generated and what is obtained from the fit is shown in Table 
15.21 The returned fit results are in very good agreement with what was used to 
generate the sample. 



Table 5.2: Results of the momentum-spectrum fit of the improved 178 pb~^ MC 
sample. The agreement between the fitted and the generated values is very rea- 
sonable. 



Parameter 


Fitted Value (nb) 


Generated Value (nb) 


Difference 


DD 


0.319 ±0.008 


0.331 


-1.5(T 


D*D 


2.392 ±0.08 


2.392 


-O.lcr 


D*D* 


4.553 ±0.038 


4.508 


1.2(7 


D*Dn 


0.766 ±0.018 


0.780 


-0.8(T 


DsDs 


0.036 ±0.005 


0.040 


-0.8cT 


DID, 


1.039 ±0.028 


0.994 


1.6a 



As an additional check on the multi-body contributions at 4170 MeV and 
4260 MeV, the two energies where either a large amount of data existed or a 
large amount of multi-body is present, one can fit the missing-mass spectrum to 
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Figure 5.30: Momentum-spectrum fit for tlie improved 178 pb^^ MC sample for 
—>■ K^TT^. The MC sample was treated just as the data and used as a cross- 
check of the momentum-spectrum fit procedure. 



139 



2500 



2000 
> 

0) 

m 

SJ500 

(A 

c 

3 
O 

^000 



500 



▼ 1xMC Sample 
— Total fit 




▼ DATA 

Total Fit Result 

DD 

D*D primary 

D*D 7i:° channel 

D*D Y channel 

D*D* 71° channel 
D*D* Y channel 
D*D7i 



0.4 0.6 

P (GeV/c) 



Figure 5.31: Momentum-spectrum fit for tlie improved 178 pb^^ MC sample for 
K^n^Ti^. The MC sample was treated just as the data and used as a 
cross-check of the momentum-spectrum fit procedure. 
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Figure 5.32: Momentum-spectrum fit for tlie improved 178 pb^^ MC sample for 
07r+. The MC sample was treated just as the data and used as a cross- 
check of the momentum-spectrum fit procedure. 
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obtain the cross section for multi-body. The fit results for D*^ D^n^ with 
— > K^n'^, and D*^ — > D+tt*^ with — > K^n'^n^ plus an additional charged 
pion at 4170 MeV are shown in Figs. 15.331 and 15.341 respectively. The fit results 
for D^n^ with D° K-n+, D° K-n+n\ or R-n+n+n- plus an 

additional charged pion at 4260 MeV is shown in Fig. 15.351 Using a spin-averaged 
phase-space model MC sample to determine the corresponding efficiencies the 
cross sections can be determined. The results from this method are shown, along 
with the momentum-spectrum fit results, in Table 15.31 The agreement between 
these two different methods supplies more confidence in the results that are being 
presented. 



Table 5.3: Comparison between the momentum-spectrum results and fits to the 
missing-mass spectrum for multi-body at 4170 MeV and 4260 MeV. The agree- 
ment between the two different methods instills confidence on our handling of 
multi-body events. 



Parameter 


Energy (MeV) 


Missing-Mass 
Value (nb) 


Momentum-Spectrum 
Value (nb) 


D*Dtt{D*^) 


4170 


0.34 ±0.06 


0.44 ±0.01 


D*D7i{D*+) 


4170 


0.42 ±0.06 


0.44 ±0.01 


D*D7i{D*^>) 


4260 


0.62 ±0.16 


0.64 ±0.09 


D*D*7i{D*^) 


4260 


0.23 ±0.05 


0.32 ±0.07 
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Figure 5.33: Missing-mass spectrum using D° ^ K n'^ at 4170 MeV. The fit 
function consists of a Gaussian and an Argus function. 
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Figure 5.34: Missing-mass spectrum using — > K n'^n'^ at 4170 MeV. The fit 
function consists of a Gaussian and an Argus function. 
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Figure 5.35: Missing-mass spectrum using — > K~7r^, — * K'TT'^n^, and 
K^n^n^n^ at 4260 MeV. The fit function consists of two Gaussians, one 
for each peak, and an Argus function. 
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Chapter 6 



SYSTEMATIC S 

6.1 Exclusive Cross Sections via Momentum 
Spectrum Fits 

Using the published CLEO-c 56 pb^^ analysis [l5l US], at all energies we apply 
a 0.7%, 0.3%, and 1.3% systematic for tracking, pion PID, and kaon PID re- 
spectively. In addition, we apply a 1% systematic to the luminosity [58]. Also, 
the error on the respective branching ratios are 3.1% for K^7f~^, 3.9% for 

[IHl SS], 0.7% for D*+ L)°7r+ [2] and the respective errors 
associated with Ds decays pi HI] . 

Besides the aforementioned systematic errors, the additional area where sys- 
tematic errors will be present is in the cross section shape assumed in EVTGEN [57j . 
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The systematic error resulting from the shape of the cross section was determined 
by adjusting the two-body exclusive cross sections from their nominal shapes (Fig. 
15. isl and Fig. 15.161) . The adjustments were determined by the existing data points 
as well as taking some extremes to understand the effect of the shape of the cross 
section assumed in EVTGEN. The systematic errors were investigated on the large 
data-set at 4170 MeV since statistical errors are minimized. Lastly, since an ad- 
justment to the shape of single cross section can effect the other cross sections in 
the fit, the total systematic error for an individual cross section is determined by 
adding in quadrature all cross section effects. 

6.1.1 DD 

Two adjustments were investigated in regards to the shape of the cross section. 
First, a step was added starting at 4015 MeV and extending to 4120 MeV where 
the cross section was doubled. Since the cross section in this range is larger, as 
compared to the nominal assumption of a flat cross section across all energies, the 
probability distribution will reflect this change. Events which radiate a photon 
for this energy range are now more probable and the fit results for this assumed 
cross section will be larger as compared to the nominal assumption. The next 
adjustment investigated is a decreasing cross section as a function of increasing 
energy, the slope is 2 nb/GeV. Following similar logic as before, if the cross section 
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at lower energies is larger, as compared to the cross section at 4170 MeV and the 
nominal flat cross section then the fit result will report a larger result for the DD 
and DgDs cross sections as is shown in Table [6?T1 Similar results are obtained for 
a cross section that is smaller at lower energies, as compared to the cross section 
at 4170 MeV. Since increasing and decreasing the cross sections, at 2 nb/GeV, 
are extreme adjustments, whereas the box is a conservative adjustment, the two 
are averaged and used in the determination of the systematic errors. The former 
is considered extreme because the data does not suggest this type of behavior. 



Table 6.1: Relative systematic errors for adjustments to the DD cross section. 



Event Type 


Box 4015-4120 MeV 


Slope (m = 2 nb/GeV) 


Average 


DD 


+2.1% 


+4.2% 


+3.2% 


D*D 


-0.5% 


-0.5% 


-0.5% 


D*D* 


+0.1% 


+0.1% 


+0.1% 


D*D7i 


+0.3% 


+0.2% 


+0.3% 



6.1.2 DW 

Two adjustments were investigated in regards to the shape of the cross section. 
First, the cross section was adjusted to take out the kink at 4060 MeV allowing 
the cross section to decrease constantly and continuously from 4015-4170 MeV. 
Since the cross section is increasing one would expect the fit result to reflect this 
change and increase. However, since there is an overlap of both DD and D*D* , 
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the simple explanation is no longer straight forward. The next adjustment was to 
assume a flat cross section until 4060 MeV at which energy the cross section then 
increased to 4015 MeV. This assumption results in a poor fit to the momentum 
spectrum, most notable in the K^ir^n^ distribution as seen in Fig. 16.11 

Therefore, because the latter assumption produces a poor fit, systematic errors of 
the former are used. Since there is overlap between the various two-body states 
present at this energy changes to the D*D effects to the other cross sections result 
which are evident in Table 16. 2[ 



Table 6.2: Relative systematic errors for D*D. 



Event Type 


Continuous Decrease From 4015 MeV 


Flat to 4060 MeV 


DD 


+0.6% 


+0.8% 


D*D 


-1.1% 


-5.2% 


D*D* 


+0.4% 


+3.0% 


D*Dn 


+3.0% 


+2.0% 



6.1.3 DW* 

For these events we investigated changes to the shape of the cross section in 
addition to changes in the respective hehcity amplitudes. First, the cross section 
for D*D* was adjusted such as to add a dip in the plateau region of the nominal 
cross section. In a similar situation as D*D, the addition of other possible two- 
body states along with the addition of multi-body a simple expectation is difficult 
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Figure 6.1: The fit result after adjusting the D*D cross section from its nominal 
value to one that is flat between 4060-4170 MeV. The change results in a poor fit 
most notable at 600 MeV/c. 
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to ascertain. Another adjustment was to remove the kink at 4030 MeV and to 
allow the cross section to decrease continuously to zero at 4060 MeV. The effects 
of both adjustments on the D*D* along with the other cross sections is shown in 
Table 

Table 6.3: Relative systematic errors for adjustments in the shape of the cross 
section for D*D*. 



Event Type 


Dip 


Decrease to Zero at 4070 MeV 


DD 






D*D 




+0.2% 


D*D* 


+2.4% 


-2.4% 


D*Dn 


-1.5% 


+9.1% 



The helicity amplitude affect the angular distributions of the D*D* final state 
and therefore has an effect of the resulting daughter D (See Appendix for a dis- 
cussion of angular distributions). This shows up, most notably, in the momentum 
spectrum resulting from the vr decay channel, between 500-600 MeV/c. Taking an 
extreme case we adjusted the coefficient in front of the cos for the angle between 
the D in the rest frame of the D* with respect to the momentum of D* in the lab 
from its nominal value of 0.8 to 2.0. The result is shown in Fig. 16.21 which clearly 
shows that such a change results in a bad fit. The effect on the cross sections is 
shown in Table 16.41 
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Figure 6.2: The fit result after adjusting tlie angle between the D in the rest 
frame of the D* with respect to the momentum of D* in the lab from its nominal 
value of 0.8 to 2.0. The fit is poor most notably between 500-600 MeV/c. 
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Table 6.4: Relative systematic errors for adjustments in the helicity amplitudes 
for D*D\ 



Event Type 


Angle {a = 2) 


DD 




D*D 


+0.3% 


D*D* 


-1.0% 


D*Dn 


+6.0% 



6.1.4 D*Dtt 

Rather than assuming a flat cross section as a function of energy, an adjustment 
that allows the cross section to increase as the energy increases. Since the cross 
section is now smaller at lower energies, as compared to the cross section at 4170 
MeV, the fit result will return a larger result. Because the momentum range of 
the multi-body event is large, between 0-500 MeV/c, such a change in the cross 
section will result in modest changes to the fit results which are shown in Table 

ESI 



Table 6.5: Relative systematic errors for adjustments made to the shape of the 
cross section for D*D7r. 



Event Type 


Increasing Slope (m = 4 nb/GeV) 


DD 




D*D 




D*D* 




D*Dn 


+2.7% 
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6.1.5 DWiT 

Since there is only one data point for 0*3*71 a systematic of 25% is assigned. 

6.1.6 Total Systematic Errors: Exclusive Momentum Fit 
Method 

The total systematic errors for the exclusive cross sections by fitting the momen- 
tum spectrum are shown in Table 16.61 

Table 6.6: Total relative systematic errors for the exclusive momentum fit method. 



Event Type 


Total 


DD 


4.5% 


D*D 


3.4% 


D*D* 


4.7% 


D*Dn 


12% 


D*D*7i 


25% 



6.1.7 Systematic Errors: L>,L>„ D^Ds and D^Dl 

For DgDs, D^Ds and D^D*, a weighted sum technique is used to combine the 
eight Dg decay modes and obtain the cross sections. This technique was used to 
minimize the error, both statistical and systematic. Since multi-body type events 
involving a D^-meson are not expected by isospin which was confirmed by not 
observing evidence of multi-body events in the momentum spectra of Dg ^tt"*" 
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at any of energies, the technique described in Sect. I4.1l is used. Besides the errors 
in regards to the branching fractions, shown in Table 16. 8[ are errors resulting 
from the technique. It should be mentioned that because of the recently updated 
Dg branching fractions [44] the weights used in the calculation are not optimal 
because they were determined using the PDG 2005 [2] values. However, the 
changes in the branching fractions have a minor, if no, effect on the final result 
by virtue of the technique. To demonstrate this we look at the results for Dg at 
4170 MeV. By comparing both the individual mode results and the weighted sum 
result. Table 16. 7^ one sees that the agreement between them is quite good and 
the updated branching fractions reduce the associated systematic error without 
effecting the quoted value for the Dg cross sections stated in Sect. 14. 1[ 



Table 6.7: A comparison between the D*Ds cross sections at 4170 MeV deter- 
mined by individual modes using the recently updated branching fractions [H] 
to that of the weighted sum technique using branching fractions based on [2]. 
The agreement between the individual modes and the weighted sum stresses the 
fact that the result from the weighted sum is not effected by the changes in the 
branching fractions. The only effect of the updated branching fractions is to 
decrease the systematic errors. 



Mode 


Cross Section (nb) 




0.85 ±0.05 




0.99 ±0.05 




0.90 ±0.04 




0.96 ±0.03 


Weighted Sum 


0.92 ±0.01 



The selection criteria, Mbc and AE, for the various exclusive cross sections 
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Table 6.8: The Dg branching fractions including the updated branching fractions, 
in 10-2. 



Modes 


Branching Fraction 


07r+, 10 MeV cut on the Invariant K+K' Mass [44J 


1.98 ±0.15 


K*<^K+,K*'^ K-TT- ig 


2.2 ±0.6 


7]7r+,r] ^ 77 [21111] 


0.58 ±0.07 


rip^, T] — i> 77, p"*" 7r+7r° [2j 


4.3 ± 1.2 


rj TT"*", ?7 ^ 7r+7r-?7, ^ 77 [21 IHj 


0.7 ±0.1 


7] p^,v[ — s> 'K~^7T~ri,ri 77, — * tt+tt" [2] 


1.8 ±0.5 


^ p+ ^ 7r+7r° [2] 


3.4 ±1.2 


KsK+,Ks TT+TT- [21111 


1.0 ±0.07 



result in added systematic errors. Having adjusted the AE cut window by 10 MeV 
the cross sections for D^Ds changed by a maximum of 3% and therefore 3% is 
applied to these cross sections at all energies. Next, by adjusting the Mbc cut 
window by 10 MeV for D*^D^ the cross sections changed by 2.5% and therefore 
2.5% is applied to these cross sections at all energies. The charge in Mbc was 
asymmetric because the photon decay mode is included and allows for to be 
spread out over a wide range of momenta. Lastly, by adjusting Mbc by 30 MeV for 
D*+jj*- ^]2e result changed by 5% therefore 5% is apphed to these cross sections 
at all energies. The table of systematics for the exclusive method is shown in 
Table EB 
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Table 6.9: Total relative systematic errors for the Ds cross sections. 



Event Type 


Total 


DsDs 


5.6% 


DID, 


5.3% 


d:d: 


6.8% 



6.2 Inclusive D Method 

Using the published CLEO-c 56 pb^^ analysis [IHl at all energies we apply 
a 0.7%, 0.3%, and 1.3% systematic for tracking, pion PID, and kaon PID re- 
spectively. In addition, we apply a 1% systematic to the luminosity ^58j. The 
largest systematic in regards to this method is assigned to the fit function and 
the yields obtained by a gaussian plus a polynomial of degree two. Changing the 
background function to a linear function causes a 2% shift, which is the largest, 
IS seen m the yields of D^ and D+. At energies above the DgD* threshold a 
7% shift, which is largest, is seen in the D, yields and 2% at all energies below 
threshold. In all cases, except 4170 MeV, the change in yields is < Icr. Lastly, 
the error on the respective branching ratios are 3.1% for D^ K^tt^, 3.9% for 
D+ K-n+n+ gSllin] and 6.4% for D+ K'K+n+ [B]. Therefore, we quote 
a systematic of 4.3% for D^ and 5.1% for D~^ at all energies. For Dg an extra 
systematic error of 4%, associated with the selection criteria, is applied. The de- 
termination of the latter systematic error was taken from the investigation of the 

exclusive cross sections for D*^D^ and DfDj. Therefore, Dg is quoted with a 
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systematic error of 8.6% for energies below 4100 MeV and 10.8% for those above. 
The systematic error associated with the branching fractions, as well as the PID 
and tracking, should be improved in the near future with the updated D and Dg 
branching fractions. 

6.3 Hadron-Counting Method 

Since the method that was used here is identical to method used to determine 
the cross section of e+e" ^(3770) hadrons at E^^ = 3770 MeV [iQllSOllST], 
most of the systematics will be identical. One difference that was investigated 
was the systematic that is associated with the hadronic event selection criteria, 
since this could vary with energy. By applying tighter cuts, THAD \5T\ , the 
maximum change in the cross at any energy was 4.5% and so is taken as the 
hadronic event criteria systematic for all energies. In a similar fashion, as in the 
reference [HI EHl EI] , a conservative error, to account for the interference with the 
continuum and because they are based on theoretical calculation, of 25% will be 
associated with ip{2S), J/ip, and '?/^(3770) subtraction at each energy. Since the 
cross section varies for these backgrounds, so does the associated error. Errors 
which are common between the energies are shown in Table 16.101 whereas the 
energy dependent errors are shown in Table 16.111 including the total. 
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Table 6.10: Summary of various systematic errors that are common amongst the 
scan energies for the hadron counting method. 



Source of error 


relative errors in 10 ^ 


QED /Two-Photons subtractions / suppressions 


0.3 


BeamGas/ Wall/Cosmic subtraction 


0.5 


Track quaUty cuts 


0.5 


Luminosity 


1.0 


Continuum scahng 


2.1 


Hadronic event selection criteria 


4.5 


Total common systematic error 


5.1 



Table 6.11: Summary of the energy dependent systematic errors, in 10 ^, for the 
hadron counting method. 



Energy (MeV) 




^{2S) 


^(3770) 


Total relative Systematic Errors 


3970 


3.1 


1.30 


0.44 


6.1 


3990 


2.3 


1.13 


0.34 


5.7 


4010 


1.8 


0.88 


0.25 


5.5 


4015 


1.6 


0.80 


0.23 


5.4 


4030 


1.1 


0.54 


0.14 


5.2 


4060 


1.1 


0.61 


0.15 


5.3 


4120 


0.99 


0.53 


0.12 


5.2 


4140 


0.91 


0.51 


0.12 


5.2 


4160 


0.89 


0.48 


0.10 


5.2 


4170 


0.88 


0.49 


0.10 


5.2 


4180 


0.88 


0.49 


0.11 


5.2 


4200 


0.91 


0.50 


0.12 


5.2 


4260 


1.53 


0.80 


0.19 


5.4 
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Chapter 7 

Radiative Corrections 



The Born or tree-level cross section is obtained by correcting the observed cross 
sections for the effects of initial-state radiation (ISR). In high-energy electron- 
positron annihilation experiments, the incoming particles can radiate. The ra- 
diated photons can be quite energetic, thereby changing the effective center-of- 
mass collision energy appreciably. Therefore, the annihilation energy is not always 
twice the beam energy and the observed cross sections obtained in the experiment 
correspond not to a single energy point but instead to a range of energies. 

The needed correction factors were calculated using two alternate procedures. 
First, by following the method laid out by E.A. Kuraev and V.S. Fadin [56] which 
states that the observed cross section aobs{s) at any energy -y/i can be written as 



1 




(7.1) 
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where cr^ is the Bom cross section as a function of the effective center-of-mass 
energy squared and Sg// = s(l — x), with = xi^^beam- The function, f{x,s) is 
defined as follows: 

fix, s) = + ) + -t In — + 277^ t(l - - 

■'^ ' ' ^ TT^ 3 2^ 4 24^3 ml 4'^ ^ 2' 

+ -[4(2 - x) In - - — ^ ln(l -x)-6 + x], (7.2) 
8 XX 

where t = ^[In(^) — 1]. Eq. 17.21 is only used to first order in t. 

The other method, by G. Bonneau and F. Martin [52], states that the observed 
cross section can be written in terms of the Born cross section as follows: 

2a 2E 13 
o-obs = cTfi 1 H 1(2 In l){\nxmin + 77; + 

^{l-x + -- + -}]. (7.3) 

x^ 2' (t(s) ' 36 6 ^^ ^ ' 

Since the integral is infrared divergent, the Inxmm in Eq- 17.31 corresponds to the 
low-energy cutoff. 

The only difference between the two methods is that G. Bonneau and F. 
Martin include the vacuum polarization for the electron, i.e. the electron bubble: 



2q; r 5 1 , s , ,„ 



e 



We confirmed this by subtracting Eq. 17.41 from Eq. 17.31 and applying to the simple 
test of a constant cross section as a function of center-of-mass energy. 

Note that Eqs. 17.21 and 17.31 are dependent on the tree-level cross section, not 
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only at s, but at all energies below s. Therefore, the ISR correction relies on a the- 
oretical model or on already-radiatively-corrected experimental data to describe 
the shape of the cross section at all relevant energies. 

It was decided to approach determining the shape of the cross section at lower 
energies with two different methods and to test on two different sets of already- 
radiatively-corrected data. The two data sets are from R measurements in the 
region above cc threshold made by BES [IS] and Crystal Ball (CB) ^U\. The 
two shape methods were a simple linear-interpolation procedure applied to each 
data set, and a fit consisting of a sum of Breit-Wigners. Each of these methods 
was applied to the data sets after the uds background had been subtracted using 



The results after radiatively correcting the inclusive cross section are shown 
in Fig. 17. 1[ The difference in the corrected cross section between the two shape- 
determining methods is mainly due to the interpolation method being influenced 
by the "jitteriness" of the data. Since this method is a simple linear interpolation, 
the fluctuations in the data play an important role, especially in the region nearest 
to the point that is to be corrected. Even so, the agreement is quite good, as 
seen in Table 17.11 The agreement between the two radiative correction methods, 
BM and KF, is excellent, giving us confidence that both have been implemented 
correctly. The difference between the methods is taken to be the systematic error 

for the ISR correction, which gives a 4% systematic error on the correction factor 
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for all energies. Since all the methods are in good agreement, the KF interpolation 
method on the CB data has been chosen for the nominal results. The corrected 
inclusive cross section for the energies investigated are shown in Table 17.21 and 
graphically in Fig. 17.21 

A measurement of -Rcharm, sometimes referred to as Rd, can be made by 
dividing acharm by the QED tree-level cross section for muon pair production, 
Eq. II. 2[ In addition to -Rcharm, one needs to determine Ruds, where u,d, and s 
refer to the contribution of the light quarks to the measurement of R. By fitting 
e+e^ hadrons, Fig. 11.41 between 3.2 and 3.72 GeV, with a ^ function we arrive 
at Ruds = 2.285 ± 0.03. The value of R is then computed as the sum of these 
contributions: 

R — Ruds "I" -Rcharm C''-^) 

the results for R are shown in Table 17.31 and graphically in Fig. 17.31 
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Table 7.1: The radiative correction factors used in obtaining the tree- level total 
charm cross section results. The two different shape methods, calculation pro- 
cedures, and the two data sets used are shown. Note that the corrected cross 
sections include the correction due to the vacuum polarization, as defined in Eq. 
17.41 as well as the bubbles due to /i and r leptons. 



Ecm (MeV) 


BESkf 
Inter. 


BESbm 
Inter. 


BESkf 
Fit 


CBkf 
Inter 


CBbm 
Inter 


CBkf 
Fit 


3970 


0.84 


0.82 


0.86 


0.83 


0.81 


0.85 


3990 


0.92 


0.94 


0.85 


0.86 


0.86 


0.83 


4010 


0.85 


0.84 


0.83 


0.82 


0.80 


0.82 


4015 


0.82 


0.79 


0.83 


0.85 


0.84 


0.82 


4030 


0.84 


0.82 


0.83 


0.83 


0.82 


0.83 


4060 


0.86 


0.86 


0.86 


0.92 


0.93 


0.89 


4120 


0.93 


0.94 


0.95 


0.93 


0.93 


0.92 


4140 


0.98 


1.00 


0.93 


0.95 


0.96 


0.91 


4160 


0.95 


0.96 


0.91 


0.96 


0.97 


0.93 


4170 


0.94 


0.96 


0.92 


0.95 


0.96 


0.92 


4180 


0.95 


0.96 


0.94 


0.93 


0.94 


0.97 


4200 


1.00 


1.02 


1.02 


1.02 


1.05 


1.02 


4260 


1.17 


1.22 


1.16 


1.13 


1.17 


1.13 
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Figure 7.1: The radiatively corrected inclusive cross sections using all methods 
are compared to the BES and CB uds-subtracted R data. Note that the corrected 
cross sections (including BES and CB) include the effect of vacuum polarization 
as defined in Eq. 17.41 as well as the bubbles due to /i and r leptons. 
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Figure 7.2: The radiatively corrected inclusive cross sections using the KF 
terpolation method on the CB data. The error bars include the systematics 
discussed in the text. 
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Figure 7.3: R measurement obtained from these inclusive charm cross sections, 
radiatively corrected, using the KF interpolation method on the CB data. The 
error bars include the systematics as discussed in the text. 
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Table 7.2: The total charm cross section as obtained by the inclusive method after 
being radiatively corrected. The errors are statistical and systematic respectively. 



MeV 


a{e+e~ DDX) nb 


3970 


5.92 ±0.22 ±0.28 


3990 


6.89 ±0.25 ±0.33 


4010 


8.60 ±0.22 ±0.41 


4015 


8.97 ±0.42 ±0.43 


4030 


13.10 ±0.35 ±0.63 


4060 


10.85 ± 0.28 ± 0.55 


4120 


9.82 ±0.30 ±0.49 


4140 


9.58 ±0.22 ±0.48 


4160 


9.48 ±0.15 ±0.48 


4170 


9.56 ±0.07 ±0.48 


4180 


9.36 ±0.21 ±0.47 


4200 


7.30 ±0.25 ±0.38 


4260 


3.73 ±0.08 ±0.20 
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Table 7.3: R, after being radiatively corrected, as determined by the inclusive 
method. The errors are statistical and systematic respectively. In this method 
we use Ruds — 2.285 ± 0.03 as determined by a - fit to previous R measurements 
between 3.2 and 3.72 GeV. 



E,,„ Mc>V 


rrfc+c — ■'hddron s) 
a(c ■ (.: ^/t 1 /< ) 


3970 


3.36 ±0.04 ±0.05 


3990 


3.55 ±0.05 ±0.06 


4010 


3.88 ±0.04 ±0.08 


4015 


3.95 ±0.08 ±0.08 


4030 


4.74 ±0.07 ±0.12 


4060 


4.34 ±0.05 ±0.10 


4120 


4.21 ±0.06 ±0.10 


4140 


4.18 ±0.04 ±0.10 


4160 


4.18 ±0.03 ±0.10 


4170 


4.20 ±0.01 ±0.10 


4180 


4.17 ±0.04 ±0.10 


4200 


3.77 ±0.05 ±0.08 


4260 


3.06 ±0.02 ±0.04 
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Chapter 8 



Interpretations and Conclusions 

8.1 Comparisons 

8.1.1 Comparisons to Previous Measurements 

A comparison of the measurements reported in this thesis of a ■ B{D^ Kit) 
and cr ■ B{D^ — > Ktttt) with previous measurements is shown in Fig. 18.11 The 
results from the CLEO-c scan agree quite nicely with previous experiments, and 
are much more precise. 

In addition to the observed branching ratio times production cross section for 
and D^, both Mark III and BES have made similar measurements involving 
the Dg. These results, BES [19] at center-of-mass energy of 4030 MeV and Mark 
III 120] at center-of-mass of 4140 MeV, are shown along with the results of this 
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Figure 8.1: The production cross section times branching ratio for —>■ K^t[^ 
and — ^ K't^^t:^ as a function of energy for this analysis, as compared to 
previous measurements. The results from the CLEO-c scan agree quite nicely 
with previous experiments. 
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Table 8.1: Measurements of observed branching ratio times production cross sec- 
tions. The first uncertainty in each case is statistical and the second (where given) 
is systematic. 





T^iT\ prcfA/ ( A/Tp\/ 1 

1 Ji i\' I y I iv±c V / 


l^^iVnpriTTipnf 








11.2 ±2.0 ±2.5 pb 


4030 


BES 


6.81 ±2.48 pb 


4030 


This Analysis 


■ B{Dt - 07r+) 






26 ± 6 ± 5 pb 


4140 


MARK III 


36.1 ±4.42 pb 


4140 


This Analysis 


a^..^, ■ B{Dt - K<^K+) 






24 ± 6 ± 5 pb 


4140 


MARK III 


23.1 ±4.9 pb 


4140 


This Analysis 


■ B{Dt - K*^'K+) 






22 ± 6 ± 6 pb 


4140 


MARK III 


21.7 ±4.0 pb 


4140 


This Analysis 



analysis in Table 18. 1[ 

During the final preparation of this dissertation the Belle Collaboration [61] 
presented measurements of the cross sections for D^D*^ and D*^D*^ in the 
center-of-mass energy range from 3.7 to 5.0 GeV using ISR events produced in 
e~^e~ annihilation at the T(4S'). Their results are only for the charged D pairs, so 
their results must be multiplied by a factor of 2 for comparison with the measure- 
ments presented in this thesis. Also, they present, up to higher-order radiative 
effects. Born-level cross sections. Based on a previous BaBar analysis Belle 
estimated these higher-order corrections to be on the order of a percent. These 
small corrections are well within the systematic errors, which are ~10%, and are 
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therefore ignored. To do a comparison, I radiatively corrected the exclusive D*D 
and D*D* from Tables ISTTUl and 18.111 following the method described in Sect. [71 
where the shape of the cross section needed for the integral term was taken from 
Fig. 15.151 As can be seen by comparing Belle's result (Fig. 18. 2p to our result (Fig. 
18. 3p the agreement is quite good, keeping in mind the factor of two difference. It 
should be noted that their results cover a wider range of Ecm, but are less precise 
than those present here. 

In addition, the BaBar experiment [63] has presented data on the DD final 
state, also using the ISR technique in e+e~ annihilation at the T{4S). BaBar's 
results are complicated by the fact that they combine several modes which have 
not been efficiency-corrected. Therefore, a direct comparison to their results is 
of limited value, because, the BaBar efficiency is acknowledged not to be flat as 
a function of the invariant mass. However, they do see a dip in the number of 
observed DD events in the neighborhood of 4015 MeV, which is in agreement 
with our observation. 

Taking the partial widths calculated at 4160 MeV by T. Barnes [28], which are 
reproduced in Table [T751 and translating them into percentages, a comparison can 
be made to our results. This is shown in Table 18.21 The disagreement between 
the measurements and Barnes's predictions are quite large, especially for DD, 
D*D and D^D,. 

The analysis results for the cross sections two-body states (Figs. 18.41 18.51) are 
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3.8 4 4.2 4.4 4.6 4.8 5 

Figure 8.2: The exclusive cross sections for e^e^ D^D*^ (top) and e~^e~ 
D*+D*- from Belle [61]. The dashed lines correspond to the ?/'(4040), ^/'(4160), 
Y(4260), and ^(4415). 
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Figure 8.3: The exclusive cross sections for e~^e~ DD* (top) and e"'"e~ — >■ 
D*D*, after applying radiative corrections, for the work presented in this dis- 
sertation. The dashed lines correspond to the ■0(4040), '?/'(4160), Y(4260), and 
^(4415). There is a factor of two difference between these results and those of 
Belle in Fig. EM 
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Table 8.2: Percentages of event types present at center-of-mass energy of 4160 
MeV for this analysis as compared to the T. Branes predictions |28j . 



Center-of-Mass 
Energy (MeV) 


DD 


D*D 


D*D* 






4160 (This Analysis) 


3.9 ±0.5 


28.2 ± 1.8 


54.1 ±3.6 




9.7 ±0.6 


4159 (Barnes) 


21.6 


0.5 


47.3 


10.8 


18.9 



in qualitative agreement with Eichten's [27j prediction, Figs. Il.6l and ll.7l A more 
quantitative comparison with Eichten's model awaits his updated calculations 
using the correct charmed-meson masses. 

Recently, there has been work on the theoretical side in regards to the cross 
sections in the neighborhood of the D*D* threshold. Using very preliminary 
CLEO-c cross section results presented at the Flavor Physics and CP Violation 
Conference (FPCP) [64j, Voloshin et al. ^65n66j, using very general properties of 
amplitudes near a threshold, obtained resonance parameters assuming the cross 
sections between 3970 and 4060 MeV result from a single resonance. In their inves- 
tigation they determined that the preliminary cross sections can not be described 
by a single resonance, but in fact hint at the presence of two, with one of the res- 
onances being narrow. This can be found in their analysis when they struggled 
to describe the data by a single resonance which resulted in a /NDF = 17.6/8. 
They determined that if they ignored the data points for the cross section at 4010 
and 4015 MeV for DD as well as D*D* at 4015 MeV a statistically significant 
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Figure 8.4: Results for the non-strange D-meson cross sections. The error bars 
include both statistical and systematic errors. 
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Figure 8.5: Results for the Dg-meson cross sections. The error bars include both 
statistical and systematic errors. 
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result is obtained, /NDF = 3.0/5.0 Fig. 18. 6[ The reason for the poor fit when 
including all data points is almost solely due to the data point at 4015 MeV in 
DD, where the value for the cross section is at a minimum. They suggest the 
existence of a new narrow resonance to describe this "dip" in the cross section. 
Following the same fit procedure as they described |66], results are updated to 
reflect the cross sections presented in this dissertation. These updated fits, shown 
in Fig. 18.71 still agree with their initial assessment, that is because of the large 
dip in the DD cross section a single resonance is insufficient in describing the 
cross sections near the D*D* threshold. The resonance parameters, with DD 
and D*D* at 4015 MeV left out the fit, are shown in Table [SJl 



Table 8.3: The resulting fit parameters from a single resonance fit to the cross 
sections, presented in this dissertation, around the D*D* threshold. The DD 
and D*D* data points were excluded in the fit resulting in x^/NDF = 3.8/6. If 
all points were used the fit is quite poor, reflected in the resulting a x^/NDF = 
36.3/8, this is entirely because of DD at 4015 MeV. 



Fit parameters 


Result from Fit 


Wo 


4013 ± 4 MeV 


To 


66 ± 8 MeV 


w 


10.4 ± 1.7 MeV 




15.4 ± 1.3 MeV 


r 

^ ee 


1.9 ± 0.7 keV 
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Figure 8.6: The fit results from Ref. [66] wliicli used preliminary CLEO-c re- 
sults. The data points that are excluded in the fit are shown with circles {DD at 
4010 and 4015 MeV and D*D* at 4015 MeV). The fit are to dimensionless rate 
coefficients Rk where Rk oc {Ri (top left) corresponds to DD, R2 (top right) 
corresponds to DgD^, R^ (bottom left) corresponds to DD*, and R4 (bottom 
right) corresponds to D*D*). 
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Figure 8.7: Fit results, assuming a single resonance [66], using the cross sections 
presented in this dissertation. The fit are to dimensionless rate coefficients 
where Rk oc {Ri (top left) corresponds to DD, R2 (top right) corresponds to 
DgDs, i?3 (bottom left) corresponds to DD*, and i?4 (bottom right) corresponds 
to D*D*). 



181 



8.2 Conclusion 

The CLEO-c scan was extremely successful in accomplishing the goals it put forth 
and as a result has provided a wealth of interesting results. 

The location in center-of-mass energy that maximizes the yield of Ds was 
determined by this analysis to be at 4170 MeV. The maximum occurs not with 
DgDs events, but rather with D^Dg events. It is with this information that CLEO- 
c has already been able to make significant progress, on the decay constant and 
hadronic branching fractions of the Dg meson [HJ [67], with the data collected 
during the scan and in subsequent running at 4170 MeV. 

The total charm cross section between 3.97 GeV and 4.26 GeV has been 
measured and presented. These results have been determined using three differ- 
ent techniques, the exclusive D-meson, inclusive D-meson, and inclusive hadron 
method, and are shown in Table 18.41 and graphically in Fig. 18. 8[ The agreement 
between these three different methods is very good and provides a powerful cross 
check on the results. Radiative corrections were determined and applied to the 
inclusive D meson total charm result. The corrected cross sections are consistent 
with previous experimental results [HI [60] and are more precise (see Fig. 17.21) . 

The composition of the total charm cross section, in terms of DD, D*D, D*D*, 
DgDs, D*Ds, DgD*, D*Dn and D*D*it, at all energies that were investigated, 
has been determined and presented. Results are shown in Tables [8.101 and [8.111 in 
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addition to graphically in Figs. I8.4[ I8.5[ and 18. 9[ There is qualitative agreement 
with the potential model put forth by Eichten et al. [27| for all modes excluding 
the multi-body production, which was not predicted. The existence of a peak 
in the D*D and DgDs channels at the D*D* threshold, along with the observa- 
tion that there is a minimum in the DD, can be interpreted as a possible new 
narrow resonance [66] • Unfortunately, the available data sample is insufficient to 
investigate this latter idea further. 

The D*D* cross section shows a plateau region above its threshold, whereas 
the D*D peaks right at the threshold, in agreement with the recently presented 
results by the Belle Collaboration [61] . 

Also, events containing D mesons plus an additional pion, referred to as multi- 
body, have been shown to exists in the energy region investigated. Figs. 15. 4[ 15. 5[ 
and 15.61 show the cross sections that have been determined. 

Using the data sample collected at 4260 MeV we searched for insight into 

the nature of the F(4260) state. By applying the techniques developed with the 

larger data sample collected at 4170 MeV we find that, besides the presence of the 

D*D*TT final state and larger overall production of multi-body, there is essentially 

nothing to distinguish the charm production at this energy from that at our 

lower-energy points. Although studies are still ongoing, we have preliminarily 

concluded that our data sample is insufficient to determine the contribution of 

states including excited charmed meson states (e.g. DiD) that are expected under 
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some y(4260) interpretations |68j . 

As for the future, it is unlikely that time will be available to allow further 
investigation of the Y(4260) and region around the D*D* to be done at CLEO- 
c. However, the next stage in the BES experiment, BES-III [51], is planning 
to investigate in-depth the region around D*D* and should potentially be able 
to add insight to this region of interest. In addition to BES-III, both Belle and 
BaBar have shown that they too can investigate this region and we are impatiently 
awaiting their updated results. 



Table 8.4: Comparison of the total charm cross section using the 3 different 
methods. First error is statistical and the second systematic. 



Energy 


Exclusive 


Inclusive 


Hadron 


(MeV) 


D-meson (nb) 


D-meson (nb) 


Counting (nb) 


3970 


4.83 ±0.19 ±0.15 


4.91 ±0.18 ±0.16 


4.91 ±0.13 ±0.30 


3990 


5.85 ±0.23 ±0.19 


5.93 ±0.21 ±0.19 


5.87 ±0.14 ±0.34 


4010 


7.10 ±0.14 ±0.23 


7.05 ±0.17 ±0.23 


7.21 ±0.12 ±0.40 


4015 


7.94 ±0.41 ±0.26 


7.62 ±0.34 ±0.25 


7.88 ±0.18 ±0.43 


4030 


10.60 ±0.34 ±0.27 


10.87 ±0.28 ±0.37 


11.30 ±0.15 ±0.59 


4060 


10.16 ±0.36 ±0.27 


9.98 ±0.26 ±0.34 


9.98 ±0.14 ±0.53 


4120 


8.95 ±0.37 ±0.25 


9.13 ±0.28 ±0.31 


9.43 ±0.15 ±0.49 


4140 


9.22 ±0.29 ±0.26 


9.11 ±0.22 ±0.30 


9.58 ±0.24 ±0.50 


4160 


9.33 ±0.20 ±0.26 


9.10 ±0.15 ±0.30 


9.62 ±0.17 ±0.50 


4170 


9.03 ±0.04 ±0.25 


9.09 ±0.07 ±0.30 


9.45 ±0.09 ±0.49 


4180 


8.67 ±0.27 ±0.24 


8.70 ±0.20 ±0.29 


9.07 ±0.12 ±0.47 


4200 


7.42 ±0.35 ±0.20 


7.45 ±0.26 ±0.25 


8.37 ±0.14 ±0.43 


4260 


4.27 ±0.16 ±0.14 


4.20 ±0.10 ±0.14 


4.34 ±0.16 ±0.23 
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Figure 8.8: Comparison of the total charm cross section using the 3 different 
methods. This is a graphical representation of Table 18. 4[ 
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Figure 8.9: Results for the multi-body cross sections. The error bars include both 
statistical and systematic errors. 
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Table 8.5: Efficiencies (units of 10 ^) at each scan-energy point for selection of 
through various decay modes in the three exclusive event types. 



Mode 


E„„ MeV 






D*D* 




3970 


56.2 ± 1.7 


- ± - 


- ± - 




3990 


55.2 ± 1.7 


- ± - 


- ± - 




4010 


57.0 ± 1.7 


- ± - 


- ± - 




4015 


57.8 ± 1.4 


- ± - 


- ± - 




4030 


58.3 ± 1.5 


- ± - 


- ± - 




4060 


56.7 ± 1.5 


- ± - 


- ± - 




4120 


57.6 ± 1.9 


55.8 ± 1.3 


- ± - 




4140 


50.6 ± 1.4 


58.2 ±0.9 


- ± - 




4160 


54.3 ± 1.9 


55.9 ± 1.3 


- ± - 




4180 


52.1 ± 1.9 


55.3 ± 1.3 


- ± - 




4200 


52.9 ± 1.9 


53.9 ± 1.4 


- ± - 




4260 


62.0 ±4.4 


58.1 ±3.1 


56.4 ± 1.5 


7777+ 


3970 


39.2 ±2.3 


- ± - 


- ± - 




3990 


39.3 ± 2.3 


- ± - 


- ± - 




4010 


37.9 ±2.2 


- ± - 


- ± - 




4015 


42.0 ± 1.9 


- ± - 


- ± - 




4030 


53.7 ±2.1 


- ± - 


- ± - 




4060 


39.7 ±2.0 


- ± - 


- ± - 




4120 


33.1 ±2.4 


38.3 ± 1.8 


- ± - 




4140 


37.1 ± 1.8 


40.7 ± 1.3 


- ± - 




4160 


26.6 ±2.3 


38.3 ± 1.8 


- ± - 




4180 


37.9 ±2.5 


38.8 ± 1.8 


- ± - 




4200 


41.1 ±2.6 


34.9 ± 1.8 


- ± - 




4260 


42.6 ±6.0 


39.2 ±4.2 


40.0 ±2.0 




3970 


42.7 ± 1.4 


- ± - 


- ± - 




3990 


42.1 ± 1.4 


- ± - 


- ± - 




4010 


42.9 ± 1.4 


-±- 


-±- 




4015 


41.0 ± 1.1 


-±- 


-±- 




4030 


43.8 ± 1.3 


- ± - 


- ± - 




4060 


43.1 ± 1.3 


- ± - 


- ± - 




4120 


39.6 ± 1.6 


40.5 ± 1.1 


- ± - 




4140 


39.8 ± 1.1 


41.8 ±0.8 


-±- 




4160 


41.6 ± 1.6 


41.3 ± 1.1 


-±- 




4180 


36.7 ± 1.5 


42.8 ± 1.1 


- ± - 




4200 


42.7 ± 1.6 


42.2 ± 1.1 


- ± - 




4260 


43.8 ±3.7 


43.4 ±2.6 


41.8 ± 1.2 
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Mode 


Ec™ MeV 








K*K+ 


3970 


39.1 ± 1.2 


- ± - 


- ± - 




3990 


40.0 ± 1.2 


- ± - 


- ± - 




4010 


42.5 ± 1.2 


- ± - 


- ± - 




4015 


41.3 ± 1.0 


1 

— ± — 


1 

— ± — 




4030 


42.0 ± 1.1 


- ± - 


- ± - 




4060 


41.5 ± 1.1 


- ± - 


- ± - 




4120 


40.1 ± 1.4 


40.8 ± 1.0 


- ± - 




4140 


39.2 ± 1.0 


42.0 ± 0.7 


- ± - 




4160 


42.7 ± 1.4 


42.0 ± 1.0 


- ± - 




4180 


43.3 ± 1.4 


42.6 ± 1.0 


- ± - 




4200 


37.7 ± 1.4 


41.6 ± 1.0 


- ± - 




4260 


41.2 ± 3.2 


35.0 ± 2.2 


41.4 ± 1.1 




3970 


19.2 ± 0.7 


- ± - 


- ± - 




3990 


17.1 ± 0.7 


- ± - 


- ± - 




4010 


16.8 ± 0.7 


- ± - 


- ± - 




4015 


16.9 ± 0.6 


— ± — 


— ± — 




4030 


24.0 ± 0.7 


- ± - 


- ± - 




4060 


16.5 ± 0.6 


- ± - 


- ± - 




4120 


16.7 ± 0.8 


16.9 ± 0.5 


- ± - 




4140 


17.1 ± 0.6 


16.0 ± 0.4 


- ± - 




4160 


14.9 ± 0.7 


16.1 ± 0.5 


- ± - 




4180 


16.3 ± 0.8 


16.2 ± 0.5 


- ± - 




4200 


16.9 ± 0.8 


17.1 ± 0.6 


- ± - 




4260 


15.4 ± 1.8 


18.7 ± 1.3 


17.4 ± 0.6 




3970 


31.9 ± 2.1 


- ± - 


- ± - 




3990 


33.2 ± 2.2 


- ± - 


- ± - 




4010 


30.9 ± 2.1 


- ± - 


- ± - 




4015 


/I 1 "1 

34.8 ± 1.8 


1 

— ± — 


1 

— ± — 




4030 


32.4 ± 1.9 


-± - 


-±- 




4060 


26.7 ± 1.8 


-± - 


-±- 




4120 


21.2 ±2.1 


28.8 ± 1.7 


- ± - 




4140 


27.6 ± 1.6 


25.8 ± 1.1 


- ± - 




4160 


29.2 ±2.4 


29.9 ± 1.7 


- ± - 




4180 


29.9 ±2.4 


32.8 ± 1.7 


-±- 




4200 


29.5 ±2.4 


29.8 ± 1.7 


-±- 




4260 


16.1 ±4.5 


35.0 ±4.1 


28.2 ± 1.8 
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Mode 


Ecm MeV 








r)'p+ 


3970 


16.0 ± 1.1 


- ± - 


- ± - 




3990 


12.9 ± 1.0 


- ± - 


- ± - 




4010 


13.4 ± 1.0 


- ± - 


- ± - 




4015 


13.7 ±0.8 


- ± - 


- ± - 




4030 


13.0 ±0.9 


-±- 


-±- 




4060 


11.6 ±0.8 


-±- 


-±- 




4120 


11.7± 1.0 


14.2 ±0.8 


- ± - 




4140 


12.0 ±0.7 


12.5 ±0.5 


- ± - 




4160 


11.1 ± 1.0 


11.8 ±0.7 


- ± - 




4180 


15.0 ± 1.2 


12.1 ±0.7 


-±- 




4200 


14.0 ± 1.1 


13.3 ±0.8 


-±- 




4260 


14.2 ±2.0 


14.9 ± 1.9 


15.0 ±0.9 




3970 


13.6 ±0.7 


- ± - 


- ± - 




3990 


13.4 ±0.7 


- ± - 


- ± - 




4010 


13.4 ±0.7 


- ± - 


- ± - 




4015 


12.9 ±0.6 


-±- 


-±- 




4030 


15.5 ±0.7 


- ± - 


- ± - 




4060 


12.6 ±0.6 


- ± - 


- ± - 




4120 


13.7 ±0.8 


12.5 ±0.6 


- ± - 




4140 


14.2 ±0.6 


11.8 ±0.4 


- ± - 




4160 


14.9 ±0.8 


12.3 ±0.5 


-±- 




4180 


15.6 ±0.9 


13.4 ±0.6 


-±- 




4200 


15.8 ±0.9 


13.6 ±0.6 


- ± - 




4260 


15.2 ±2.0 


15.5 ± 1.4 


11.0 ±0.6 



Table 8.6: Numbers of events at — 3671 MeV for various event types which 
are used for continuum subtraction in the Hadron Counting Method to determine 
the total charm cross section. 



Off-Res Data 


244663 ± 495 


B.W. Tail of i/j{2S) 


7201 ± 493 


Radiative Return to J/ ip 


10240 ± 157 




9227 ±55 




3058 ± 649 




184 ±21 


Efficiency Corrected Num. of qq 


378420 ± 1934 
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Table 8.7: Numbers of events at the scan energies for various event types which 
are used in the Hadron Counting Method to determine the total charm cross 
section. 



Ecm = 3970 MeV 




Number of Raw Hadronic Events 

Scaled qq 
Radiative Return to ip{2S) 
Radiative Return to J/ip 
Radiative Return to ■0(3770) 

Number of Beam Gas Events 


60811 ± 247 
38041 ± 233 
2352 ± 11 
1059 ± 8 
332 ± 1 
2915 ± 98 
340 ± 133 

32 ±5 
439 ± 102 


Efficiency Corrected Num. 


18926 ± 485 


E^rn = 3990 MeV 




Number of Raw Hadronic Events 
Scaled qq 
Radiative Return to il'{2S) 
Radiative Return to J/ip 
Radiative Return to ■0(3770) 

Number of Beam Gas Events 


54818 ± 234 
32649 ± 201 
1869 ± 9 
889 ±6 
264 ± 1 
2365 ± 84 
471 ± 140 

28 ±5 
341 ± 122 


Efficiency Corrected Num. 


19696 ± 458 


Ecm = 4010 MeV 




Number of Raw Hadronic Events 
Scaled qq 
Radiative Return to iIj{2S) 
Radiative Return to J/ip 
Radiative Return to ■0(3770) 
r+r~ 
e+e~ 

Number of Beam Gas Events 


97629 ± 312 
54635 ± 327 
2929 ± 13 
1429 ± 10 
408 ±2 
4128 ± 144 
635 ± 212 

46 ±8 
565 ± 173 


Efficiency Corrected Num. 


40539 ± 679 
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Ecm = 4015 McV 




Number of Raw Hadronic Events 


26436 ± 163 


Scaled qq 


14240 ± 96 


Radiative Return to ip{2S) 


753 ±4 


Radiative Return to J/ip 


371 ±3 


Radiative Return to '^(3770) 


106 ± 1 




1137 ±39 




165 ± 55 


r r 


12 ±2 


Number of Beam Gas Events 


152 ±87 


Efficiency Corrected Num. 


11581 ± 267 


E = 4030 MeV 




Number of Raw Hadronic Events 


62354 ± 250 


Scaled qq 


29017 ± 180 


Radiative Return to ip{2S) 


1460 ± 7 


Radiative Return to J lib 

1 1 


733 ±5 


Radiative Return to ■?/'(3770) 


197 ± 1 




2252 ± 78 




593 ± 145 




25 ±4 


Number of Beam Gas Events 


285 ± 127 


Efficiency Corrected Num. 


33955 ± 457 


E.^ = 4060 MeV 




Number of Raw Hadronic Events 


64117 ±253 


Scaled qq 


31305 ± 193 


Radiative Return to ip{2S) 


1467 ± 7 


Radiative Return to J/ip 


797 ±6 


Radiative Return to ■0(3770) 


194 ± 1 




2582 ± 87 




277 ± 108 




26 ±4 


Number of Beam Gas Events 


345 ± 139 


Efficiency Corrected Num. 


32803 ± 455 



195 



Ecm = 4120 MeV 




Number of Raw Hadronic Events 


51863 ± 228 


Scaled qq 


25694 ± 161 


Radiative Return to ip{2S) 


1027 ±5 


Radiative Return to J/ip 


552 ±4 


Radiative Return to ■0(3770) 


126 ± 1 




2210 ± 74 




462 ± 118 




22 ±4 


Number of Beam Gas Events 


280 ± 131 


Efficiency Corrected Num. 


26016 ±411 


E = 4140 MeV 




Number of Raw Hadronic Events 


91780 ± 303 


Scaled qq 


45202 ± 220 


Radiative Return to ipilS) 


1705 ± 7 


Radiative Return to J/ip 


954 ±6 


Radiative Return to '0(3770) 


221 ± 1 




3585 ± 100 




739 ± 162 




38 ±5 


Number of Beam Gas Events 


518 ± 163 


Efficiency Corrected Num. 


46910 ± 548 


E.^ = 4160 MeV 




Number of Raw Hadronic Events 


190764 ± 437 


Scaled qq 


93833 ± 398 


Radiative Return to ip{2S) 


3467 ± 11 


Radiative Return to J/ip 


1876 ± 11 


Radiative Return to -0(3770) 


400 ±2 




7761 ± 192 


e+e" 


1417 ±301 




78 ±9 


Number of Beam Gas Events 


1019 ±212 


Efficiency Corrected Num. 


97750 ± 879 
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Ecm = 4170 McV 




Number of Raw Hadronic Events 
Scaled qq 
Radiative Return to ip{2S) 
Radiative Return to J/ip 
Radiative Return to V'(3770) 

Number of Beam Gas Events 


3301243 ± 1817 
1641982 ± 9325 
59506 ± 239 
33037 ±252 

7001 ± 28 
138524 ± 4746 
24731 ± 7419 
1361 ±227 
5256.47 ±205.89 


Efficiency Corrected Num. 


1683110 ± 15884 


E = 4180 MeV 




Number of Raw Hadronic Events 
Scaled qq 
Radiative Return to ip{2S) 
Radiative Return to J/ip 
Radiative Return to ■?/'(3770) 

Number of Beam Gas Events 


103834 ± 322 
52037 ±310 
1806 ± 8 
1009 ± 8 
221 ± 1 
4744 ± 156 
840 ± 234 

43 ±7 
613 ± 131 


Efficiency Corrected Num. 


51410 ± 662 


E.^ = 4200 MeV 




Number of Raw Hadronic Events 
Scaled qq 
Radiative Return to ip{2S) 
Radiative Return to J/ip 
Radiative Return to ■0(3770) 

e+e" 

Number of Beam Gas Events 


49296 ± 222 
25592 ± 161 
852 ±4 
470 ±4 
109 ± 1 
2159 ±75 
291 ± 97 

21 ±4 
329 ± 131 


Efficiency Corrected Num. 


23509 ± 397 
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= 4260 MeV 




Number of Raw Hadronic Events 


183192 ± 427 


Scaled qq 


117111 ± 520 


fVdQldtlVt! rtcLUIIl i/y^^ZO J 


Q/l 77 -1- 1 1 


Radiative Return to J/ip 


1811 ± 12 


Radiative Return to ijj{3770) 


424 ±2 




10239 ± 267 




2178 ± 426 




97± 12 


Number of Beam Gas Events 


906 ± 262 


Efficiency Corrected Num. 


56787 ± 1066 
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P (GeV/c) 

Figure 8.10: Sideband-subtracted momentum spectrum for K^tt^ at 3970 

MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.11: Sideband-subtracted momentum spectrum for K^tt^tt^ at 

3970 MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.12: Sideband-subtracted momentum spectrum for ^ K vr^ at 3990 
MeV. The data is the points with error bars and the histograms are MC. 
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P (GeV/c) 

Figure 8.13: Sideband-subtracted momentum spectrum for K^tt^tt^ at 

3990 MeV. The data is the points with error bars and the histograms are MC. 
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P (GeV/c) 

Figure 8.14: Sideband-subtracted momentum spectrum for K^tt^ at 4010 

MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.15: Sideband-subtracted momentum spectrum for K 7r+7r+ at 

4010 MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.16: Sideband-subtracted momentum spectrum for ^ K 7r+ at 4015 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.17: Sideband-subtracted momentum spectrum for K 7r+7r+ at 

4015 MeV. The data is the points with error bars and the histograms are MC. 
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P (GeV/c) 

Figure 8.18: Sideband-subtracted momentum spectrum for K^tt^ at 4030 

MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.19: Sideband-subtracted momentum spectrum for K 7r+7r+ at 

4030 MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.20: Sideband-subtracted momentum spectrum for ^ K vr"'' at 4060 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.21: Sideband-subtracted momentum spectrum for K 7r+7r+ at 

4060 MeV. The data is the points with error bars and the histograms are MC. 
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P (GeV/c) 

Figure 8.22: Sideband-subtracted momentum spectrum for K''tt~^ at 4120 

MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.23: Sideband-subtracted momentum spectrum for K^tt^tt^ at 

4120 MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.24: Sideband-subtracted momentum spectrum for ^ K 7r+ at 4140 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.25: Sideband-subtracted momentum spectrum for K 7r+7r+ at 

4140 MeV. The data is the points with error bars and the histograms are MC. 
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MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.27: Sideband-subtracted momentum spectrum for K^tt^tt^ at 

4160 MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.28: Sideband-subtracted momentum spectrum for K''tt~^ at 4180 

MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.29: Sideband-subtracted momentum spectrum for K 7r+7r+ at 

4180 MeV. The data is the points with error bars and the histograms are MC. 



218 




P (GeV/c) 

Figure 8.30: Sideband-subtracted momentum spectrum for K^tt^ at 4200 

MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.31: Sideband-subtracted momentum spectrum for K 7r+7r+ at 

4200 MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.32: Sideband-subtracted momentum spectrum for ^ K 7r+ at 4260 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.33: Sideband-subtracted momentum spectrum for K^tt^tt^ at 

4260 MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.34: Sideband-subtracted momentum spectrum for — ^ 07r+ at 3970 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.35: Sideband-subtracted momentum spectrum for — ^ 07r+ at 3990 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.36: Sideband-subtracted momentum spectrum for — ^ 07r+ at 4010 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.37: Sideband-subtracted momentum spectrum for — ^ 07r+ at 4015 
MeV. The data is the points with error bars and the histograms are MC. 



226 



4.5 

4t 
3.5 

> 

o 3 

00 

« 2.5 
c 

3 

1.5 
1 

0.5 

a 



Total fit 

Dashed MC events with E,^>1 MeV 
Solid MC events with E,<1 MeV 



I ' ' ■ ■ 



0.2 




0.4 0.6 
P (GeV/c) 



0.8 



Figure 8.38: Sideband-subtracted momentum spectrum for — ^ 07r+ at 4030 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.39: Sideband-subtracted momentum spectrum for — ^ 07r+ at 4060 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.40: Sideband-subtracted momentum spectrum for — ^ 07r+ at 4120 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.41: Sideband-subtracted momentum spectrum for — ^ 07?"*" at 4140 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.42: Sideband-subtracted momentum spectrum for — ^ 07r+ at 4160 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.43: Sideband-subtracted momentum spectrum for — ^ 07r+ at 4180 
MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.44: Sideband-subtracted momentum spectrum for 07r+ at 4200 

MeV. The data is the points with error bars and the histograms are MC. 
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Figure 8.45: Sideband-subtracted momentum spectrum for 07r+ at 4260 

MeV. The data is the points with error bars and the histograms are MC. 
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Table 8.8: Observed cross sections for various center of mass energies as determined by the cut and count method 
described in the text. All error are statistical and cross sections shown are the weighted sums or averages for the 
event types. 



Cross Section 


3970 MeV 


3990 MeV 


4010 MeV 


4015 MeV 


4030 MeV 


4060 MeV 




3.85 pb-^ 


3.36 pb-i 


5.63 pb-^ 


1.47 pb'i 


C= 3.01 pb-i 


3.29 pb-i 


a{DsDs) nb 


0.102 ±0.026 


0.133 ±0.031 


0.269 ± 0.030 


0.250 ±0.059 


0.174 ±0.036 


0.051 ± 0.028 


(t{D*Ds) nb 


-± - 


-±- 


-±- 


-±- 


- ± - 


- ± - 


a{D*D*) nb 


- ± - 


- ± - 


- ± - 


- ± - 


- ± - 


- ± - 


a{DD) nb 


0.194 ±0.027 


0.146 ±0.027 


0.085 ± 0.020 


0.037 ±0.033 


0.410 ±0.037 


0.615 ±0.041 


a{D*D) nb 


4.289 ±0.117 


5.105 ±0.136 


6.285 ±0.116 


6.999 ±0.236 


6.200 ±0.163 


3.687 ±0.129 


a{D*D*) nb 


-± - 


-±- 


-±- 


0.153 ±0.035 


3.745 ±0.151 


4.779 ±0.189 


(T{Charm) nb 


4.585 ±0.123 


5.384 ±0.142 


6.639 ±0.121 


7.439 ± 0.248 


10.529 ±0.228 


9.132 ± 0.234 



Table 8.9: Observed cross sections for various center of mass energies as determined by the cut and count method 
described in the text. All error are statistical and cross sections shown are the weighted sums or averages for the 
event types. 



Cross Section 


4120 MeV 


4140 MeV 


4160 MeV 


4180 MeV 


4200 MeV 


4260 MeV 




C= 2.76 pb"^ 


£= 4.87 pb-i 


C= 10.16 pb-i 


C= 5.67 pb-i 


C= 2.81 pb-i 


C= 13.11 pb-^ 


a{D,D,) nb 


0.026 ±0.026 


0.025 ±0.020 


-0.008 ±0.012 


0.007 ±0.016 


0.015 ±0.022 


0.034 ±0.009 


o{D*D,) nb 


0.478 ± 0.064 


0.684 ±0.059 


0.905 ± 0.044 


0.889 ±0.059 


0.812 ±0.082 


0.047 ±0.022 


a{D:D:) nb 


- ± - 


- ± - 


- ± - 


- ± - 


- ± - 


0.440 ± 0.027 


a{DD) nb 


0.466 ± 0.039 


0.328 ±0.026 


0.267 ±0.017 


0.249 ± 0.022 


0.295 ± 0.032 


0.105 ±0.011 


a{D*D) nb 


2.481 ±0.122 


2.243 ±0.090 


2.100 ±0.061 


1.951 ± 0.079 


1.615 ±0.107 


1.407 ±0.047 


(t{D*D*) nb 


5.073 ±0.213 


4.868 ±0.165 


5.153 ±0.117 


4.620 ±0.147 


3.730 ±0.197 


1.142 ±0.059 


aiCharm) nb 


8.524 ± 0.258 


8.148 ±0.200 


8.417 ±0.141 


7.716 ±0.179 


6.467 ± 0.242 


3.175 ± 0.084 



Table 8.10: Observed cross sections for various center of mass energies as determined by fits to the and 
momentum spectra. The Dg cross sections are determined using a weighted sum techneque discribed in Sect. 14.11 
The first error is statistical and the second is systematic. 



Cross Section 


3970MeV 


3990MeV 


4010MeV 




/L- — O.oOpU 


— o.oupu 


>L- — O.UOpL) 


rr( f) n nh 




u.ioo HI u.uoi m u.uu/ 


U.ZUy HI U.UOU HI U.UlO 


rr( D* f) \ nVi 

ayu^Ug) no 


-1- -1- 

— IL — IL — 


-L -L 

— m — m — 


-L -L 
— HI — HI — 


rr( D* n*"! nVi 


-L -L 

— m — m — 


-I- — + — 


— -1- — -1- 

— HI — HI — 


[lyiy j liu 


fi ooQ _|_ fi noQ _|_ fi m n 


fi OOQ _|_ n _|_ fi m n 


fi 91 1 -u n o'i'i -u n nriQ 

HI u.uoo HI u.uuy 


rr( n* D\ nVi 
U LJ j liU 




^ ZLQD -1- n 991 -\- r\ ^R7 


« «On _L fl 1 QQ _L fl OOC 
U.U.^U HI U.IOO HI yj.ZiZiO 


u\^j-y u J nu 


-L -L 

— m — m — 


-L -L 
— HI — HI — 


— m — HI — 


U\^U UTl J UU 


-L -L 

— m — m — 


-L -L 
— HI — HI — 


-L -L 
— HI — HI — 




— m — m — 


-1- — -1- — 
— HI — HI — 


_ 4- _ -1- _ 
— m — HI — 


U\\^ILUjI lllj liU 


^.000 HI u.iyo m u.io^ 


^ RAf, -\- n 998 -1- n 1 87 

tJ.O^U HI U.^^O HI U.lOl 


7 1 nn -1- n 1 zld -1- n 99fi 

1 .lUU HI U.l^U HI \J.ZiZi\J 


Cross Section 


4015MeV 


4030MeV 


4060MeV 




£= 1.47pb-i 


C= 3.01pb-i 


C= 3.29pb-i 


a(D,Z},) nb 


0.250 ±0.059 ±0.014 


0.174 ±0.036 ±0.010 


0.051 ±0.028 ±0.003 


(r{DlDs) nb 


- ± - ± - 


- ± - ± - 


- ± - ± - 


a(D:Z}:) nb 


- ± - ± - 


- ± - ± - 


- ± - ± - 


(t(L'D) nb 


0.038 ±0.020 ±0.002 


0.530 ±0.078 ±0.024 


0.890 ±0.091 ±0.040 


a{D*D) nb 


7.443 ± 0.394 ± 0.256 


6.500 ±0.257 ±0.221 


4.400 ±0.205 ±0.150 


a{D*D*) nb 


0.213 ±0.076 ±0.009 


3.400 ±0.211 ±0.160 


4.680 ±0.258 ±0.220 


a{D*D7i) nbw 


- ± - ± - 


- ± - ± - 


0.144 ±0.094 ±0.017 


a{D*D*7T) nb 


- ± - ± - 


- ± - ± - 


- ± - ± - 


a (Charm) nb 


7.941 ±0.406 ±0.257 


10.604 ±0.344 ±0.274 


10.165 ±0.356 ±0.270 



Table 8.11: Observed cross sections for various center of mass energies as determined by fits to the and 
momentum spectra. The Dg cross sections are determined using a weighted sum techneque discribed in Sect. 14.11 
The first error is statistical and the second is systematic. 



Cross Section 


4120MeV 


4140MeV 


4160MeV 


4170MeV 




/L- — Z. ( UpU 


JL, — ^.0 ( pU 


— iu.iupu 


JL, — i ( o.ypu 




fi n9fi -1- n n9fi -i- n nm 


U.UZO ZL U.UZU ZL U.UUl 


— n nns -1- n m 9 -1- n nnn 

— u.uuo ZL u.uiz ZL u.uuu 


u.uo^ Zl U.uuo zl u.uuz 


rr( D* f) ^ nl-i 

a[UgUg) no 




u.uo^ ZL u.uoy ZL U.UOU 


u.yuo ZL u.uii ZL u.u^o 


fiQiR-ufinii -i-n n4Q 

u.yiu ZL u.uii ZL u.u^y 




— -1- — -1- 

— ziz — ziz — 


-L -L 

— ZL — ZL — 


— zn — zn — 


-L -L 

— ZL — ZL — 




fi R1 Q -u n nsfi -1- n n98 

U.UIO ZL U.UOU ZL \J.\JZjO 


fi Q77 _|_ fi fizLQ -1- n m 7 

u.o 1 1 ZL u.u^y ZL u.ui ( 


fl Q«7 _|_ fl f|Qc; _|_ fl m 7 
u.ou 1 ZL u.uoo ZL u.ui / 


fl Q«fi _i_fi fiin-i-n nifi 

U.OUU ZL U.UIU ZL u.uiu 


rr( D* ^^ nh 


9 Q«n _|_ fi 1 Q9 _|_ fi 1 fiT 


9 79« _|_ fi 1 zLD -1- n DQ'? 
z. 1 zu ZL u.i^u ZL u.uyo 


9 «90 _|_ fl flQO _|_ fl f|OQ 

z.uzo ZL u.uyo ZL u.uoy 


9 c;«f| _|_ fl fi'^n -1- n (187 

Z.OUU ZL U.UOU ZL U.UOi 


(J\LJ U J llL) 


/I CQfi _L n 970 _|_ fi 997 


A QQQ -u n 998 -1- n 9^?^ 
^.yyy zl u.zzo zl u.zoo 


nzLp; -1- n i -1- n 9'?7 
o.u^o ZL u.ioy ZL u.zo( 


A 79fi _L fl fi'^n -1- n 999 

^. 1 ZU ZL U.UOU ZL u.zzz 


rr( D* D-tt') nVi 
(71 _Ly7l J llL) 


u.u^o ZL u.uoo ZL u.uuo 


n zii 9 -u n ns7 -1- n dzlq 

U.^lZ ZL U.UO ( Zl u.u^y 


fl QCQ -U fl fifin -1- n 047 
u.ooy ZL u.uuu ZL u.u^/ 


u.^^u ZL u.uzu ZL u.uoo 


n( Tl* r)*'Tr\ nVi 
u I ly Lj i\ j liu 


— -1- — -1- — 

— ZL — ZL — 


— -1- -1- — 

— ZL — ZL — 


_ -1- _ -1- _ 

— ZL — ZL — 


_ -U _ -1- _ 

— ZL — ZL — 


CIO 1 itil 1 11 j 11 U 


qp:9 _|_ fi oec i n 9p:i 
ZL U.OUU ZL Vj.^OL 


Q 99Q _|_ fi 9Q9 _|_ fi 9«n 
c/.Zi^o ZL u..6yz ZL U..6UU 


q Q9« _|_ fl 9fin -1- n 9fi'? 
y.ozu ZL u.zuu ZL u.zuo 


Q fl'? -1- n 044 -1- n 9^0 

y.uo ZL u.u^^ ZL u.zou 


Cross Section 


4180MeV 


4200MeV 


4260MeV 






C= 5.67pb-i 


£= 2.81pb-i 


C= 13.11pb-i 




a[DgDg) nb 


0.007 ±0.016 ±0.000 


0.015 ±0.022 ±0.001 


0.047 ±0.022 ±0.003 




(T{D*Dg) nb 


0.889 ±0.059 ±0.047 


0.812 ±0.082 ±0.043 


0.034 ±0.009 ±0.002 




a{D:D:) nb 


- ± - ± - 


- ± - ± - 


0.440 ±0.027 ±0.030 




a{DD) nb 


0.376 ±0.047 ±0.017 


0.361 ±0.066 ±0.016 


0.180 ±0.022 ±0.008 




a{D*D) nb 


2.507 ±0.127 ±0.085 


2.100 ±0.169 ±0.071 


2.102 ±0.080 ±0.071 




a{D*D*) nb 


4.318 ±0.201 ±0.203 


3.394 ±0.256 ±0.160 


0.506 ±0.069 ±0.024 




a{D*Dn) nb 


0.575 ±0.092 ±0.069 


0.735 ±0.129 ±0.088 


0.638 ±0.093 ±0.077 




a{D*D*7T) nb 


- ± - ± - 


- ± - ± - 


0.322 ±0.067 ±0.080 




a{Charm) nb 


8.672 ±0.266 ±0.236 


7.417 ±0.350 ±0.201 


4.269 ±0.162 ±0.138 





Appendix A 

Helicity Formalism and Angular 
Distributions 



In the spirit of obtaining the best possible MC sample a considerable amount of 
time was spent on determining the final-state angular distributions for each of 
the event types and the subsequent decays of the starred states. There are three 
different and distinct types of events that need to be considered: 

• 1- ^ O-Q-, as in e+e" ^ -f* ^ DD 

• 1- ^ 0-1-, as in e+e" ^ 7* ^ DD* 

• I- ^ as in e+e" ^ 7* ^ D*D* 

In addition to the above, there are also the subsequent decays: 
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• 1- ^ 1-0", as in D* -fD 

• 1" O^O", as in D* nD 

Each of the above cases, except 1~ 0~0~, since it is already understood 
from '?/'(3770) running (for completeness its distribution is sin^6'), is described 
individually below using the helicity formalism as discussed by J.D. Richman in 
Ref. HO] 

A.l e+e- ^ 7* ^ DD* 

The quantum numbers, as stated above, are 1^ — > 0^1^. Since the reaction in 
question, and all reactions that will be discussed henceforth, is an electromagnetic 
interaction, parity must be conserved. As proved in Ref. [TOj, parity conservation 
can be stated as follows: 

^A^A^, = i-iy-' PrPoPD'A^Xo-x^. (A.l) 

By Eq. lA.ll and the quantum numbers of the reaction, it is seen that party 
conservation requires Aqq = —Aqq and Aiq = —A_iq. Therefore, Aqq must be set 
to zero and so the D* will always be produced with helicity of A = ±1. Next, the 
virtual photon will have a spin projection along the beam axis. This is because 
when Ecm » electrons and positrons will only couple to each other if there 
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combined helicity is odd. Using the above information we get the following final- 
state angular distribution: 

dN 



J2 \A^o\MdI,^-i\' + \Dlj^i\' - Dl,^,D]C,^_, - Dl,^_,D]^^,) (A.2) 

where Dl^^^ = e-^^^'<^d\i^^{e)e'^ and D\j^_^ = e-'^^-'^dlj__^{e)e-''f'. Substituting 
and realizing that only = 1 needs to be looked at, since (i^^/ = d-'_^,_^, one 
arrives at by ignoring overall constants and integrating over azimuth at, 

dN 



dcos6 



oc lAoRl + cos^^), (A.3) 



where 6 is defined as the angle between the D* and the beam axis. 

To generate the correct angular distributions for DD* events, EVTGEN's Helic- 
ity Amplitude model (HELAMP) [57j was used and the following was added to the 
decay.dec file for D^^D^: 

D_s*+ D_s- HELAMP 1. 0. 0. 0. -1. 0. : 



A.2 e+e- ^ 7* ^ DD% D* D-f 

This decay is identical to 7* — DD* except with the 7 and the D* switched. 
Thus the D* will have a spin projection along its momentum vector and the 7 

240 




o 



300- 
200^ 
100^ 

n t_i I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-1 -0.8 -0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 1 

Cos(e) 

Figure A.l: Generator-level angular distribution for D*D events. 9 is defined as 
the angle between the D* and the beam axis. 
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will be produced with helicity of A = ±1. Therefore, with Biq being the helicity 
amplitude for D* D'j 

, oc \B^o\'\A,o\\\dU{9')\' + \dl,{9')\') oc \B^o\'\A,o\'{l + cos' O') (A.4) 



d cos 9 



where 9' is defined as the angle between the D in the rest frame of the D* and 
the momentum vector of the D* in the lab. The full derivation for this decay is 
given in Appendix IA.7I 




Figure A. 2: Generator-level angular distribution for D in D* —>■ D^y for D*D 
events. 9 is defined as the angle between the D in the rest frame of the D* and 
the momentum vector of the D* in the lab. 
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A.3 e+e- ^ 7* ^ DD% D* Dtt 



Since — Xd — the angular distribution is as follows: 



dN 
dcos9 



J oc |AiorMio(^')r oc |5oo||Aio|'(sin2e') 



(A.5) 



where 6* is defined above and Bqq being the helicity amplitude for D* — > D'y. 




-0.8 -0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 

cos(e') 



Figure A.3: Generator-level angular distribution for D in D* ^ Dtt for D*D 
events. 6 is defined as the angle between the D in the rest frame of the D* to 
that of the momentum vector of the D* in the lab. 
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A.4 e+e" ^ 7* ^ D*D* 



There are seven helicity amplitudes that determine the angular distribution of 
the final state: Aqo, v4o_i, v4_io, Aio, Aqi, Ah, and v4_i_i. By the conservation of 
parity and the symmetry of the final state, Aa^.a^, = ^a^.a^. = ^-a^.-a^*, so 
the only independent amplitudes are Aqo, ^oi and An. It therefore follows that 

dN 



dcosO 



l^or(ci}oW)'+2|Air((c^}iW)' + (rfLiW)')+2|Anr(c^loW)'- (A.6) 



The three independent helicity amplitudes have not been measured. For sim- 
plicity, in generating MC we assumed that all helicity amplitudes are equal, 
^00 = ^01 = ^11- 

oc |(l-cos2e) + (l + cos2e) oc (l-^cos^e). (A.7) 



d cos 9 2 5 

The above helicity amplitude assumption is consistent with the Mark I measure- 
ment of e^e^ ^ 7* — i> D*D* angular distribution from 1977 [7T]. However, it 
should be noted that the error bar of that measurement is ~ 100%. Here, the 
error bar refers to the fit parameter a in their fit function, 1 + a cos^ 6. In order 
to generate this angular distribution for D*D* events the HELAMP model was used 
and the following was added to the decay.dec file for D*^D*^: 

D*0 anti-D*0 HELAMP 1.0 0.0 1.0 0.0 1.0 0.0 1.0 0.0 1.0 0.0 1.0 

0.0 1.0 0.0; 
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Figure A. 4: Generator-level angular distribution for D*D* events. is defined at 
the angle between the D* and the beam axis. 
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To check our assumption the hehcity amphtudes have been determined and are 
presented in Appendix B using the scan data above D*D* threshold. However, it 
should be noted that the detection efficiency does not depend on the coefficient 
of cos^ e in Eq. \K1\ 

A. 5 e+e- ^ 7* -> D*D% D* D7 

It has already been shown in a previous section that the angular distribution, 
when D* has helicity of 1, is 

dN 



oc(l + cos^^). (A.8) 



d cos 9' 

Next, it can been shown that when D* has helicity of the angular distribution 
is as follows: 



J oc (|4_i(^')r + \doi{0 )r) cos' 6) = sin'e . (A.9) 



dcos9 

Now when both helicities are allowed the angular distribution will be a combina- 
tion of these two. There are three cases to consider: 

• Both D* are produced with helicity of 1. 

• Both D* are produced with helicity of 0. 

• One of D* is produced with helicity of 1 while the other is produced with 
helicity of 0. 
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Combining all the above situations, and being mindful to keep track of the needed 
helicity amplitudes, leads to the following angular distribution: 

= isoi|2{(3iAor + + i^oor) + (lAir - i^oor - i^ioHcos^^'}, 

a cos 6* 

(A.IO) 

where O' is defined above and Biq being the helicity amplitude for D* — > D7. 
Again making the same assumption as above, that is all the helicity amplitudes 
are equal one arrives at the following: 

"^^.ocll-^cos^^'). (A.ll) 



d cos 9 5 

A.6 e+e- ^ 7* ^ D*D\ D* Dtt 

It has already been shown that the angular distribution, for production of D* 
with helicity of 1 is 

-oc (1-cos^e). (A.12) 



dcosd 

Next, it can been shown that when D* has helicity of the angular distribution 
is as follows: 

-^^\dl,{e')\'^cosH'. (A.13) 

Qj COS U 

Again combining the three possible D*D* hehcity cases, with their associated 
amplitudes, yields the following angular distribution: 

|5oo|'{(|^ii|' + l^oir) + d^ool' + l^oir - lAiiHcos^^'}. (A.14) 
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Figure A. 5: Generator-level angular distribution for in D* — > D7 for D*D* 
events. 9 is defined as the angle between the D in the rest frame of the D* to 
that of the momentum vector of the D* in the lab. 
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Applying the same assumption as above, that is Aqq = Aqi = An, one arrives at 
the following distribution: 

7 OC (l + ^cos^^'). (A. 15) 



d cos 9 2 
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Figure A. 6: Generator-level angular distribution for D in D* ^ Dn^ for D*D* 
events. 6 is defined as the angle between the D in the rest frame of the D* to 
that of the momentum vector of the D* in the lab. 



249 



A. 6.1 Momentum Distributions 

Momentum Distributions of D-mesons in e+e" — > 7* ^ D*D* events 

Notice that in all cases the angular distribution of a D meson can be written as 
follows: 

oc1 + q;cos2(^) (A.16) 



c?cos^ 

where a is a function of the various helicity amplitudes present for the decay and 
event in question. 

The momentum distribution for the decay D mesons is as follows: 

dN dN dcosO 



dPiab dcos9 dPiab 
The energy of the decay D in the lab can be written as follow: 



(A.17) 



Eiab = 7^" + PlP" cos(^) (A.18) 



thus one gets, using P^^ — Ef^j^ — M^, the following expression: 



and so: 

Cic0s(6') Plab 



dPlab (3^Po^Pii, + M^ 



(A.20) 



with p-f = 7 = and E° = ^^^^^ ^ 



Pni. ^ = Mm. and = M^.+Mj-M^ 
mass of the daugther D and other particle (either a photon or pion) present in 
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the decay, respectively. Therefore: 



'''"^ + (A.21) 



It is easily seen in the above expression the role that the angular distribution, 
indicted by the parameter a, has on the momentum distribution. Examples of 
two cases, aX E = 4030 MeV, are shown in Figs. IA.7tiA.8i 



12000 



10000 




-0.05 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Pd (GeV) 



Figure A. 7: Momentum distribution for in D*D* events. The different his- 
tograms correspond to different helicity amplitudes, and therefore different values 
for a in Eq. IA.21[ The solid histogram has a = 0.65 and a = —0.24 for the pion 
and gamma decay respectively. The dashed histogram has a = 5.0 and a = —0.71 
for the pion and gamma decay respectively. 
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Pd (GeV) 

Figure A. 8: This has the same parameters and values as Fig. IA.7I but has been 
convoluted with a gaussian of 10 MeV to account for detector resolution. This 
distribution should be compared to Fig. 18.181 
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A. 7 Complete Derivation of e^e ^ 7* ^ D*D, 
D* L>7 

To help curb confusion let us define e~^e~ — >^ 7* — > D*D, D* D'y as 7* 12 
and 1 — * 34 and define AaiA2 to be the helicity amplitudes for 7* — > 12 and -BA3A4 
to be the helicity amplitudes for 1 34. Therefore 

A(M„ Ai, A2, A3, A4) = J2 Kx'Dm,xBx,x,Ax^x, (A.22) 

Ai 

where A' = A3 — A4 and A = Ai — A2 and since A2 = A4 = 0, A' = A3 and A = Ai 
therefore: 

A{M,, Ai, 0, A3, 0) = Dl,^Dl,_^,Box,Aoo + 

D\,^Dl,^_,Box,A^,o (A.23) 

By Eq. lA.ll and the quantum numbers of the reaction, it is seen that party 
conservation requires Aqq = —Aqq and Aiq = —A_iq. Therefore, Aqq must be set 
to zero and so the D* will always be produced with helicity of A = ±1. So: 

A{M„ Ai, 0, A3, 0) = Box,A,oiDl,^Dl,_^, - D\^^D\,^^,) (A.24) 

The angular distribution is therefore: 

M-y A4 

= E E \Aio?\B,x,?\Dl^,D\,^, - D\^^D\,^_,\^ (A.25) 

A4 
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where Dl,^, = e-^^'^^S^^Me"^ . ^m,-i = e-^-<^di,^_i(^)e-^ 
D\^^ = e-''^^d\^^{e^)e'^^'^^ and D\^^^ = e^'^4c^LiA4(^4)e^^*'^^ where 9 is defined as 
the angle between the D* and the beam axis and 6^ is defined as the angle between 
the D in the rest frame of the D* to that of the momentum vector of the D* in 
the lab. By substitution of the D functions one arrives at: 

-«.(^4)rfL,i(^)^iu.(^4)cii,,_i(^))2cos(250)) (A.26) 

where 50 = ^4 — 0. Next, we sum over A4 where A4 = ±1 since it represents a 
photon for the present case. Therefore |i^oiP = l-^o-il^ meaning: 

^ = iAonsoir(Mi,(^4)rfL,i(^)r+ 

-(^^n(^4)c^L,iW^^iii(^4)c^L^_i(^))2cos(250) + 

\dU{e,)d}M,,{e)\' + 

Hdl^^{e^dl,^Md\-^{04)dl,^-M)'2<^os{2S<f>)) (A.27) 

again since we are not concerned by overall constants we set My = 1 and using 
the fact that dl^ = dl^_^ = i±f^ and d\_^ = dl^^ - i^f^ one arrives at the 
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following: 

^oc \Aiof\Boi\'^{{l + cos^94:){l + cos'^9)+sm^9sm'^94Cos{2d(j))) (A.28) 

therefore by integrating over the azimuthal angle, the D* has a (1 + cos^^) dis- 
tribution with respect to the beam axis. In addition, one sees that the D, in 
D* — > D7,will have a (1 + cos^ ^4) distribution in the rest frame of the D*. 
In repeating the above for D* — > Dn, that is setting A4 = 0, one arrives at: 

^ oc Uion5oo|^((l + cos^^) +sin2^cos(250))sin^^4 (A.29) 
ail 

therefore one sees that this time the D, in D* — > Dn, will have a sin^ 9^ distri- 
bution in the rest frame of the D*. Also, it is clear that there is also azimuthal 
angular dependence present. 
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gure A. 9: Generator level scatter plot showing the expected 
e- ^ 7* ^ D*D where D* Dn. 



dependence 
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Appendix B 

Results of Angular Distributions 
of L>* in Events 

B.l e+e" ^ 7* ^ D*D* 

The method used to reconstruct D* in D*D* events is to first find D mesons that 
populate the appropriate region in Mbc, see Table 13.61 After doing so, a pion, 
either charged or neutral, is added to the D and a cut applied to AM = jM/)* — 
M^il to purify the sample. For this study only D*^ D^n^ with —>■ K^Tt^ 
and D*^ D~^tt^ with K~tt^7t^ are used. These D modes represent the 

cleanest modes available and the 7r° decay modes for D* have an efficiency that 
is fiat as a function of cos^, which is not the case of the charged pion because of 
the opening angle in the D*^ decays. 
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This method was first apphed to Monte-Carlo and in addition to looking at D* 
in D*D* we also looked at D in DD since its angular distribution is well-known 
and therefore used as a control. The MC results for a, where a is defined in 
the general formula 1 + a cos^ 6, as a function of center-of-mass energy is shown 
in Fig. IB. 11 where the solid black lines correspond to the MC input value. The 
data result of a as a function of center-of-mass energy is shown in Fig. IB. 21 where 
the red and blue lines correspond to the average value for K^-k^ and 

K^Ti^Ti^ respectively. In both cases the top plot presents the results of 
D in DD and the bottom plot of D* in D*D*. Lastly, 6 is defined as the angle 
between the D or D* momentum vector and the beam axis, its polar angle. 

To help in our understanding it is beneficial to translate the helicity amplitudes 
into partial wave amplitudes using the Jacob- Wick Transformation: 

LS ' 

(sd^Ab.,sz5. - Ab.|5(Az5. - Ab.))M(L,S) (B.l) 

where (Sm^, Tm^ | Jm^) correspond to Clebsh-Gordon coefficients. Either repre- 
sentation is valid, however, the Partial Wave representation is somewhat more 
intuitive since it is written in terms of orbital angular momentum and spin. For 
the case at hand, e'^e~ D*D*, there are three independent partial wave ampli- 
tudes, M(1,0), M(l,2), and M(3,2), analogous to the three helicity amplitudes. 
By using Eq. IB. II the following relations are obtained: 
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Figure B.l: a, defined in 1 + a cos^ 6, as a function of center-of-mass energy from 
MC. The solid black lines correspond to the MC input value. Top plot presents 
the results of D in DD and the bottom plot of D* in D*D*. 
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Figure B.2: a, defined in 1 + a cos^ 6', as a function of center-of-mass energy from 
Data. The red and blue lines correspond to the average value for —>■ K~7r~^ 
and K^tt^tt^ respectively. Top plot presents the results of D in DD and 

the bottom plot of D* in D*D*. 
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Aoo = -\I\m{1. 0) - ^ ^M(l, 2) + ^^M(3, 2) (B.2) 

All = ^1^(1, 0) - \[^M{1, 2) + y^M(3, 2) (B.4) 

Since in Eq. IA.6I the amplitudes are always in terms of the modulus square 
Eqs. IB.2tlB.4l are more useful if written as follows: 

|AooP = V(l,0)P + ^|M(l,2)p + ^|M(3,2)|2 + 
3 15 5 

4W-^|M(l,2)||M(l,O)|cos0M,.M,o- 
V 45 

2^|M(l,O)||M(3,2)|cos0MioM3. - 

^ y||M(l, 2) 1 1 M(3, 2) I cos </)M,,M3. (B.5) 

|AioP = ^|M(l,2)p + l|M(3,2)p + 

^y||M(l,2)||M(3,2)|cos0M,,M3, (B.6) 



|Aiir = -|M(l,0)p + -|M(l,2)p + -|M(3,2)p 

-2J^|M(l,2)||M(l,O)|cos0Mi,A^io 
V 45 

+2y^|M(l,O)||M(3,2)|cos0MioM3. 
-^y^|M(l,2)||M(3,2)|cos0M,.M3. (B.7) 

where ^m^^Mio is the arg(M(l, 2)M*(1, 0)), 0MioM32 is the arg(M(l, 0)M*(3, 2)), 

and 0Mi2A/32 is the arg(M(l, 2)M*(3, 2)). By applying the above equations to Eq. 
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IA.6l one arrives at: 



^|M(l,0)|^(l-z) + -|M(l,2)|^(l--^) + 



dco^e 2' ' ' ' ^ 20' ' ' " ' 7 
^|M(2,3)r(l - 1^) - ^y||M(l, 2)1 |M(3, 2)1(3^ - l)cos0M,.M3. (B.8) 

where z = cos{6y. Now the L = 3 F-wave amphtude contains two extra powers 
of the center-of-mass momentum, relative to a P-wave, and therefore should be 
small near the threshold [73]. By setting |M(3, 2)| = in Eq. MM one arrives at 
the following: 

^|M(1, 0)1^(1 ~z) + ^|M(l,2)r(l - Iz). (B.9) 



dcos^ 2' ' ' ' ^ 20' ' ' " ' 7 
Notice that this distribution is independent of any phase. By setting |M(1, 2)^ = 
|M(1,0)|2 we get: 

ocl - — cos^ = 1 - 0.65 cos^ 0. (B.IO) 



dcos9 17 

It is seen that the translation to the partial wave basis is only beneficial for 
the case near threshold and adds another unknown to our equation at all other 
energies. Therefore, we will stick to the helicity basis and will rewrite Eq. IA.6I as 
follows: 

ocl + acos^^ (B.ll) 



dcos9 
where a is defined as: 



2|Aoi| 


P-2|An| 


2 _ 


1^0 


|2 1-26 


2|Aoi 


\' + 2\An\ 


|2 + 


l^ool 


|2 1 + 26 
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with b = 4|Aoip • '^^^ hest fit to the D* production distribution is shown in 

Fig.ES 
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Figure B.3: Plot of the production angle for D* in D*D* events for data at 
4170 MeV. The sohd hne shows the best fit to the data using the fit function 
1 + a cos^(6'), where 9 is the production angle. Based on the fit a = —0.62 ± 0.01. 
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B.2 e+e- ^ 7* ^ D*D\ D* Dtt 

It has been shown earher that the angular distribution for D-mesons in D* Dn 
decays is as follows: 

= \Boo\'{{\Au\' + lAoiP) + (lAoor + l^oiP - l^iir) cos^ O'}- (B.13) 

a cos 6* 

where 6' is defined as the angle between the D in the rest frame of the D* and 
the momentum vector of the D* in the lab and Bqq is the helicity amplitude for 
D* ^ Dn. By substituting Eqs. [RSllRTl one arrives at: 

i|M(l,0)|2 + ^|M(l,2)|2(^ + l^) + l|M(3,2)p(l + 2;2) 



cicos^' 3' ' ' " 60' ^ ' ^' ^ 2r 5' 

+2y||M(l,2)||M(l,O)|cos0Mi.M,o(^-^) 

+2y^|M(l,O)||M(3,2)|cos0M,„M3.(l-3;2;) 
+ ^yi|M(l,2)||M(3,2)|cos0M,.M3.(l-32;) (B.14) 

Now assuming that the D-mesons are produced close to threshold, that is L = 3 
F-wave amplitude is small, we can set |M(3, 2)| = in Eq. IB. 141 and arrive at the 
following: 



dN 13, , 1, , „.,^ „ / 1 



^|M(l,2)p + -|M(l,0)|2-2W^|M(l,2)||M(l,0)|cos</.02 



dcos^' 60' ' " 3' ' " V 45 



21 / 1 

+ (-|M(l,2)|2 + 6W-|M(l,2)||M(l,O)|cos0o2)cos2^^' (B.15) 
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Notice that this distribution is NOT independent of (pUioMw- setting |M(1, 2) p - 
|M(l,0)p and (pMioMii = f we get: 

, ocl + — cos^ e' = 1 + 0.64 cos^ e' (B.16) 
a cos 33 

and with 0MioMi2 = we get: 

a cos 33 _ o / J_ 

60 Y 45 

showing that, even in the case where |M(3,2)| = 0, a wide range of angular 
distributions are allowed. 

Again, we see that transforming to the partial wave basis is only useful near 
threshold, therefore we will stick to the helicity basis and will rewrite Eq. IB. 131 
as follows: 

7 ocl + a cos^^ (B.18) 



d cos 6 

where a is defined as: 



I /I |2 I I /I |2 I /I |2 1 _|_ _ Si 

' - I-^OOI +1-^01 1 -1-^11 1 _ ^+ 4b 2 . 

14 |2 I 14 12 ~ a j_ 1+a yD.l^) 
l^lll + I^Oll 2 + IT 

with a = ^^^"|2 . Using Eq. IA.21I convoluted with a Crystal Ball resolution func- 
tion, the width of the resolution function was chosen to be ~ 8.5 MeV which is 
based on MC DD events generated with a momentum of 500 MeV/c, one can 
obtain the coefficient a' . A Crystal Ball function Eq. IB. 201 is used in an attempt 
to account for the ISR effects which are present. The Crystal Ball function is 
defined as follows: 
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f{x; a, n,x,a) = N< 



[B.20) 



e , for^^ > — a 

A-{B- ^)-", for^ <- a 

where A = (^)" ■ e 2 , _B = — |a| and N, a, n, x, a are the fit parameters 
which are determined by the generated MC and then fixed in the convolution. 
The fit to the D momentum spectrum is shown in Fig. IB.4[ 




Figure B.4: The D momentum spectrum. The points with error bars is the data 
collected at 4170 MeV and the solid line is the best fit result. The fit function 
consisted of Eq. IA.21I convoluted with a gaussian resolution function. The result 
of the fit yields a = 0.80 ± 0.06. 
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Using the definitions of a and a along witli tlie normalization condition that 
l^ooP + 4|^ioP + 2 1 All I = 1, one is able to arrive at the following: 

l^ool^ = 0.314 ±0.009 (B.21) 

l^iol^ = 0.080 ±0.003 (B.22) 

l^iip = 0.183 ±0.003 (B.23) 

Since in the equations Eq. IA.6I and IB. 131 the helicity amplitudes enter as the 
modulus squared, one needs not to worry about the relative phase or the sign of 
the amplitudes. Therefore, in order to generate the observed angular distributions 
for D*D* events using EVTGEN's Helicity Amplitude Model (HELAMP) one needs 
the following: 

D*0 anti-D*0 HELAMP 0.43 0.0 0.28 0.0 0.28 0.0 0.56 0.0 0.28 0.0 0.28 

0.0 0.43 0.0; 

Comparison between the data and a generated MC sample, using the above 
HELAMP model, is shown in Figs. IB. 51 and IB. 61 
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Figure B.5: Plot of the production angle for D* at 4170 MeV. Data is shown 
along with a generated MC sample using the helicity amplitudes described in the 
text. 
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Figure B.6: Plot of the angle between the D in the rest frame of the D* and the 
momentum of the D* in the lab for D*D* events at 4170 MeV. Data is shown 
along with a generated MC sample using the helicity amplitudes described in the 
text. 
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